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MOTHZV-2: A COMPUTER SIMULATION OF 
H ELIO THIS ZEA AND HELIOTHIS VIRESCENS 
POPULATION DYNAMICS 
Users Manual 

By A. W. I-Iavtstack, Jr., J , A, Witz, J. P, Hollingsworth, R. L. Ridgway, and J. D. Lopez 1 

ABSTRACT 

Consisting of a main program and 16 subroutines, MOTHZV-2 simu- 
lates the population dynamics of Heliothis zea (Boddie) and Ileliothis 
viresoens (F.) from egg- to adult. The input includes measurable and re- 
ported factors influencing population development, such as crop phenology, 
temperature, moonlight, migration, predation, parasitism, oviposition, and 
cannibalism. The model accurately predicted field population of Heliothis 
eggs, larvae, and adults. The complete computer program and example 
input and output data are given, KEY WORDS: computer-simulated Heliothis 
populations, Heliothis vivescens, Heliothis zee, insect population, dynamics, 
insect population forecasting, MOTHZV-2. 


INTRODUCTION 

MOTHZV-2 is a computer simulation of the 
population dynamics of Heliothis zea (Boddie) 
and Heliothis viresoens (F.). It is based on 
MOTHZ-1 (5) and MOTI-IZ-2 2 but incor- 
porates considerably more detailed population 
and physiological algorithms. Stinner et al. (22) 
described a somewhat similar Heliothis model; 
however, they did not include the details of 
various biological phenomena affecting popu- 
lation dynamics. Most of the basic factors, e.g., 
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mortality, parasitism, and crop phenology, were 
entered in the context of the input data. Their 
model did, however, include a spatial grid of 
crop types so that movement of the population 
could be studied. 

MOTHZV-2 consists of a main program and 
16 subroutines. Simulation begins with input 
of either eggs or moths and can be carried 
through as many generations as occur during 
one season. No provisions are included for simu- 
lating moths entering or emerging from dia- 
pause. A bookkeeping system is used so that the 
numbers of eggs, first- to third-instar larvae, 
fourth- and fifth-instar larvae, pupae, preovi- 
positing adults, and ovipositing adults can be 
recorded for each day of simulation. Printouts 
of all input as well as other calculated param- 
eters are made along with daily population 
of simulated insects. Graphs of life tables are 
optional. 

We describe the model in the order of the 
insects’ life cycle events. The model can be used 
to simulate either H. zea or H. viresoens or both 


1 



simultaneously. If in the following discussion 
no reference to species is made, both species are 
treated alike. Many biological factors other than 
those in the model affect Heliothis populations, 
but at this time these parameters are assumed 
to be constants or are neglected completely. 
Also, in some cases, more information is needed 
on the assumptions made in the model to en- 
hance accurate prediction. 

The program is written in FORTRAN and 
has been run on IBM 370 and AMDAHL 470 
computers. It is stored on disk in compiled form 
at Texas A&M University Computer Center. 

DESCRIPTION OF MODEL 

Rate of Egg Production 

Temperature has been shown to influence 
the number of eggs a moth lays, but research 
results have not been consistent (2, S, 17). 
In general, maximum egg production is attained 
between 22.2° and 25.5° C and is reduced to 
nearly zero at 12.8° and 35° C. Figure 1 (G9) 
shows the relation between oviposition rate 
and temperature as described by 

X — 175— EGLYTM | (1) 

and P^l.l— 0.05146X, (2) 

where EGLYTM=average temperature (°F) 
during the 3-hour period 
after sunset (fig. 1:G7), 
and P l =probability of oviposition 

(assumed to be 1.0 between 
22.2° and 25.5° C) . 

Quaintance and Brues (21) and Phillips and 
Whitcomb (IS) found that most eggs were laid 
in early evening. However, no data are avail- 
able to predict what happens if temperature 
inhibits oviposition early in the evening but be- 
comes ideal later or during the day. McColloch 
(13) stated that oviposition normally occurs at 
night but that it may occur during the daytime 
on cloudy or cool days. Accordingly, the 3 hours 
after sunset were assumed to be the oviposition 
period. 

Moth age has been shown to affect fecundity 
(2, 8, 17, 21). Isely (8) and Quaintance and 
Brues (21) also found an interaction between 
age and temperature in which longevity and 
reproductive potential were greater at cooler 


temperatures (21.6°-26.7° C). Higher tem- 
peratures increased daily production of eggs but 
shortened the oviposition period, thereby re- 
ducing the overall reproductive potential. Fig- 
ure 2 shows the oviposition period of Heliothis 
moths calculated for two temperatures. A func- 
tion such as the gamma might fit these data; 
however, we used a linear approximation by 
fitting straight-line segments through the 
points determined by 3 


TRUN { = / ^ TEMAVG,)/AGE, (3) 
\j — i— AGE J 

PREOVI= (5,330) (TEMAVGf 1 ' 82 ) +0.6, (4) 


IPEAK= (3,030,000) (TRUNf 3 - 11 ), (5) 

ITOP= (240,000) (TRUNr 2 01 ) , (6) 

and IEND= (911,000) (TRUNr 2 * 3 ), (7) 

where TRUN t =average temperature during 
moth’s adult life on day i= 
AGE, 

TEMAVGi—average temperature for ith 
day of moth’s life (°F), 

AGE = age of moth on day of simula- 
tion, 

PREOVI=length of preoviposition pe- 
riod (rounded to nearest 
whole day), 

IPEAK=age (rounded to nearest 
whole day) when Heliothis 
moth attains maximum re- 
productive potential, 
ITOP=number of days (rounded to 
the nearest whole day) a 
moth remains at maximum 
reproductive potential, 

and IEND=age at which a moth no longer 
oviposits. 

When the above events are calculated, a 
linear relationship between them is assumed, 
so egg production potential (P 2 ) can be readily 
calculated. 


3 Some of the terms used in the equations in the 
model description section of this paper may be slightly 
different in spelling or meaning from the terms used 
in the computer program. This was done to simplify 
explanation of what the simulation model accomplishes 
rather than how it is accomplished. 
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Figure 1.— Flow diagram of MOTHZV-2, 
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Figure 2. — The interaction of age and average daily 
temperature on reproduction potential of Heliothis 
moths. 

Researchers have shown that phenological 
phenomena (e.g., silking date) have a signifi- 
cant influence on the probability of Heliothis 
moths' being attracted to a crop and ovipositing 
(3, 6, 8, 13, 19). This influence makes a crop 
submodel necessary. A corn model developed by 
Curry (1 ) was adapted to generate the needed 
corn phenology (fig. 1: E — 1-8, 7, 8). (The 
probable attraction of corn, P 3 , to H. virescens 
is assumed to be zero, since corn is not a report- 
ed host of this insect.) A simple sorghum model 
and a portion of SIMCOT II ( li )), a cotton 
model, were also adapted (fig. 1: E — 1, 2, 4-8). 

Nemec (15, 16) reported reduced oviposition 
by Heliothis moths during periods of full moon 
and during simulated moonlight. Oviposition 
was inhibited 3 or 4 days before full moon but 
was resumed 1 or 2 days after full moon and 
usually peaked 3 or 4 days after full moon. 
Hartstack’s data (6) showed a significant re- 
duction in rates of Heliothis population increase 
during moonlight periods. Consequently, if one 
assumes that maximum egg laying is from 2030 
to 2330 hours (approximately, the 3-hour period 
after sunset) , the relative amount of moonlight 
during that period can be calculated (fig. 3) . 
Since the moon rises about 1 hour later each 
night, the 2030-2330 period has considerably 
more moonlight before full moon than after. 
The probability of maximum oviposition (P.,) 
in relation to moonlight is stored as a data 
array and used in the simulation as needed. The 
probability is increased (not >1.0) for cloud 
cover (e.g., if cloud cover is 50%, the effect of 
moonlight is reduced by 50 % ) . 


The egg-laying probability (EGLYPR) for a 
particular day is calculated by multiplying the 
four probability factors together: 

EGLYPR= (P.) (P = ) (P 3 ) (P 4 ) . (8) 

The total number of eggs laid per moth per 
day (TOLEGG) is calculated by 

TOLEGG= (EGLYPR) (EGPEMT) , (9) 

where EGPEMT=maximum number of eggs 
a Heliothis moth could lay 
per day if all factors were 
ideal. 

EGPEMT is assumed to be 300 and 400, re- 
spectively, for H. zea and II. virescens and 
would result in maximum lifetime laying of 
2,000 to 3,000 eggs per moth (8, 21). 

Egg Mortality 

There are many probable causes of egg in- 
fertility (e.g., temperature, population density, 
number of matings, and age) . MOTIIZV-2 uses 
a constant fertility rate of 90% (fig. 1: J). 
The model does not remove infertile eggs from 
the egg population until they should have 
hatched. 

Parasites, specifically Trichogramma spp., 
are an important cause of egg mortality (11, 
13). Detailed quantitative models for the popu- 
lation dynamics of various egg parasites have 



Figure 3. — The probability of maximum oviposition as 
related to the relative intensity and phase of moon- 
light during the 3-hour period after sunset. 
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not been constructed at this time; therefore, 
estimates must be used (fig-. 1: I — 2). Four 
options are available: 

Option 1 . — Enter number of parasites. Knip- 
ling and McGuire (10) presented an empirical 
population model of parasites based on the fol- 
lowing assumptions: (a) 5,000 Trichogramma 
per acre would find 50% of the host eggs, 
(b) 10,000 Trichogmmmu per acre would 
search the same plant area twice. However, 
since searching is random, the additional 5,000 
parasites would search one-half the area pre- 
viously searched and one-half the unsearched 
area, thus leaving one-fourth the area un- 
searched. (c) The searching capability of 
Trichogramma and its efficiency in parasitizing 
host eggs are assumed to be completely inde- 
pendent of host egg density. These assumptions 
were based on a search area factor (SAREA) 
of 1.0, which is assumed to occur when the 
crop becomes attractive (6- to 12-inch height 
for corn) and to increase to about 6.0 at crop 
maturity. Equation 10 is Knipling and Mc- 
Guire’s model. 

PTRICII=1 — exp ( i—Q- 693) (TRICI-lA ( 10 \ 
^ (N) (SAREA) J K ’ 

where PTRICIi=probability of parasitism by 

Trichogramma, 

TBICH— number of parasites, 

N— number of parasites re- 
quired to obtain 50% para- 
sitism, 

and S ARE A= relative search area. 

N is assumed to be 5,000 (as postulated by 
Knipling and McGuire) for native egg para- 
sites. N was calculated to be 45,000 for reared 
Trichogramma (specifically, for those reared 
at Agricultural Research Service’s Cotton In- 
sects Research Laboratory, College Station, 
Tex.), The number of parasites must be read 
in for each clay of simulation. Figure 4 shows 
typical parasitism curves calculated by the 
model. 

Option 2 . — Enter daily percentage of para- 
sitism. If the number of parasites is not 
available, an estimated daily percentage of 
parasitism can be read in. 

Option 3 , — Enter maximum daily rate of 
parasitism. The maximum daily rate of para- 
sitism is assumed to occur in com during the 
peak silking period (90%-95% silks present), 


when many eggs are present. From the egg 
parasitism percentage the model generates an 
egg parasitism curve starting from the day 
when the crop becomes attractive; figure 5 
shows the curve for corn. Parasitized eggs are 
not removed from the egg population by the 
model until after the parasites would have 
emerged, since field counts of eggs may include 
parasitized eggs. A running total of parasitized 
eggs is kept so that the apparent field para- 
sitism can be calculated for each day for com- 
parison with field data for verification. 

Option U . — Enter number of Trichogramma. 
The rate of parasitism by Trichogramma is cal- 
culated as a function of the number and the 
age of adult Trichogramma and of environ- 



Figure 4. — Typical eg'g' parasitism curves generated 
by model. 



TIM]: (DAYS) 

Figure 5. — Parasitism or predation curve generated 
for corn (option 3). 
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mental temperature. The potential population 
of Tricho gramma emerging from parasitized 
Heliothis eggs is calculated and can be used .to 
predict parasitism rates once the population of 
Trichogramma has been initialized. 

Egg predators are handled like parasites. 
Knipling and McGuire’s (10) model, equation 
10, was adapted for egg predators. The value of 
N is assumed to be 14,500 and 45,000 for native 
egg predators and laboratory-reared Chrysopa, 
respectively (determined from field data by 
the Cotton Insects Research Laboratory). Field 
counts of predators are obtainable and can be 
used in the model. Egg predators have been 
placed by the Cotton Insects Research Lab- 
oratory into seven groups. The number of pred- 
ators in each group can be adjusted by a relative 
efficiency factor, which can be estimated from 
intuitive and factual knowledge of the preda- 
tors. The predator groups are Scymnus, other 
Coccinellidae, other Coleoptera, Orius, Geocoris, 
other Hemiptera, and Chrysopa immatures. 

The daily percentage of predation can be used 
if the number of predators is not available. 
Also, as described for parasites, the maximum 
daily rate of predation can be used, and the 
model will generate daily rates. 

The model provides for a daily input of mor- 
tality from ovicides for each species. A constant 
4.5% daily rate is assumed 'for mortalities 
caused by other agents. 

Rate of Egg Development 

Researchers have shown that temperature 
is the major factor regulating the rate of de- 
velopment of Heliothis under normal conditions 
(2, 6, 8, 12, 21). Development ceases at tem- 
peratures below 12.6° C and is interrupted at 
temperatures above 33.3° C (fig. 1: G — 1). 

The 484.9 degree-days required for a complete 
life cycle of Heliothis ( 6 ) was divided into 
stages (e.g., eggs and small larvae) . It was de- 
termined that 40.5 degree-days between 12.6° 
and 33.3° C were required for development and 
hatching of eggs. Degree-days above '33.3° C 
are subtracted from the sum. Figure 6 shows 
the development time in days required for var- 
ious temperatures. Hartstack (6) and Fye and 
McAda (2) reported a slightly longer time for 
H. virescens than for H. zea; therefore, the 
model is written so that a percentage of reduc- 
tion or increase can be used. 



Figure 6. — Development time as a function of tem- 
perature for Heliothis eggs, larvae, pupae, and adults, 


Larval Mortality 

Knipling ( 9 ) developed a population model 
for larval parasites that is, in principle, identi- 
cal to the Trichogramma model discussed ear- 
lier. The value of N for larval parasites was 
assumed to be 300, as suggested by Knipling. 
Parasitism is applied only to small larvae (first- 
to third-instar) . If numbers of parasites are not 
available, the daily percentage of parasitism or 
the maximum daily rate can be used for small 
larvae. The larval parasitism pattern is calcu- 
lated exactly as described for egg parasitism 
(shown in fig. 5). 

Knipling and McGuire’s (10) model, equation 
10, was also adapted for predators of small 
larvae. The value of N was assumed to be 9,000 
and 13,000 for native and laboratory-reared 
Chrysopa immatures, respectively (determined 
from field data by the Cotton Insects Research 
Laboratory ) . The daily rate of predation calcu- 
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lated by equation 10 is reduced in relation to 
the age of the larvae, as follows: 

PREDOR'= 

(PREDOE) exp (11) 

where PRED OR' = percentage of predation 
after adjustment for age 
of larvae, 

PREDOR— percentage of predation 
calculated by equation 10, 
find AGE— age of larvae in days. 

Larval predators are grouped by the Cotton 
Insects Research Laboratory into the same 
groups listed for egg predators. The larval 
predator numbers are adjusted by estimated 
relative efficiency factors. Daily percentage of 
predation of small larvae can be read in if the 
number of predators is not available. The daily 
percentage is reduced for age of larvae by 
equation 11. If the maximum daily rate of 
predation of small larvae is used, the larval 
predation pattern is calculated as was described 
for egg predation (shown in fig. 5), and again 
the daily percentage is reduced for age of larvae 
by equation 11. 

Mortality of larvae by cannibalism is con- 
sidered for II. zea only (fig. 1: M), Subroutine 
CROP calculates the number of “sites" avail- 
able daily for egg laying. Pieters and Sterling 
(20) reported that II. zea eggs are distributed in 
the field as described by the negative bionomial 
and that the dispersion factor (DF) of the 
negative bionomial for II. zea eggs in cotton is 
0.36. Cannibalism of II. zea is considered by 
many to be negligible in cotton but very high in 
corn, since II. zea eggs are concentrated on silks 
when these sites are available. A dispersion 
factor of 0.36 is assumed here for corn. The 
number of sites available is assumed to be the 
number of plants times the relative search area 
(SAREA) until silks appear. At that time the 
number of sites is reduced to the number of 
silks available, since most eggs are laid on silks 
when they are available. It is assumed that 
when silks become 16 days old they are dry and 
no longer attractive for oviposition. It is also 
assumed that a larvae reaches the fourth instar 
before it becomes cannibalistic and destroys all 
smaller larvae on the same site. Only a small 


percentage of small larvae will be destroyed in 
1 day, but this percentage will increase as the 
larvae develop and are more active. The num- 
ber of surviving larvae is determined by 

PRBLARj=TOTALS3j-7-SITES/, (12) 

ATEj=0.05-f 0.05 ( AGEj) , (13) 

and POP'j—POPj— (POP,) 

(PRBLARi) (ATEj), (14) 

where PRBLARi ““probability of a large II. 

zea larva being on a site on 
day i, 

TOTALS; —number of surviving large 
H. zea larvae on day i, 
SITES;— number of sites available 
for oviposition on day j 
(day j is the day when the 
eggs that are being ad- 
vanced through the simu- 
lation were laid), 

ATE i —percentage of small II. zea 
larvae cannibalized on sites 
with a large larvae on day 
i. 

POPi— the number of small H. zea 
larvae surviving on day i 
from eggs laid on day j, 
and POP'j— the small II. zea larvae sur- 

viving cannibalism until 
the next day. 

This calculation is done each day until the 
surviving small larvae laid on day j become 
large larvae. On the last day of the third instar 
the probability of more than one surviving small 
larva being on the same site is calculated, and 
the population is reduced to one per site by 
solving 

XBARj— POP h-SITES/, (16) 

PRBNEW,-!- (ssqlb®) " < 16 > 

and POP'j— (SITES/) (PRBNEW/), (17) 

where XBAR,— number of surviving small H, 
zea larvae per site on day i, 
PRBNEW/— probability that one or more 
small larvae are on a site on 
day j, 
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and DF— dispersion factor of the nega- 

tive bionomial (0.36) . 

The model provides for separate inputs of 
mortalities of small and large larvae from in- 
secticides. A constant daily rate of 4.0% is 
assumed for mortalities from natural causes. 

Rate of Larval Development 

Degree-day figures used for the average de- 
velopment of small larvae (one- to three-instar) 
and large larvae (four- or five-instar) are 81.7 
and 120.6, respectively, calculated as for eggs. 
Isely (<S') reported different development times 
for H. zea larvae held at constant temperatures 
when they were fed various crops and crop 
parts. Corn ears were apparently the most nu- 
tritious, since larval development was signifi- 
cantly faster with this diet. Our unpublished 
field data at College Station support these re- 
sults, and this factor has been included in the 
model. The nutrition factor used for corn plants 
and for cotton plants and fruit is 1.0. If larvae 
are developing on corn plants when ears are 
present, 0.65 is used. The nutrition factor for 
larvae developing on sorghum is 1.1. 

Pupal Mortality 

The only mortality inputs for pupae are a 
constant natural mortality of 3.0% per day. 
Mortalities from insecticides can also be used. 

Rate of Pupal Development 

Pupae require 123.9 degree-days for develop- 
ment. Figure 6 shows this temperature-time 
relationship. 

Adult Mortality 

Fye and McAda (2) and Sterling 1 reported 
a temperature survival relationship for Helio- 
this moths. Figure 7 shows the probability of 
survival of adult moths as related to age and 
temperature. The following equations were 
developed to describe this survival probability: 

FACTK=1.2919 

+ (0.02724) (TEMAVG) , (18) 

1 Winfield L. Sterling, personal communication. 
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and SURV= 1.0—1.7724 

/ —11.4457 \ /1n . 

6XP y(FACTK) (L) J ’ (19) 

where FACTK— temperature constant for ex- 
ponential equation, 

TEMAVG =average air temperature 
(°F) on day of simulation, 
SURV=probability of survival to L 
age, 

and L=age of moth in days. 

Survival of moths at temperatures below 
20° C is assumed to be the same as at 20° C. 
Daily mortalities from insecticides can be used 
for adult moths of both species. Snow et ah 
(24) reported that few moths lived beyond the 
sixth night after their emergence in the field. 
To simulate this high mortality, we used a daily 
rate of 15.0% for mortalities from natural 
causes. 

Migration 

Crop phenological data are used to calculate 
emigration from corn and sorghum. In corn 
and sorghum the mature date (IM) is used to 
start the emigration, which builds to a maxi- 
mum at the dry date (ID) and then remains 
constant for the rest of the simulation. The 
maximum percentage of emigration is read in. 
These moths are stored and can be used later 
as immigrants to other crops. 

Duration of Oviposition 
Fye and McAda (2) and Isely (8) showed 



ARR (DAYS) 

Figure 7, — Probability of Heliothis moth survival as 
related to age and temperature. 



that the duration of oviposition is a function 
of temperature. We used the following equation 
to estimate the duration and to terminate it: 

IGENZV6=28.25 

(0.24037) (TEMAVG), (20) 

where IGENZV6— length of oviposition pe- 
riod, 

and TEMAVG=average air temperature 
(°F) on day of simulation. 

DESCRIPTION OF 
SUBROUTINES 

Temperature (TEMP) 

All temperature and cloud r cover input and 
the calculation of degree-days are handled by 
TEMP (fig'. 1: G). There are two options for 
putting daily temperature in: (1) hourly tem- 
perature and (2) daily maximum and minimum 
temperatures. The hourly temperature takes 
precedence if the data are available. When 
maximum and minimum temperatures are used, 
hourly temperatures are estimated by assuming 
a consine-wave fluctuation between maximum 
and minimum readings at 1400 and 0200 hours, 
respectively. 

TEMP calculates the degree-days for corn as 
well as for Heliothis. For example, the degree- 
days for Heliothis are calculated by 

24 

ED a — ^ (T.v — 12.6) — 2 (T, v — 33.3) (21) 

X = l 

and ED 0 '=EDa-|-EDn.i+ . . . -(-EDa.,,, (22) 

where ED„=degree-days on day a, 

X=hour, 

T.*— hourly temperature on day a 
(°C), 

ED„'=aceumulated degree-days from 
day n to day a, imposing only 
positive values of (T* — 12.6) 
and (T.v — 33.3) be summed, 
and ?i=day number of start of simula- 

tion (Julian). 

Multiplying (T.r— 33.3) by 2 causes a reduc- 
tion in degree-days for temperatures above 
33.3° C, thereby simulating the retardation of 


the development of Heliothis by high tempera- 
tures. Also calculated in TEMP are the ac- 
cumulated degree-days for corn, sorghum, and 
cotton, the average daily temperature, the aver- 
age nightly temperature, the average tempera- 
ture for the 3 hours after sunset (oviposition 
temperature) , the time of sunrise and sunset, 
and the percentage of cloud cover during the 
oviposition period. 

Generation Length (GENLN) 

The degree-days calculated in TEMP are used 
to calculate the number of days required for 
development of the Heliothis (fig. 1: L). The 
following equation illustrates the GENLN func- 
tion: 

D*=ED'<i. n — ED'<i>40.5, (23) 

where Dr= degree-days required for develop- 
ment of eggs, 

ED'«B=number of degree-days accumu- 
lated to day a+n, 

ED'„=number of degree-days accumu- 
lated to day a, 

and «=riumber of days required for egg 
development from day a. 

Also calculated in GENLN are the lengths of 
the preoviposition and oviposition periods 
(equations 4 and 22), The development times 
are also adjusted for the percentaegs of differ- 
ence between species. 

Moonlight (MOON) 

This subroutine calculates the number of 
days, before or after, that each day of the 
simulation is from full moon. MOON then stores 
in the ADMOON array the proper moonlight 
adjustment factor, P.„ for egg laying. 

Insecticides (SPRAY) 

This subroutine reads in daily mortalities 
caused by insecticides for eggs, small larvae, 
large larvae, pupae, adults, parasites, and 
predators for each species. 

Parasites (PARINP) 

This subroutine calculates parasitism of eggs 
and larvae as described under “Egg Mortality.” 
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Predators (PRDINP) 

This subroutine calculates predation of eggs 
and larvae as described previously under 
“Egg Mortality.” 

Adult Moths (ADUI/T) 

This subroutine calculates adult mortality as 
described previously. 

Crop Phenology (CROP) 

This subroutine calculates the various pheno- 
logical stages, the relative search area, and the 
number of sites (e.g., ears) for corn. Degree- 
days as previously described are used to make 
these predictions. 

Cannibalism (CANBAL) 

This subroutine calculates the probability of 
cannibalism of H. zea larvae. 

Migration (MIGRAT) 

This subroutine calculates the probability of 
moths migrating away from a particular crop. 

Field Population (FLDPOP) 

This subroutine converts light-trap catches to 
field populations expressed as moths per acre 
<,5, 7 ) . The light-trap catches are adjusted fol- 
low night temperatures. Light-trap catches for 
5 years (1970-74) were reduced significantly 
before June 10. The effect of temperature was 
laige early (April) and decreased to nearly zero 
by June 10. Hardwick (4) also found a signifi- 
cant correlation between low night tempera- 
tures and light-trap catches. Figure 8 shows the 
effect of temperature on trap catch on four 
dates (Julian). We used a quadratic equation 
to calculate the adjustment for any particular 
day. The coefficients of the quadratic (A) and 
the base factor (F=26.7° C) are calculated by 
four linear equations whose variable is time 
(Julian date): 


Al— 812.4 — 5.8 (IDY) , 

(24) 

A2=28.6753+0, 20281 (IDY) , 

(25) 

A3=0.252667— 0.001762 (ID Y) , 

(26) 

F=136.9— 0.86667 (IDY) , 

(27) 


F'=A1+A2(NA)+A3(NA) 2 , (28) 

and TRAPIN=F~F', (29) 

where IDY= Julian date, 

F'= temperature factor, 
NA=average temperature from 1800 
to 0600 hours (°F), 

and TRAPIN=adjustment factor for trap 
catch. 

This adjustment is limited to the College 
Station, Tex., climate or a similar one. The 
time factor would have to be changed for areas 
where early-season temperatures differ signifi- 
cantly. These equations are only valid to clay 
153 of the simulation, and all temperatures be- 
low 12.8° C are assumed to be 12.8° C. 

Graph (GRAPH) 

This subroutine may be used if a graph of the 
output or input is needed. 

Heliothis Population Input (EGGINP) 

This subroutine reads in the Heliothis popu- 
lations for initiating the simulation. These pop- 
ulations can be either eggs or moths. White eggs 
(eggs up to 1 day old) or total eggs can be used 
as input. The eggs are read in daily. If total 
eggs is used, EGGINP divides total egg num- 



Fiourb 8.— Relative insect inactivity as a function of 
temperature at different times of the year com- 
pared to the inactivity on day 164. 
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her by 0.31 to obtain an estimated white-egg 
count and moves the count to 3 days earlier. 

Light-trap catches tabulated by sex or on a 
total basis can be used as input. Subroutine 
FLDPOP converts the catch to an estimated 
number of moths per acre. Male and female 
moth data are preferred if available because the 
sex ratio is calculated and used later by CROP 
to estimate egg-laying probability. 

Double Three-Day Moving Average 
(SMOOTH) 

This subroutine applies a double 3-day mov- 
ing average to insect population data. 

Trichogramma (TRKGRM) 

This subroutine is a model of the population 
dynamics of Trichogramma. Development of 
Trichogramma in Heliothis eggs is controlled 
by temperature, and potential parasitism is con- 
trolled by number and age of adult Trichogram- 
ma. These relationships are described by 

PTEMPi= 1 72 — TEMAVGi | (30) 

and PTEMP Z =1.56 — 0.07 (PTEMPO , (31) 

where PTEMP^the degrees average daily 
temperature is from 72° F, 
imposing that if PTEMP! 
>22, PTEMPi=22 and if 
PTEMP! <8, PTEMP=8, 
PTEMP 2 =the activity factor of adult 
Trichogramma attributable 
to tempei’ature, 

and B t = — 0.693 h-[(N) (SAREA)], (32) 

and B 2 = 1— exp [PARSmBO ] 

(PTEMP 2 ) (1— POS), (33) 

where Bi=index of potential parasitism, 
B 2 =probability of parasitism, 
POS=mortality of parasites from in- 
secticides, 

SAREA=relative search area, 
and N=number of Trichogramma adults 
required to parasitize 50 % of eggs 
in 3 days (15,000). 

The parasitized eggs are subject to death by 
predators and other natural agents and by in- 
secticides until the Trichogramma emerge. It 


is assumed that the parasites will live 10 days, 
that two female parasites will emerge from 
each egg, and that the number of eggs para- 
sitized will increase as a function of search area 
(to account for eggs missed) . The fertility rate 
of adult Trichogramma is read in as array 
FERATE. 

Function Distribution (DISTR) 

DISTR calculates the development distribu- 
tions of cohorts of eggs, larvae, and pupae. A 
normal distribution is assumed (22) . The stand- 
ard deviations (SDEV) are calculated in sub- 
routine GENLN as a function of temperature 
(2 ) . Array GAUS is the cumulative distribution 
function for eight standard deviations centered 
on the mean of the normal probability distribu- 
tion. 

VALIDATION OF THE MODEL 

Heliothis populations on corn, sorghum, and 
cotton were sampled at least two times a week 
during low populations and three or four times 
a week during peaks. The samples were taken 
by examining all plants in 6 (4 feet of row at 
six random locations in the field. 

Simulations were run from April 1 to July 
24, 1974, which covers the first three Heliothis 
generations. The first two H. zea generations 
at College Station oviposit primarily on corn 
and sorghum, and the third generation oviposits 
on cotton. It was assumed that the first two 
generations of H. virescens would oviposit pri- 
marily on wild hosts and that the third genera- 
tion would move into cotton. Since no wild hosts 
were sampled, it was assumed that all eggs 
found on corn and sorghum were H. zea, an 
assumption verified by laboratory rearing of 
collected larvae. No counts of parasites or pred- 
ators were made; therefore, option 3 was used 
in subroutines PARINP and PRDINP (fig. 5) . 
Two simulations were run and compared with 
field-collected data. 

Light-trap input . — Counts of H. zea males 
from two light traps (each having a 40-watt 
black-light lamp) located in a cornfield were 
used as model input from April 1 to May 20, 
1974. These moths were assumed to have 
emerged from diapause. The next two genera- 
tions were predicted and compared with actual 
field counts. Figure 9 shows the light-trap input 
converted to total eggs (April 1-May 20) , the 
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Figure 9, — Predicted and actual field counts of eggs. 

Input to model, April 1 to May 20, 1974. 

predicted number of 'total eggs (May 21- June 
29), and the number of eggs counted in the 
field. Actual hourly temperature was used for 
the entire simulation. The predicted June peak 
was 1 day late and the predicted number of eggs 
was low, the latter caused primarily by the 
light-trap model, probably one of the weakest 
parts of MOTHZV-2. 

Figure 10 compares the predicted third gener- 
ation of ovipositing moths migrating into cot- 
ton and the actual field count of eggs in 
cottonfields closest to the cornfield. The pre- 



Figure 10. — Predicted number of moths migrating from 
corn and actual field counts of eggs in cotton. 
Input to model, April 1 to May 20, 1974. 


dieted moth peak was 3 days ahead of field egg 
peak, as expected, since eggs would require 2 or 
3 days to hatch. 

White-egg input . — Counts of white eggs in 
the cornfield from April 1 to May 20 were used 
as input, and predictions of the next two gener- 
ations were made. Figure 9 shows the white- 
egg input and the predicted second generation. 
The predicted peak was 1 day earlier than the 
actual peak, and the predicted egg peak was 
larger than the actual field peak. (Field scouts, 
bear in mind, are less than accurate. Moreover, 
many eggs may have been removed by predators 
before the count.) The predicted third gener- 
ation moth peak was also 3 days before the 
actual egg peak in cotton. 

Figures 11 and 12 show predicted and actual 
counts of small and large larvae. The dates of 
these peaks are very close; however, the magni- 
tudes of the peaks vary similarly to the egg 
peaks, as discussed. Further validation of as- 
sumptions concerning development periods and 
mortalities could alleviate differences in mag- 
nitude. It is felt that the field measurements 
were in error, even though very competent 
entomological help was used, because it is im- 
possible to find every egg or larva and because 
field counts were not made every day. A factor 
that may affect the size of the small and large 
larval peaks is the subjective judgment of which 
group to put them in. The model apparently is 



Figure 11. — Predicted and actual field counts of small 
larvae. Input to model, April 1 to May 20, 1974. 
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Figure 12,— Predicted and actual field counts of large 
larvae. Input to model, April 1 to May 20, 1974. 

dividing them at an earlier age (1 day off) than 
does the field scout. 

Similar results were obtained when data from 
previous years at College Station and Pearsall, 
Tex., were used. Simulation of H . zea in sor- 
ghum was also accomplished with similar 
results. 


DEFINITION OF TERMS 


ADMOON 

ADTS 


ADTSMG 

AGE 

AP 

APARA 

APOSON 

B 


BOLLS 

C 

CLOUD 

DATX 

DEATH 

DEATIIY 

DEGDAY 


Main Program 

Moonlight adjustment calculated for each 
day. 

Number of adults migrating on each day 
of simulation plus the number that mi- 
grated and survived previously. 

Average number of ovipositing adults that 
have migrated and survived. 

Number of days moth has been ovipositing. 
Proportion of the eggs not hatched. 
Storage of number of eggs parasitized 
during simulation. 

Storage of number of eggs sterilized from 
insecticide during simulation. 

Used in calculation of probability of 
Trichogramma parasitism (option 4, 
KODE 5). 

Calculated cotton bolls per acre. 

Mortality probabilities. 

Percentage of cloud cover read in for each 
hour. 

Array of read-in data to be plotted with 
output. 

Natural mortalities of all stages of insects. 
Mortality of adults the previous day. 
Array of degree -days required for each 
stage of Heliothis development. 


EARS 

ED 

EDCORN 

EGGFRT 

EGLYCR 

EGLYPR 

EGLYTM 

EGP 

ELPDTA 

EQELAD 

FACT A A, 
FACTBB 
FACTD 

FM 

FOODFK 

HEADS 

IDX 

XEND 

IFIRST 

IFISTT 

IFRUIT 

IGEN 

IGENV 

IGENZ 

IGKODE 

IH 

I LAST 

ILASTC 

ILASTT 

IPEAK 

IPX 

ITOP 

IV 

IXBAR 

JFIRST, 

JLASTC 


Calculated corn ears per acre. 

Array of degree- days calculated for Hello - 
this development and cotton growth. 

Army of degree -days calculated for corn 
and sorghum development. 

Percentage of Heliothis eggs that are 
fertile. 

Probability that a female will oviposit be- 
cause of a crop factor. 

Probability that a female will oviposit be- 
cause of crop and temperature factors. 
Average temperature (°F) during 3-hour 
period of maximum ovi position. 

Maximum number of eggs a Heliothis 
moth will lay in 1 day. 

Probability that a moth is at maximum 
egg production. 

Relative number of female moths consid- 
ered to be ovipositing at maximum rate. 
Factors for averaging migration of moths 
from fields of the same crop. 

Calculated light-trap factor to correct trap 
catch to moths per acre. 

Array of moonlight adjustment factors, 
Nutritional factor affecting development 
time of Heliothis, 

Calculated heads of sorghum per acre, 

Used for plotting DATX. 

Number of days a moth will produce eggs. 
Day number of first day of simulation 
(Julian) . 

First day of temperature data to be read 
in and degree-days calculated (Julian), 
Day number first silk appears in corn 
(Julian) . 

Array containing calculated average de- 
velopment times for all stages of FI. zea 
and FI. virescens , 

Array containing calculated average de- 
velopment times for all stages of IL vires- 
cens. 

Array containing calculated average de- 
velopment times for all stages of H . zea . 
Plot code. 

Code for FI. zea. 

Day number of last clay of simulation 
(Julian) . 

Day number of last day of input to model 
(Julian) , 

Last day of temperature input and degree- 
days calculated (Julian). 

Number of days it takes an ovipositing 
moth to reach maximum egg* production. 
Code for double 3-day moving average of 
migrating moths, 

Number of days moth will be at maximum 
egg production. 

Code for FI. virescens. 

Mean length of pupal development time in 
days. 

Variables to flag feedback of moths. 
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K.KODE 

KODE 

K.RUN 

KTRAP 

LLPOP 

LX 

NITMIN 

NN 

NPLOT 

PAKSIT 

PCTCRP 

PCTMG 

PDERTE 

PDLRTE 

POAPOP 

POP 

POS 

PP 

PPOP 

FRED 

PREDOR 

PRMGRT 

PTEMP 

S 

SAREA 

SDEV 


SEXRTO 

SITES 

SLPOP 

SQUARE 

SRVIVE 

SS 

T 

TCLOUD 


Index of type of input to model (eggs or 
moths) , 

Code read in that flags the option wanted. 
Code of a particular simulation. 

Index of type of input to model (eggs or 
moths) . 

Temporary storage of field population of 
1 -day-old large larvae. 

Mean length of development time for small 
larvae in days. 

Average nighttime temperature (°F). 
Variable used to flag day that simulation 
is at I FRUIT or past. 

Plot code. 

Calculated percentage of parasitism of 
eggs and larvae for each day of simula- 
tion. 

Percentage of the area in corn, sorghum, 
cotton, pasture, and other crops. 

Maximum percentage of the adult popula- 
tion that would migrate out of the crop on 
any one day. 

Predator rating for egg predators. 
Predator rating for larval predators. 
Temporary storage of field population of 
1-day-old ovipositing adults. 

Temporary storage of Hcliothis eggs. 
Read-in array containing daily mortality 
of Hetiothis, predators, and parasites at- 
tributable to insecticide treatments. 
Temporary storage of field population of 
1-day-old larvae. 

Temporary storage of field population of 
1 -day-old pupae. 

Number of predators per acre of the seven 
groups of predators described in text. 
Calculated daily percentage of predation 
of eggs and larvae for each day of simula- 
tion. 

Calculated percentage of adult population 
that will migrate on any one day. 

Used in calculation of probability of 
Tricho gramma parasitism (option 4, 
KODE 6). 

Scale factor 1 for plots. 

Relative plant area searched by Tricho- 
gramma (arbitrary scale). 

Array of calculated standard deviations of 
development time for eggs, larvae, and 
pupae attributable to temperature. 

Sex ratio of adults (males + 1) / (females 

H-l), 

Number of sites available for oviposition. 
Temporary storage of field population of 
1-day-old small larvae. 

Calculated cotton squares per acre. 
Number of fertile eggs remaining in popu- 
lation after all mortalities. 

Scale factor 2 for plots. 

Total survival probability of small larvae. 
Percentage of cloud cover during time of 
maximum oviposition. 


TEMAVG 

TINF 
TIME A 

TOTAL 

TOT MGS 

TOTO 

TP 

THUN 

TT 

TPP 

XBAR 

XSD 

FACTK 


ATE 


ATEUP 

CANLAR 

POP 

PRBLAR 

PRBNEW 

S 

XBAR 

XPOP 


A20, 
A21, 
A22 
A30, 
A31, 
A32, 
A33 
A40, 
A41, 
A42 
AM20, 
AM21, 
AM22 
AM30, 
AM31, 
AM32, 
A M3 3 


Average air temperature on day of simula- 
tion (°F). 

Total infertile eggs. 

Array of expected or measured high and 
low temperatures. 

Array for storing total Heliothis popula- 
tions of various stages. 

Total number of ovipositing adults that 
migrate and survive. 

Percentage of eggs parasitized during sim- 
ulation. 

Number of eggs parasitized on any one 
day. 

Average temperature during oviposition 
period (°F). 

Total probability of something surviving. 
Number of eggs sterilized on any one day, 
Day number that the mean number of eggs 
hatch (Julian). 

Standard deviation of egg hatch. 

ADULT 

Temperature effect on adult longevity. 

CANBAL 

Percentage of small II. zea larvae canni- 
balized on sites with a large larvae on a 
given day. 

Total number of small II. zea larvae can- 
nibalized. 

Number of small II. zea larvae canniba- 
lized by other small larvae. 

Number of surviving small II. zea larvae. 
Probability of a large II. zea larva being 
on a site on a given day. 

Probability of a small H. zea larvae being 
on a site on a given day. 

Number of sites available for II. zea 
larvae. 

Average number of small larvae per site, 
Total number of small larvae cannibalized 
from this cohort. 

CROP 

Coefficients for calculating egg-laying 
probability on corn during silking period. 

Coefficients for calculating egg-laying 
probability on corn before first silk ap- 
pears. 

Coefficients for calculating egg-laying 
probability on corn after silking period. 

Coefficients for calculating egg-laying 
probability on sorghum after boot stage. 

Coefficients for calculating egg-laying 
probability on sorghum before boot stage. 
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DD 

Degree- days required for corn or sorghum 
to develop from IA to peak silk or head. 

F 

DDA 

Degree-days from corn and sorghum emer- 
gence to IA. 

FD 

DDE 

Degree- days from first square to first 
boll in cotton. 

FI 

DDF 

Degree-days from I A to IF in corn and 

IN 


sorghum. 

IPX 

DDM 

Degree-days from first boll to maturity in 



cotton. 

R 

DELTAX 

Percentage of crop development per day. 


D1STS 

Percentage of eggs hatched. 


IA 

Day number when crop becomes attractive 



to Helio this (Julian). 

Al, 

IE 

Day number when crop emerges (Julian), 

A2, 

IF 

Day number when first fruit appears 

A3 


(Julian) . 

EDICK 

IM 

Day number when crop becomes fully ma- 



ture (Julian). 

IGENX 

IP 

Day number when crop planted (Julian). 


ICCROP 

Set to value of KODE 11 (crop). 

IKODE 

NATRAC 

Day number when crop becomes attractive 
to II el to t his (Julian). 


PLANTS 

Number of plants per acre. 

IN 

S 

Factor used to calculate normal distribu- 
tion for development of fruit. 

IX 

TP 

Factor used to calculate normal distribu- 
tion for development of fruit. 

P 

YIELD 

Estimated yield of cotton bales per acre. 

X 

YIELDB 

Calculated peak number of bolls. 


YIELDS 

Calculated peak number of squares. 



Base factor used in temperature adjust- 
ment of light-trap catches (2G.7° C), 
Average distance (feet) a moth can fly in 
1 night. 

Factor used in temperature adjustment of 
light-trap catches. 

Insect species code. 

Code for locating light-trap data for 
double 3-day moving average. 

Effective radius (feet) of light trap. 

GENLN 

Coefficients for calculation of standard 
deviations in H cl lot his development times. 

Accumulated degree-days at start of Helio - 
this development. 

Used in calculation of Helio this develop- 
ment times. 

Flags generation length calculations so 
that calculation of multiple simulations of 
same weather data will not be repeated, 
Insect species code. 

Used in calculation of Helio this develop- 
ment time, 

Percentage of difference in development 
time of IL zea and H. virescens. 

Standard deviations for 1 day. 

GRAPH 


DISTR 

GAUS, Arrays of accumulated normal distribu- 
GAUS1, tions used to vary development time of 
GAUS2 rieliotkis. 

EGGINP 

B Number of H. zea eggs in field on first day 

of simulation. 

BFMALE Number of II. zea females. 

BMALE Number of H. zea males, 

C Number of TL vircscens eggs. 

CFMALE Number of II. virescens females. 

CM ALE Number of II. vircscens males. 

IDD Used for calculation of missing egg data, 

L Set to code of type of input (eggs or 

moths) * 

X Adjustment factor for total egg input. 

XDL1, Used for calculation of missing egg data, 
XDL2 

XEGG Used for adjustment of nutritional factor 
for sorghum, 


B 

Used for calculation of variable format 

DIGIT, 

DIGIT1, 

DIGIT2, 

Arrays used for plotting, 

INSECT 

Array used for plot headings and identifi- 
cation. 

ISTAGE 

Array used for plot headings and identifi- 
cation. 

ISTG 

Used to calculate variable scale factors 
and printout. 

K 

Used for calculation of variable format. 

KSTG 

Used to calculate variable scale factors 
and printout. 

NSCALE 

Used to calculate variable scale factors 
and printout. 

SCALE 

Used to calculate variable scale factors 
and printout. 

sxv 

Code for II, vires ecus. 

sxz 

Code for IL zea. 

sxzv 

Code for H. zea and II. virescens , 

SX1 

Used to calculate variable scale factors and 
printout. 

YFA 

Variable format used to plot data. 

VFB 

Variable format used to plot data. 


FLDPOP 


MIGRAT 


Al, 

Coefficients used in temperature adjust- 

XD 

A2, 

ment of light-trap catches. 


A3 


XM 

E 

Efficiency of light trap. 



Day number crop becomes completely un- 
attractive to Ilcli othis (Julian), 

Day number crop becomes mature (Jul- 
ian), 
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A 

E 

J 


MOON sp 

Set to day number of first full moon SPID 
(Julian). X 

Day number of simulation (Julian). Y 

Phase-of-inoon code. 2 


Number of spiders per acre (specifically 
for data collected at Pecos, Tex.). 

Number of spiders per acre. 

Pj*obabiiity of predation of eggs. 
Probability of predation of small larvae. 
Probability of predation of large larvae. 


PARINP 


SMOOTH 


B, 

FX 

L 

N 

OLARPA 

OTHERS 

PA 

PEND 
PDOWN 
PM A XI 
PM I NI 
POTHRS' 

PTRICH 

TRICH 

X 

Y 


ALB 

CALLOP 

CHA 

CHL 

CHRYSP 

HEMP 

IFF, 

ILL 

NE 

NL 

OTHEGG 

OTHLAR 

PCHRYS 

POTHRS 


Used to calculate percentage of parasitism. 

Set to parasitism code (option). 

Number of parasites per acre required to 
obtain 50 Yr parasitism. 

Number of naturally occur ting parasites 
of He} io this larvae per acre. 

Number of naturally occurring parasites 
per acre. 

Percentage of maximum parasitism of 
eggs and larvae. 

Last day of parasites in crop (Julian). 
Start of downtrend of parasites (Julian). 
First peak of parasites (Julian). 

Start of buildup of parasites (Julian). 
Probability of parasitism by naturally oc- 
curring parasites. 

Probability of parasitism by Tricho- 
y ram mu. 

Number of released Trichoyramma per 
acre. 

Probability of egg parasitism, 

Probability of larval parasitism. 

PRDINP 

Number of lady beetles per acre (specifi- 
cally for data collected at Pecos, Tex.). 
Number of callops per aero (specifically 
for data collected at Pecos, Tex,). 

Number of Chrysopa adults per acre 
(specifically for data collected at Pecos, 
Tex,), 

Number of Chrysopa larvae per acre 
(specifically for data collected at Pecos, 
Tex,) . 

Number of released Chrysopa per acre. 
Number of Hemiptera per acre (specifi- 
cally for data collected at Pecos, Tex.) 
Used for calculating number of predators 
when actual counts are not available. 
Number of egg predators needed per acre 
to obtain 50<7< predation when search 
area — 1.0, 

Number of larval predators needed per 
acre to obtain 50 ( A predation when search 
area = 1.0. 

Number of naturally occurring egg preda- 
tors per acre. 

Number of naturally occurring larval 
predators per acre. 

Probability of predation by released Chry- 
sopa. 

Probability of predation by naturally oc- 
curring predators. 


IP 

IXF 

S 

SM 


Factor used in double 3-day moving aver- 
age. 

First day of smoothed data. 

Factor used in double 3- day moving aver- 
age. 

Equivalent array used to calculate double 
3-day moving average. 


SPRAY 

P Percentage of mortality of Ileliothis stages, 

predators, and parasites on any one day. 

TEMP 

DAYLN Number of hours of daylight on each day 
of simulation. 

ET Minimum temperature ( fl F) for fleliothis 

development. 

IECJG First hour of maximum oviposition after 

sunset. 

T1 Minimum temperature (°F) for corn and 

sorghum development. 

T2 Maximum temperature (°F) for corn and 

sorghum development. 

TEM Effective hourly degree-days for Heliotkis. 

TEMOVR Number of degrees hourly temperature 
(*F) is over 02. 

THP Actual or expected daily high temperature 

(° F). 

TLP Actual or expected daily low temperature 

TF). 

TM Actual or calculated hourly temperature 

T F). 


TRKGRM 

FERATE Fertility rate of Trichoyravima . 

IAGE Age of Tricho gramma adults in days. 

It Mean development time for Trichoyramma 

in days. 

TS Keeps track of number of parasitized 

He Hot his eggs, 

X Mortality factor. 

INPUT REQUIREMENTS 

Device 05 or 
SYSIN: KODE 
FORMAT (2014) 


KODE Col. No. Description 

1 4 0 Leaves population input the 

same as in previous run. 
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1 Reads in moths (total num- 
ber) . 

2 Beads in eggs (1-clay-old). 

3 Beads in moths (male and fe- 
male) . 

4 Beads in eggs (all ages). 

6 Migration used as input. 

6 Reads in total eggs with miss- 
ing data. 

7 Beads i .11 0- to 1- day-old eggs 
with missing data. 

2 8 0 Temperature stays the same 

as in previous run; genera- 
tion length stays same. 

1 Reads in hourly temperature. 

2 Temperature stays the same 
as in previous run; genera- 
tion length can be changed. 

3 Reads in daily maximum and 
minimum temperature. 

Note: Always read in temperature for at least 1 day 
before and 30 days after time of simulation of 
crop or insect. 

3 12 0 No moonlight effect on insect. 

10-12 xxx Day number of first full 

moon in January (add 1 for 
leap year) . 


36 0 No adjustment in trap catch. 

1 Adjusts trap catch to moths 
per acre. 

2 Uses same adjustment as pre- 
vious run. 

38-40 xxx Day number that model 
should run to. 
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44 1 Corn. 

2 Sorghum. 

3 Cotton. 


48 0 Leaves predator input the 

same as in previous run. 

1 Reads in maximum predation. 

2 Reads in number of predators 
per day. 

3 Reads in precentage of preda- 
tion on each day. 

4 Reads in number of predators 
with incomplete data. 

52 0 Leaves natural mortalities 

the same as in previous run. 

1 Reads in natural mortalities, 
egg fertility, and the maxi- 
mum number of eggs per 

moth per day. 


4 

16 

0 

No mortalities from insecti- 

14 

56 

0 

Leaves crop phenology the 




cides. 




same as in previous run. 



1 

Reads in insecticide mortality. 



1 

Reads in crop phenology. 



2 

Leaves input the same as in 
previous run. 

16 

60 

0 

Leaves migration the same as 
in previous run. 

5 

20 

0 

Leaves input the same as in 
previous run (0 to 1st run). 



1 

Calculates migration. 



1 

Reads in maximum para- 

16 

62 

0 

Log plot of data. 




sitism. 



1 

Linear plot of data. 



2 

Reads in number of parasites 
on each day. 



2 

No plot, 



3 

Reads in percentage of para- 
sitism for each day. 


63-64 

XX 

Number of plots wanted. 



4 

Roads in number of Tricho - 
U rmmna and generates fu- 
ture populations. 

17 

68 

0 

1 

No cannibalism, 

Calculates cannibalism of lar- 








vae. 

0 

24 

0 

Migration not calculated. 







1 

Migration away from crop. 

18 

72 

1 

This run is for H . zea and H. 



2 

Migration into crop as input. 



2 

viresccns , 

This run is for H , virescena 

7 

26-28 


Not used. 



3 

only. 

This run is for H. zea only. 

8 

32 

0 

No smoothing of data. 

19 

76 

0 

No cloud cover input, clear 



1 

Smooths trap catch by a dou- 




sky assumed. 




ble 3-day moving average. 



1 

Input cloud cover and use. 
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Device 05 or 

SYSIN: Mortalities, egg fertility, and egg 
laying 

KODE { 13) =1 or more 
FORMAT (2 (6F4.3, F4.0) ) 


Column No. 

Description 

1-4 

H.zea: Natural mortality of eggs ( f /c-r 
100). 

5-8 

H.zea: Natural mortality of 1- to 3- 
instar larvae. 

9-12 

H.zea: Natural mortality of 4- and 5- 
instar larvae. 

13-16 

H.zea: Natural mortality of pupae. 

17-20 

H.zea: Natural mortality of adults. 

21-24 

H.zea: Egg fertility. 

25-28 

H.zea: Maximum number of eggs per 
female per day. 

29-32 

H. virescens: Natural mortality of eggs. 

33-36 

H. virescens: Natural mortality of 1- to 
3-instar larvae. 

37-40 

H. virescens : Natural mortality of 4- 
and 5-instar larvae. 

41-44 

H. virescens: Natural mortality of 
pupae. 

45-48 

H. virescens: Natural mortality of 
adults. 

49-52 

H. virescens : Egg fertility. 

53-56 

H. virescens: Maximum number of eggs 
per female per day. 


Device 05 or 

SYSIN: Days of population data on device 11 
KODE (1) -lor more 
FORMAT (214) 

Column No. Description 

1-4 First day of initial population data to be 
read in, 

5-8 Last day of initial population data to be 
read in. 

Device 05 or 

SYSIN: Days of temperature data on device 15 
and device 16 
KODE: (2) =1 or 3 
FORMAT (214) 

Column No. Description 

1-4 First day of temperature data to be read 
in and degree-days calculated. 

5-8 Last day of temperature data to be read 
in and degree-days calculated. 

Note: Reading switches automatically to device 16 
when end-of-file encountered on device 15. Suf- 
ficient days of input should be allowed for crop 
growth and multiple-run situations, 


Device 05 or 

S Y SIN : Crop phenology 
KODE (14) =1 

FORMAT (416, F6.0, 16, 5F3.2) 


Column No. 

Description 

1-6 

Day number of planting date of crop, IP. 

7-12 

Day number of emergence date of crop, 
IE. 

13-18 

Day number when first fruit appeared, 
IF. 

19-24 

Day number when crop became mature, 
IM. 

25-30 

Number of plants per acre. 

31-36 

Day number when crop becomes attrac- 
tive to moths, IA. 

37-39 

Percentage of the area in corn. 

40-42 

Percentage of the area in sorghum, 

43-45 

Percentage of the area in cotton. 

46-48 

Percentage of the area in pasture, 

49-52 

Percentage of the area in other crops, 

Device 05 

or 

SYSIN: 

Maximum parasitism 

KODE (5) =1 or 4 

FORMAT 

(2F6.3) 

Column No. 

Description 

1-6 

Maximum parasitism of eggs (%-^100), 
Set to zero if KODE (5) —4. 

7-12 

Maximum parasitism of 1- to 3-instar 
larvae (%-i-100). 

Device 05 

or 

SYSIN: 

Maximum predation 

KODE (12) — 1 

FORMAT 

(2F6.3) 


Column No. Description 

1-6 Maximum predation of eggs (%-rlOO). 
7-12 Maximum predation of 1- to 8 -instar 
larvao (%-rl00). 

Device 05 or 
SYSIN: Plot control 
KODE (16) less than 200 
FORMAT (3612) , 1 card required for each plot 
wanted (10 curves per plot maxi- 
mum) 

[Put a 1 in column indicated for desired plot] 


Column number 


Scale 1 

Information plotted 

H. 

zea 

h\ 

virescens 

Both 

spp. 

1 

Total eggs 

2 

26 

50 

1 

Small larvae (1- to 

3 -instar) 

4 

28 

62 

1 

Large larvae (4- to 
5-instar) 

6 

30 

64 
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[Put a 1 in column indicated for desired plot] 


Scale 1 

Information plotted 

PL 

zea 

Column number 

IP. 

viresce ns 

Both 

spp. 

1 

Pupae 

8 

32 

56 

2 

Preovipositing 
adults 

10 

34 

58 

2 

Ovipositing adults . . 

12 

36 

60 

2 

Input data (as 

adults) 

14 

38 

62 

2 

Migrating adults 

16 

40 

64 

1 

Parasitized eggs .... 

18 

42 

66 

2 

Predicted ovipositing 
adults equivalents - 

20 

44 

68 

1 

Other variable 
scale (1) 

22 

46 

70 

2 

Other variable 
scale (2) 

24 

48 

72 


] 1 based on maximum eggs; 2 based on maximum 


preovipositing adults. 


13-18 PI. zea : Number of females. 

19-24 H. viresccns : Number of males. 

25-30 PP. virescens: Number of females. 


B. KODE (1) =2: Egg input (0- to 1-day-old) 
FORMAT (9x, 13, 2FG.O) , 1 card for each day 


Column No. 
10-12 
13-18 
19-24 


Description 
Day number. 

//. zca : Number of eggs. 

IP. virescens : Number of eggs. 


C. KODE < 1 ) — 3 : Moth input (initial population = 
2 times number of males, sex 
ratio calculated ) 

FORMAT (10, 4F6.1), card for each day 


Same as A above 


*D. KODE (1)— 4: Egg input (all ages) 
FORMAT (9x, 13, 2FG.0) 


Device 05 or 
SYSIN: KRUN 
(always required) 

FORMAT (13) 

KRUN designates the status of additional runs to be 
made. Values greater than aero indicate that additional 
runs are to be made, and a value of zero (blank) indi- 
cates that the present input is the final run to be made. 

If the migration outputs from several runs are to be 
averaged for input to a second crop, special sequential 
values should be coded as follows: for corn, 11, 12, 
13, ... 19; for sorghum, 21, 22, 23, , , . 20; for cotton, 
31, 32, 33, . . . 39. 

KRUN indicates the field number associated with the 
data following the KRUN card. The value of KRUN 
for the first data set entered is assigned automatically 
by the program depending on the crop type selected. 
For a series of three corn fields and two sorghum fields 
followed by a cotton field the input sequence would be 
as follows: corn data 1 with KRUN = 12, corn data 2 
with KRUN = 13, corn data 3 with KRUN — 21, sor- 
ghum data 1 with KRUN — 22, sorghum data 2 with 
KRUN — 31, cotton data with KRUN = 0. The average 
migration from the three corn fields would then be 
added to the average migration from the sorghum and 
used as input to the cotton with the proper selection of 
constants on the KODE cards, If independent multiple 
runs are to be made, values of KRUN from 1 to 9 should 
be used. No averaging of output will take place, 

Device 11 : Moth or egg input 

A. KODE (1)— 1: Moth input (initial population = 
males + females) 

FORMAT (IG, 4F6.1), 1 card for each day 

Column No. Description 

1-G Day number. 

7-12 PI, zea : Number of males. 


Same as B above 


E. KODE (1)=5: Migration of H . zea moths from 
corn or sorghum or both; moth 
input of H . virescens if KODE 
(IS) — 1 or 2 

FORMAT (IG, 12x, 2FG.1) 

Same as A above; enter zero for H . zea. 


F. KODE (1) — C: Egg input (all ages) with miss- 
ing data 

FORMAT (18x, 2FG.0, 13), 1 card for each day 

that data are avail- 
able; computer will 
calculate missing 
days. Cards must be 
in order by date. 


Column No. 
19-24 
25-30 
37-39 


Description 

IP, zea: Number of eggs. 

II. vh'escens: Number of eggs. 
Day number, 


G. KODE (1) — 7: Egg input (0- to 1-day-old) with 
missing data 


Same as F above 

Device 12: Insecticide mortalities 
FORMAT (16, 14F4.3), 1 card for each day 

there is mortality 


Column No. 
1-G 
7-10 
11-14 
15-18 
19-22 
23-26 
27-30 


Description 
Day number. 

Il.zea: Parasites (%-rlOO). 

H.zea: Predators. 

PI. zea: 1- to 3-instar larvae. 
PI. zea : 4- and 5-instar larvae. 

PI. zca: Pupae. 

PI. zca: Adults. 
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31-34 

H, zea : Eggs, 


35-38 

H. virescens: 

Parasites. 

39-42 

H. virescens : 

Predators. 

43-46 

PL vires cam : 

D to 3-instar larvae. 

47-50 

H. virescens: 

4- and 5-instar larvae. 

51-54 

PL virescens: 

Pupae, 

55-58 

PL. virescens : 

Adults. 

59-62 

PL virescens : 

Eggs. 

Device 13: 

Parasites 



A. KODE (5) =2 
FORMAT (16, 3F6.0) 

Column No. Description 

1-C Day number. 

7~12 Number of Tricho gramma released. 
13-18 Number of natural egg parasites, 
19-24 Number of natural larval parasites. 

B. KODE (5) =3 
FORMAT (IG, 2F6.3) 

Column No. Description 

1-G Day number. 

7-12 Parasitism of eggs ( r /r -r 100) . 

13-18 Parasitism of larvae ( r /c -s-100). 

C. KODE (5) —4 
FORMAT <!G t FG.O) 

Column No, Description 

1-6 Day number. 

7—12 Number of Tricfiogramma. 


Device 14: Predators 


C. KODE (12) =4 

FORMAT (6X, GF5.0, 13) (Special format for 
handling 1974 Pecos, Tex., data.) 

Column No. Description 

7-11 Number of spiders. 

12-16 Number of adult lady beetles. 

17-21 Number of other Hemiptera. 

22-26 Number of Chrysopa larvae. 

2 7-3 1 N u m ber of 67 * rysopa adults. 

32-36 Number of Coleoptera. 

37-30 Day number, 

Devices 15 and 16: Air temperature 

A. KODE (2) = 1: Read from device 15 

FORMAT (2X, 13, 3X, 24F3.0) : Measured hour- 
ly temperature, 1 card per day 

Column No. Description 

3-5 Day number. 

9-11 Average temperature from 0000 to 0100 

hours, 

12-14 Average temperature from 0100 to 0200 

hours. 


78-80 Average temperature from 2400 to 0100 
hours. 

B. KODE (2) =3 or end of file on device 15: Read 
from device 16. Sufficient data 
corresponding to final day speci- 
fied for temperature input on 
SYS IN must be provided, 
FORMAT (13, 2F3.0) : Measured high and low 
for each day, 


A. KODE (1:2) =2 

FORMAT (9X, 13, GF6.0, 6X, 2F6.0, 6X, F6.0) 
(Special format for handling 1973-74 
Frio County, Tex., data.) 


Column No. Description 

1-3 Day number, 

4-6 Measured high temperature for day. 

7-9 Measured low temperature for day. 


Column No. 
10-12 
13-18 
19-24 
25-30 
31-36 
37-42 
43-48 
55-GO 
67-72 
73-78 


Description 
Day number. 

Number of Scyrmtus. 

Number of other Coccinellidae. 
Number of other Coleoptera. 
Number of Onus. 

Number of Geocoris. 

Number of other Hemiptera. 
Number of Chrysopa im matures. 
Number of spiders. 

Number of released Chrysopa, 


B. KODE (12) =3 

FORMAT (16, 2F6.2) 


Column No. 
1-6 
7-12 
13-18 


Description 
Day number. 

Predation of eggs ( -I-IOO) , 
Predation of larvae ( -^100) . 


Device 17: Cloud cover 
KODE (19)— 1: First and last day read in 
correspond to temperature 
data. On end-file condition 
program uses last data read 
in for remainder of days of 
input requested. Final card 
should show desired cloud 
conditions accordingly, 
FORMAT (2X, 13, 3X, 24F3.0) 

Column No. Description 

3-5 Day number. 

11 Cloud cover from 0000 to 0100 hours. 

12-14 Cloud cover from 0100 to 0200 hours. 


78-80 Cloud cover from 2400 to 0100 hours. 
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Device 18: Extra data to be added to print 
plot for comparison with model 
output. 

FORMAT (6X, 16, F6.0) 


Column No. 
7-12 
13-18 


Description 
Day number. 

Data value. 


- LCAC f'GCtILE 


* ** * Ji \ ******** * 

* gc sir p * 

►IlXFCUT \cu CATA * 

* SL^UIKEURVI 5 * 

* * * * * * **.**•*«« 


CCN inrti HPLT • 


♦ IMtlAL ** 

♦ PtlPLLATICN * 

♦ ISPL1 * 

***** 4*414******* 


12 CO 2 4 1 * 


CITE 

4*4 1*4 L * .* 


C4ILY 
flA si TI !« 
IAPL1 


4 * 444 ** 444 * 4*444 


ei • . 

.« *. 

*rnem-o 


C3 «, 

.+ *, 

4 4 

XCDEt'*»*l 

**» ,4* 

* NC 


...... 

• *« 

Ktoe<5i=o 


» *.*C2* 44 (4 (4 4* 4 03 *, 

CEVJCO-14 COAFM .* ♦. 

« NC .* « 

rmv pfepaticmk. xcoE<izt*o 

IMPLT * 4. .* 


TFPPECATLRE • 
IA PL T 4 


.« Gl *•*. 
.^ENn-nF-FTtF*. 
«. CN « 

'".OEVICf-lS.* 


YES .» 

1. ...... ♦*. KCDt- <21-1 


6 3 


♦OEVICE-17 COAA-i* 

• noinn 

* cnvEF 


noinn ciclc *x.. *,* xcnMniM 

'* ni t:P tNPlT 4 


***** J2**4 ******* 

♦ DEVKE-te C04fA* 

i”'1hiWo ;5r i x 


KODEm.l 

I, ,4* 

» 

...x. 


H3 


EXTRA 4. 

Pinr data ** 
ftFQUFSTEC*-* 


*4 4*K3 444 *44* ** 

* « 

* E NO CF t AT A * 

* * 
4 * ***** •* ««»■** 
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worniva - load module 

INPUT REQUIREMENTS 


SI'S IN DEVICE- U DEVICE-12 

• 4 4 4Al •" •*«**» wMAItUMM** ""AS"******* 

*f? AS ] C CAT* ANH * * INITIAL * * DAILY * 

♦ CCNTFCL !U=UT * * MQPULAT | ON • * INSECTJCIOL * 

» • * t NPU7 * * * 


■PI 


INPUT KOBE 


• INPUT PALES t • YES « 

FEMALES 

* » 


X 

* t» 

04 »* 

* *4 

«' 

NODE U 1 = 1 . * 

a * 

* • • * 

♦ a l * 

* NU 


X 

* INPUT DAILY * 
MONT AL I TIES 


J< 


■ | « t « t " Mt*t 

• vcrtauties, * 

EGG F« M T ] LI Tr i X 
* EGG LAY # 

V 


• C2 

♦ ’ 

MODE! 1 3 > = 0 


x! 

X 


It" "«>*•»* 

► DAYS OF * 
POPULATION DATA Xmim 
*cn n=vrc=--ii * 


Nfl 


02* * »♦ 

*» *t 

KCDF.I U-0 !* 




"*.**»£ I ItlMMtlM 
CAYS JF T E YD 
* PATACN * 
DEVICE- 15 ANT> X. 
* CfVJCC-16 * 

0 



NO • ♦ «. 

«♦». KODC(a)-c ,* 
*• or a ► * 

*• * * 


It Ml llllMM l« MMI IMllXi 

X 

F2 " 

i» «• 

*C«CP PHEKCLCGV * ND ► * 

CATA X KDDE(lAl-0 

Y • *. I 

* »l |t 

tilt MMMiil *• i * 

• * YES 


i** •••Cl ** 

»«AA PARASITISM * YES . 

FCfl cGGS AND X* •*«■•»■ *» 

« small larva * 

r 



K0t>E<5 I = I 

*• OR * 



;• 



J 




a *, 


***« .*C3* ******* 


C4 *. 




, t t. 


* 

* 

YES ,* *. 


INPUT EGGS 

X. ■ . 



* 

9 

*. OT 4 . * 


* 


*. .* 




*« . * 


f 


* NO 

X 

• 


■ 

4 

* * 



X 




* 



f 14 

DA *4 




1 * *» 


• INPUT MALLS, 

* 

YES .* *. 


FEMALES, calc 


""*4 K.0t)EU1=3 t* 


* SP ic RAT rn 

9 

*• . • 


-* 


* . 4* 




*• . * 


• 

1 


* NO 

X 




* p 



X 




a *. 


t*t«,E]*>tt*M* 

* 9 

E4 *. 


• OBTAINS 2a. A 

* 

* * *• 


* input rnew 

• 

res a* *. 


* MIGRATION OF 

* x* * • 



• PRE VIOLS RUN 

• 

*• * * 


« 

» 

* • • ♦ 



• • 

*• # * 


* 

9 

1 


* NU 

• 

9 


9 

X 


■9 



*»* 



♦INPUT YtPESCENS* 



ADULTS Ip 




* NF.FOtn 

* 

9 


♦ 


9 

9 

9 

X 







X 


MlM*G3Ut**«m,) 

G*' 






•EGG INPUT VITH 

9 

YES i* *< 


MISSING DATA 

Km 





*. OR 7 a * 


Y 


** 




*. t* 


* 


* 

!. 





m2 »» 

• MAX FPICATION * YES • ♦ *. *** 

CF ECCS AND x I, KCOE, (12 1 - I «« 

♦ SMALL L Hv* * *, ,a 


«* Mh3*»« *4*4* A 
* END OF CAT A ON « 
* OEVrCE-Jl * 


* AO 


X 

♦♦♦♦C5**«*»**** 
*ENO O' 5 DATA ON ♦ 
* DEVJCE-12 * 

< * 
«**• *♦**"""• 


01. FORMAT (2014) 

CU FORMAT 
IZI 6F4«3»F4»0J J 

Dl» FORMA T ( 2 [ A ) 


• * « 

J2 », 

I* •» 

* PiLCT CCMFQL * YES «t a* 

C A K C ! S J - C N 5 PER KQ0E<1«] .« 

* FLCT • *..LT. 200 .* 

0 *t ,i 

MtillHItlM *# ,t 

• * NO 



X 


• l"**K|«" a* .4. ♦*• 

* * 

K C L N a i 

' * * 

Jf 

MMHIiimi 


•***K£*** * **"• 

*ELO OF INPUT ON* 
;* sys IN - * 
* * 
* it** 


El« FORMAT (214) 

Ft* FORMAT 

(41 6.F6.0 ♦ IG.5F3. 2) 

Git FORMAT! 2F6* 31 
Hi* FORMAT! 2F6.3) 

Jl. FORMAT! 3612) 

K|» FORMAT! 13) 

S3* FORMAT! 16, 4F6. 1) 

C3. FDRMAT(9X« t3t 2F6*0) 
03* FORMAT ( I 6 • 4F6* 1 ) 


0S. FORMAT! !&• 14F 4. 3) 
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THZV2 - LC AH MODULE 
INPUT REQUIREMENTS 


* GAILY * 

* PARASITISM * 

* ♦ 
+ *+«*«* ******** 


* ft + ft * * g ^ ft * ** 4 * * * ft* * 

* INPUT NUMBER OF* 
EGG AND LARVAL X, 
* PAR AS 1165 * 

ft 

************* 


82 ♦ . 

*♦ •. 

KCDE C 51 »2 


DEVICE-15 

1 * 4 * Ml ******** 

* HOURLY * 

* TEMPERATURE * 

* * 
<» 44 t*M**«**«* 

* 

X 



* IKPUT HOUR IV * 
TEMPERATURE 

* 9 

A 


****44 ********* 

* HOURLY CLOUD * 

* COVER * 


INPU T HOURLY 
CLCUO COVER 


* 9*****499949 


44*999C\44**4*9**99 

* INPUT PERCENT * 
EGG AND LARVA X. 
* PARASITISM * 
ft 

99*94944 4*9*4 


9 *9*4 9$ 14449444*49* 

♦INPUT NUMCER OF* 

TR I CHDQR AMHA X. 

* * 

, ft 

944444**4**49 

........ ,x* 


C2 ♦. 

.* ». 

KCCE ( 5 1 =3 

’ft. .** 

* NC 


YES .♦ *• 

KCOE I 5) *4 t .4 

"*♦ * 4 * 


**■4 *C3 * **4499** 

~'" c oV 


* e ~ c oie.Eiif 5 ch • 


*EN0 OF DATA ON * 

* DEV ice - 17 * 

* * 
»**« *********** 


*«**£ 1 *»****+* * 

* ENO OF DATA ON * 

* CEV1CE-13 * 

* * 
*************** 


oevicc-u 

•* 4 *e: 

: # MK 


3 *MU***t 

Y FAX-KIK + 
PE RATI RE * 


•EXTRA DATA FOR =* 
* PLOTTtNG * 


f!EV ICE- 1 A 

♦♦♦♦F l***ft*ft>**» 

*C A I L Y PREDATION* 
* * 
*************** 


i*4* 4 4fj4**44**4‘4*4 


IKPUT MAX-KIN ♦ 
TEMPERATURE 


4 4 4 ***** * 4 49* 


944*4499***** 


* 4*99 4 Q 1*4 *4 4*44 4 44 


EDATORS BY X. 
* CLASS ♦ 
ft 

444*444*4*494 


X 

G2* ** *, 


4. .9 
* NC 


• •••G3***Mft 

•END CF CATA 
* DEVICE-16 


44 4 4 4* 94 *4 * *4 ft* 


***«C 4 ** ******* 
•ENO OF CATA ON * 
* OEV ICE-18 * 

**944***9*49949 


• •ft*** HIM* ******** 

* INPUT PERCENT * 
PRECATION X* 
* * 
ft 

*444449944499 


ft* ****jl*ft* ******** 


44*4444944444 


YES *♦ *, 

KCOE I I 2 I * 2 .ft 
«* 


ft. .♦ 

* NC 


J2 *, 

. * *. 

YES .ft *« 

*,.,*. KCCEM2I-4 ♦* 


Bl. 

FORMAT! 16, 

1 3F6 

.01 

Cl, 

FORMAT 1 16 1 

f 2F6 

.31 

01, 

FORMAT C 1 6 i 

i F6 * 

01 


Gl, FOR FflTO CC* DATA 
FORMA T|9Xi I3*6F6.0i6XrZ 
F6«0r6X|Fb»0| 

HI, FORMAT C 16, 2F6 . 2 ] 


****Kl********* 

* END OF CATA ON * 

* 0 EV I CE -14 * 

* * 
*************** 


I 6 X, 6 F 5.0 , I 3 ) 

FofiNATIZX, T3t 3X. 2 AF 3 .0) 
F3, FCRKATtI3.2F3.0l 

FORMAT (2Ki 1 3*3X* 2A F3 ,0> 
FA. FORK AT (6Xit6 T F6»0l 
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PROGRAM LISTING 


1 

2 

3 

4 

5 

6 
7 
3 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
2 3 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 
45 

47 

48 

49 

50 

51 
£2 

53 

54 

65 
56 
K 7 

55 

59 

60 
61 
62 

63 

64 
£5 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 


MODEL OF THE HELIOTHIS ZE A AND VIRESCENS BY A . W . HART STACK 
AGRICULTURAL ENGINEERING BLDG. USDA ARS COLLEGE STATION* 


D I ME NS 1 0 N 
DIMENSION 

dimension 

DI MENSION 
D I MENS ION 
DI MENSI ON 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 


1 

C 

9 

C 

19 

C 

29 

C 

39 

C 

C 

c 


MQTHZV 2 
ROOM 231 
TEXAS 

REAL NITMIN 

REAL SLPOP (370)* LLPOP(370)t PPaP(370), PQAP0P(370) 

COMMON KODE (2 0) .FACTO , SDEV( 370 . 4 ) . NI T M I N ( 3 7 0 ) . ED < 37 0 ), EDCORN ( 3 00 ) , 
ASQUARE! 3 70) * BOLLS! 370 ) • = ARS ( 37 0),HE4DS( 37 0) .S<4) , POP <37 0 ) *E GLY TM < 3 
C70 ) . I FIRST , I LASTC . IH, IV . A DMO ON ( 37 0),SAREA(370). TCLOUDC 370 ) , I LA ST . 

D KR UN , SE X RTO( 3 70,2) , EGPEMT ( 2 ) ,PCTCRP<5) , I GEN < 37 0 ,7 , 2 ) , TEMAV G ( 3 7 0 ) , 

E EQELAD ( 3 70 ♦ 2 ) , A DT SMG ( 370,2, 2 i, TOTAL! 370,9,2) ,PO$( 370 ,7 , 2) , PAR SI T 
8(370.2), PREDOR ( 37 0 ,2 ) .PRMGRT (3 70 ,2 ) , I A, I M,L, I , DE ATHY * SURV , EGL YPR{ 

G 37 0, 2) .PLANTS .KCROP, I S , S I TE S ( 3 69 > , P , DF , EGLYC R ( 3 70 ,2 ) *CANL AR , J J , KK , 
HNNtXPOP, N,PP * IN,K,SRV IVE, X, IPX, IGKODE! 12, 3 ) , OAT X( 370 ) ,TP,DEATH< 5, 2 
I ) , XBAR , X SD * D I STS, I FI STT, I LASTT , I E » PR6L A R , PR BNE W 
DIMENSION 5S< 2) 

ZZE RO (2590 ) , ZONE! 1 628 0 ) ♦ ZTWO ( 74 0 ) , ZT HREE ( 1 P54 ) 

EGGFRT! 2 ) , E GP( 2 ) , T OT MG S ( 37 0 , 2 ) .CLOUD (24) 

FM( 30 ) 

PDE RTE( 7) fPOLRTE! 7) ,PRED(7) 

TME A ( 80 0 ) 

DE A UL T (2 5 ) 

OEGDA Y ( 6,2 ) 

ADT S ( 370 ,2 ) 

PC T MG ( 2 ) 

SITEOC! 370) 

EQUI VALENCE ( 2 ZERO (1 ) , $ L POP ( 1 ) ) , < ZZERO ( 37 1 ) , LLPOP ( 1 ) ) • 

D ( ZZERQ{ 741 ) , PPQP( 1 ) ) , ( ZZEROd 11 1 ) ,P0AP0P<1) ) , 

E ( ZZEROU4B 1 ) ,S ITEOCM ) ) , ( ZZ ER 0 ( 1 8 5 1 ) , TO TMGS ( 1*1)) 

EQUIVALENCE ( ZONE ( 1 ). EQE LAD ( 1 . 1 ) > 

EQUIVALENCE ( ZTWO ( 1) . ADMOON ( II ) 

EQUIVALENCE ( ZTHREEU ) , S QUA R E ( I ) ) 

DO 1 1=1,370 
TCLGUD ( I ) = 0 .0 

INITIALIZE TO ZERO VARIABLES NOT IN COMMON. 

DO 9 1=1 ,2590 
ZZERO(I) =0.0 

INITIALIZE TO ZERO VARIABLES EQELAD THRU PRMGRT IN COMMON. 

DO 19 1=1,16280 
ZONEU ) = C.O 

INITIALIZE TO ONE VARIABLES AD MOON THRU SAREA IN COMMON. 

DO 29 1=1 ,740 

ZTWO! I ) =1.0 

INITIALIZE TO ZERO VARIABLES SQUARE TH3 U POP IN COMMON. 

DO 39 1=1 ,1854 
ZTHREE! I )= 0 • 0 

KRUN-0 MEANS THE INITIAL SIMULATION. 

ADDITIONAL SIMULATION. IFIRST-DAY NO, 

-DAY NO OF FINAL DAY OF SIMULATION. 

IOX=0 
KR UN = 0 

KODE-i-20 OPTIONS. 

34 PE AD (5,500) (KCDE{ I ) .1=1 ,20) 

500 FG RM AT (2014) 

ILAST=K0DE(10 ) 

IH-CODE for HELICTH1S ZE A- IV CCDE FOR H ELIOTHI5 V I RES CENS . 

I H — 1 
I V=2 

IF(KODE< 18). EC. 2) I H = 2 
IF(KODE< 18) .EG. 3) IV=1 
IF(KRUN.EG.O) GO TC 7 
IF ( KODE( 13 ) .EQ. 0) GO TO 37 

EGGS.SMALL LARVAE.LARGE L AR V AE . P UP A E. 
TrEAD<! T* ,? F , E t GS> 1 WPCMT-EGGS PER FEMALE MOTH. * 

506 FOR^ATt 2(lF4. ( J * ’ 1 = 1 ’ 5 ‘ ‘ FGGFRT( J EGPEMT! J > . J =1 . 2) 

DO 38 I = I H , IV 
EGP( I) =EGPEMT (I ) 

EGPEMT ( I )=EGP EM T( I )/2 • 

IF (KODt ( 1 ) .EQ .0 ) GO TO 60 
EGG IMP- SUBRQUTI NE FOR INPUTING INITIAL 
VARIOUS POPULATION STAGES. ILASTC-LAST 
I .^-, F i£I£“ L !. GHT TR * P CONVERSION FACTOR. 

INTO CROP. SEXRTO-SEX RATIO MALES-FE MALES . 

OVERALL AREA IN EACH CROP S I MJL AT ED . 

RE AO ( 5*5 CO ) 1FIRST.IL ASTC 
CALL EGG INP 

IF(KRUN.EQ.O) GC TC 62 
KF=IFIRST-1 
00 61 I=KF,1LAST 
SLP0P( I ) =0. 

LLPOP ( t )=0 . 

PPOP ( I ) *0 . 


KR US ~ 1 
F IRST 


MEANS GREATER THAN 1- 
DAY OF SIMULATION. I LAST 


38 

37 


POPULATIONS. TOTAL-ARRAY OF 
DAY OF INITIAL INPUT DATA READ 
ADTSMG-ADULTS MIGRATING 
PC TCRP— PERCENTAGE OF 


60 
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8 3 

84 

85 

86 

87 

88 
8 ? 
SO 

91 

92 
S3 

94 

95 

96 

97 
95 
99 

l 00 
1 01 
1 02 
l 03 
1 04 
1 05 
1 06 
1 07 
1 C 8 
1 09 
1 10 
i n 
112 
113 
1 14 
1 15 
1 16 
117 
1 IS 
1 19 
1 20 
IS 1 
122 
1 23 

124 

125 
1 26 
1 27 
1 28 
1 20 
1 30 
121 
1 32 
l 33 
1 34 
1 35 

126 
137 
1 30 
1 39 
140 
1 41 
142 
1 43 
144 
1 45 

146 

147 
1 48 
1 49 
1 50 

151 

152 
1 53 
l E4 
l 55 
1 56 
1 57 
1 58 
1 59 
1 60 
1 61 
l 62 
1 63 
1 64 
165 
1 66 
l 67 


PQAP OP < I )t= 0. 

C AT X (I ) = 0 • 

DO 61 K = IH ,1 V 
TQTMGS (I,K )-0 .0 
E QEL AD < I ,K ) -O .0 
TOTAL (I , S , K> =0 . 0 
TOT AL ( I . 9* K)asO . 0 
DO 61 J=I , 6 

61 TOTAL ! I . J ,K ) =0 .0 

62 CONTINUE 

IF! KODE! 2 ) *EO • 0 • QR« KODE! 2) • EQ.2 ) GO TO 2 
C READS IN TEMPERATURE 

C TEMP-SUBROUTINE USED TO READ IN TEMP AMD OTHER VARIABLES* 

C ED-ACCUM ALATED DEGREE DAYS USED FOR INSECT + COTTON- T E MAV G~ AVE RA GE 
C TEMP PER DAY * EO CO RN- A CCU MA L AT ED DEGREE 0 AYS FOR CORN* 

C EGLYTM-TEMP DUPING THE PERIOD WHEN THE MOTH LAYS EGGS. 

C NITMI N-AV ERAGE TEMP DURING THE NIGHT- TCL D UD- A YEP AGE CLCUD C OVF R 
C DURING THE PERIOD WHEN THE MOTH IS LAYING EGGS. 

READ ( 5,500) IF I S T T , I L A S T T 
CALL TEMP 

C GEN L I N- S UR R OUT I NE USED FOR CALCULATING DEVELOPING TIMES FOR ALL 
C STAGES OF INSECT. IGEN- ARRAY OF DEVELOPMENT TIMES FOR VARIOUS STAGES 
C OF INSECT. IGENZ- ARRAY OF DEVELOPMENT TIMES *0P HE L I Q T H I S ZEA. 

C IGENV-AftRAY FOR DEVELOPMENT TIMES FOR HELIGTHIS VIR. 

C SOEV-STANOARD DEVIATION OF DEVELOPMENT TIMES VARIOUS STAGES. 

2 CALL GEN LN 

C CROP-SUBROUTINE FOR CALCULATING CROP PHENOLOGY. £ GL YPR- EGG LA Y 
C PROBABILITY. 

C PL4 NT $“ NO . OF PLANTS PER ACRE. 

C SQUARL-NO. CF SQUARES PER ACRE. 

C BOLLS- NO OF 30LLS PER ACTE. EARS- NO OF EARS PER ACRE. 

C HEADS- NO OF SORGHLM HEADS PER ACRE. KCROP- KODE FOP CROP. 

C I A— CODE FOR 1ST ATTRACTNESS OF CROP. I M - DAY WHEN CROP PEACHES 
C MATURITY. 

C SAREA-SEARCK AREA FOR PARASITES AND PREDATORS. IS-DAY NO WHEN CROP 
C IS MAX ATTRACTANCE. SITES-NO. OF SITES FDR EGG LAY ING. 

C DF-DI SPERSI ON FOR NEGATIVE BINOMIAL. EG L Y C R- AT T RA CT ANCE OF CROP FOR 
C EGG LAY. 

CALL CROP 

IF ( KODE( 5 ) .EO .0 ) GO TO 33 

C PAR I M P-SUBR GU TINE FOR CALCULATING PAR AS I T I Z AT I ON . 

C PARS I T-PtR CEN T P AR A S I T I Z A T I 0 N CALCULATED. 

32 CALL PAR I N P 

33 IF( KODL{ 12) .E Q. 0) GO TO 35 

C PRO IMP-SUBROUT INE FOR PREDATION INPUT. PREDOR- PERCENT PREDATION 
C CALCULATED. 

CALL PRD INP 

35 CONTINUE 

5 IF< KODE. ( 3) .EQ . 0* GC TO 3 

C MODN-SUBROUTI NE FOR CALCULATING REDUCTION OF EGG LAY DUE TO 
C MOONLIGHT. 

C ADMOON-PERCENT REDUCTION OF EGG LAY OF MOONLIGHT CALCULATED. 

CALL MOON 

3 CONTINUE 

C SPRAY-SUBROUTINE FOR INPUTING INSECTICIDE MORTALITY* 

C P OS-PERCENT MORTALITY DUE TO INSECTICIDE. 

CALL SPRAY 

63 IF {KOOE< 9) .EQ .0 ) GO TO 36 

C F LD POP- SUOR OUT INE FOR CONVERTING LIGHT TRAP CATCH TO MOTHS D ER ACRE. 

C FACTO- CAL CUL A TED FACTOR FOR CONVERTING LIGHT TRAP CATCH. 

CALL FLDPOP 
GU TO 64 

C FACTD-EQUALS 1 IF INPUT IS NOT LIGHT TRAP CATCHES. 

36 FACTD^l . 0 

64 CONTINUE 

C NEXT SERIES OF CAROS PRINT OUT READ IN AND CALCULATED INPUT D4TA- 
C FOR USE IN SIMULATICN. 

IF(KODE(18).EQ.l)WRITE(6, 70 8) 

IF { KODE< 19 ) *£ C.2) WRIT E<6 ,710 ) 

IF ( KODE < 18) .EQ. 3) WRITE! 6,707) 

707 FORMAT ( * INPUT FOR THE MODEL HAS BEEN HELI3THIS ZEA ONLY* ) 

700 FORMAT!* I NP UT FOR THE MODEL HAS BEEN HELIOTHIS ZEA AND VIRESCEN* 
X I 

710 FORMAT!* INPUT FOR THE MODEL HAS BEEN HELIOTHIS VIRESCEN ONLY*) 

WR ITE( 6, 707 0) I FIR ST, ILASTC 
7070 FORMAT!* *,* INPUT READ F R QM * ,14* * TO *,I4) 

4 WR I TE( 6,75) FACTO, KRUN »E G G F RT (1 ) , EGGFRY <2 )* EGPC1 )♦ EGP (2 ), PL ANTS 

75 FORMAT!* • , 2X , • FA C TD= • , F 6 .3 , 2X * RUN N 0 . - * , I 3 » • Z FA E6GFRT -• ,F6 . 3 , * VI 
XRS EGG^RT-* ,F6, 3,/ ,* ZEA EGGS/MOTH **,F6.0,* VIRS EGGS/MO TH = *, F 6# 0 
X. • PLANTS /ACRE-* ,F6.0, /) 

WRITE! 6, 77 00 > IE, I A, IS, I M 

7700 FORMAT!* •.•CROP PHENOLOGY: IE—* , 13, *,1A=*,I3, f , IF= * ,13, » , IM=*, 13/ 
C ) 

WR ITEC6 ,76 ) 

76 FOPMATC* DAY NO. EGGPAR 1 -3 PAR EGGPRD 1-3PRD TCLOUD NITMIN TeMAVG 
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HEADS SAREAM 


1 68 
1 69 
170 
1 71 
1 72 
1 73 
1 74 
175 
1 76 
1 77 
178 
1 79 
1 80 
1 31 
1 82 
1 83 
1 84 
1 85 
1 86 
1 67 
188 
1 39 
1 90 
1 91 
i 92 
1 9*3 
1 94 
1 95 
1 96 
1 97 
1 98 

1 99 
200 
201 
202 

2 03 
2 04 
205 
2 06 
207 
200 
2 09 
210 
211 
212 

213 

214 
2 15 
2 16 
217 
219 

219 

220 
221 
2 22 
223 
2 24 
225 
2 26 
227 
2H8 
2 29 

230 

231 
2 32 
233 
2 34 
2 35 
236 
2 37 
2 38 
2 39 

240 

241 
2 42 
2 43 

244 

245 

246 
2 47 

248 

249 
2 50 
2 51 
2 52 


2 ADMOON EGLYTM SQUARES BOLLS EARS 
DO 70 I - I F I RS T * 1 LA ST 
WR I TE ( 6# 79) I « ( PAR S I T( I . J ) * J = 1 *2 ) * ( PRED OR ( I . J ) . J— 1 .2 ) ,T CLOUD ( I) , NIT 
2 MIN (I ) , TEMAVG ( X ) ,ADMOON< I ), EGLYTM! 1 > , SQUARE ( I ) ,BOL_ S<t),EARS(I),Hz 
3AD S( I > t SAREA ( I ) 

79 FORMAT < • », I 7 ,4F7 . 3. 3F7. 1.F7.3.F7. 1.4P7,0«F7.1) 

70 CONTINUE 

DO 72 IN =IH,IV 
WR HE < 6 » 603) 

603 FORMAT!' ' ,// .3 OX , 1 NATURAL MORT AL IT Y • * / , ■ EGGS 1-2 LARVAE 4- 5 L 

XARVAE PUPAE* ,6X, ‘ADULTS* ) 

WRITE (6 .604 ) (DEATH! I . IN) , 1M, 5) 

604 FORMAT! 1 2F 10. 3) 

WR IT E ( 6 ♦ 77 ) 

77 FORMAT!* i DAYNO. FCSPAR POS PRD POS1-3 POS4-5 POSPUP POSADT POSEGG 
2 EGG 1-3 4-5 PUP ADP ADO E-E EGLYPR SGLYCR PRMGRT SEXRTOM 

IF < I N.EQ.l ) W R I T E < 6 .605) 

IF! IN .EG .2) WR ITE( 6,606) 

DO 72 I=IFIRST, ILAST 

WR I Ti ( 6.78)1 . <POS( I *J .IN) ,J=a ,7) » < I GEN! I , J . IN ) , J=1 . 7 ) , EGLYP 

2R( I .IN), EGLYC R( I , IN) , PRMGRT ( I » IN ) . SE XR TO ( 1 . I N ) 

78 FORMAT!* * , I 7 .7 F7 . 3 , 7 14 , 3 F7 . 3 , 2 F7 . 2 ) 

72 CONTINUE 

WR ITE ! 6 ,01 ) 

81 FORMA T ( ' KK I SITES PRBLAR PR BN EW * ) 

720 CONTINUE 

C IFRUI T-ESTI MATEO TO BE 9 DAYS PRIOR TO THE LAST SILKS TO APPEAR. 

C (CORN ONLY) 

IFRU I T = 1 M-9 

IF ( KODE { 1) .EQ.O) KODE ( 1 > -KKOOE 
KKODE = K 3 DE ( 1 ) 

C START OF MAIN PQFU LAT T’CN LOOP. 

DO 200 IN =IH,I V 

C FOODFK-NUTRI ENT VALUE OF CROP CORN LEAVES - 1 * 0-M l LO 1.1/ CORN EARS .65 

C C OT TON- 1 • 0 FOR ZEA ALL CROPS 1.0 FOR VIRESCEN 

C FOODFK-VARI ABLE USED TO INCREASE OR REDUCE LENGTH OF LARVAE 
C DEVELOPMENT. 

FOOD F K -1 .1 

IF ( KCRQP .EG . 1 ) FOO DFK = 1 . 

IF ( I N • EQ .2 ) FCODFK-1 .0 
I F (K CROP .EQ *3 ) FOODFK^l.O 

C NN- VARI ABLE USED TO FLAG WHEN DAY OF SIMULATION IS IFRUIT OR 
C GREATER (CORN ONLY). 

NN-0 

TSIT0C=0 . 

JFIRST = IFIRST 
JL AST C= ILAST C 
C MAIN POPULATION LOCP. 

DO 200 K = IFIR ST, ILAST 
521 CONTINUE 

C K TRAP-CODE NAME OF TYPE OF INPUT. 

KT RAP=KD DE ( 1 ) 

IF (K0DE<6) * EQ .2 )GO TO 51 
IFtK.GT. ILASTC) GO TO 50 
GO TO <52,520 ,52,520, 8, 520. 520 ) .KTRAP 
GO TO 8 

C POP NO 0-1 CAY OLD EGGS. E QEL AD-NUMO ER OF EQUI VALENT EGG LAYING 
C ADULTS. 

C TOTAL-! K, 7, IN) INPUT DATA CONVERTED TO MOTHS PER ACRE. 

C K-DAY NUMBER OF SIMULATION. IN-SPECIES OF INSECT. 

C STATEMENT NO 51 LSED CNLY WHEN INPUT IS MIGRATING MOTHS. 

51 POP1K-1 )—( EQELAD! K-l, IN) + (TOTAL( K— 1,7 « IN)*EGLYPR( K— 1 .IN) ) ) * EGPEMT ( 
1 I N ) * A D MO C N ( K- 1 ) 

GO TO 0 

C STATEMENT 50 USED ONLY WHEN DAY OF SIMULATION IS PA 5T LAST DAY 
C OF INPUT. 

C GENERATED EQELAD IS USED AS INPUT, 

5 0 POP ( K— 1 ) =EQE LAD (K-l , IN* *EGPEMT < IN) * ADMOON ( K - l ) *E GL YPR ( K- 1 , I N ) 

GO TO 8 

C STATEMENT 52 IS USED WHEN INPUT IS LIGHT TRAP CATCHES. 

52 POP! K-l ) =T0TAL< K-l ,7, IN ) *EGP£MT ( IN ) *EGL Y CR (K-l , IN ) 

GO TO 8 

C STATEMENT 520 IS USED WHEN IN^UT IS EGGS. 

520 POP( K- 1 ) - TOTAL! K-l .7, IN) *EGPEMT (IN) 

C APOSON-STOR AGE OF NUMBER EGGS STERILIZHD FROM INSECTICIDE/ 

C APARA-STORAGE OF NUMBER OF EGGS PARASITIZED. 

8 CONTINUE 
APOSON=0 
APAR A — 0 
TOTPAR=0 
TOTPOS-O 

IF ( POP( K-l ) • L T* 1. ) GO TOIIOO 

C XBAR-DAY NUMBER THAT THE MEAN NUMBER OF EGGS HATCHED. XSD- SDE V- ( K , 1 > 
C STANOARD DEVIATION OF EGGS HATCHED DUE TO TEMPERATURE# 

C OISTS- THE PERCENTAGE OF EGGS HATCHED. 
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253 
25 4 
2 55 

256 

257 
2 58 
2 59 
260 
261 
262 
263 
2 64 
2 65 
266 

267 

268 

269 

270 
2 71 
2 72 
2 73 

274 

275 
2 76 
277 

2 78 

279 

280 
2 01 
2 Q2 
283 
2 84 

285 

286 

267 

268 

289 

290 

291 

292 

293 
2 94 

295 

296 

297 

298 

2 99 

300 

301 

302 

303 

3 04 

305 

306 

307 

308 
3 09 
310 
3 11 
3 12 

313 

314 
3*5 

316 

317 

318 

319 

320 

321 

322 

323 
3 24 
3 25 
326 
3 27 
3 28 
3 29 
3 30 
331 
3 32 
3 33 
3 34 
3 35 
3 36 
337 


XBAR-K+I GEN< K *1 *IN)-1 
XSD=SQ£ V( K*1 ) 

D 1 STS =0 . 

C SRVIVE-NU OF FERTILE EGGS REMAINING IN THE POPULATION AFTER ALL 
C MOR TALI TVS. 

C AP— THE PERCENTAGE OF EGGS NOT HATCHED. 

SR VI VE~POP< K-l ) * EGGFRT (IN) 


AP= 1 • 

C BEGINNING OF EGG LCOP. 

DO 1 00 I ”K .1 LAST 
IF (SRVIVE.LT • 1 • ) GO TO 1100 
OPREDORd *1 > *(1.-P05C 1.2 . IN)) 

C T T= THE PROBABILITY THAT EGGS WILL SURVIVE NATURAL MORTALITIES AMD 
C PREDATION. TP=THE NO. OF EGGS P AR A3 1 T I 5 E D- AS SUM ES IF DVICID5 APPLIED 
C THAT IT WOULD ALSO KILL PARASITE EGG. 

C TPP- THE NO. OF EGGS STERILIZED BY OVICIDE 
C PARASITISED EGGS ARE ASSUMED TO REMAIN IN FIELD 10 DAYS 
C INFERTILE EGGS ARE ASSUMED TO REMAIN IN FIELD FOP 3-4 CAYS 
C TINF-TQTAL INFERTILE EGGS IN POPULATION. 

10 TT = ( 1 * — C ) * ( 1*— DEATH (1 .IN)) 

POP< I )=POP(I- 1 ) *TT 

IF < KODE( 5) . NE .4 ) GO TO 1120 

PT EM P“ A6S( 72 # — TEMA VG( I ) ) 

IF ( PTEMP .GT .2 2. ) PT EMP=22 
IF ( PTEMP.LT. 0.) PTE MP=?S 
PT EMP - 1 . 56- .0 7* PTEMP 
8=-. 693/(15O00.*SAREA< I ) ) 

0= (1 .-EXP (PA RSI T( I , 1 > *3 ) )*PTEMP*( 1 . - PG S C I . 1 , I N ) ) 

TP = (POP( n-AP/^RA) *B*< 1 . — POS (1*7. IN) ) 

CALL TRKGRM 
GO TO 1121 
1120 CONTINUE 

B=PARSIT (I . 1) *< 1. -POS ( I . 1 .IN)) 

TP= ( POP ( I >-APAR A) *B*( 1 .-POS ( I«7 * IN ) ) 


) * (1 . -P0S< I . 7 t IN ) ) 


CALL TRKGRM 
1121 CONTINUE 

SR VI VE=SRVIVE*TT*( 1 .-B 
APARA=< A PARA* TT)+ TP 
TPP-POP ( I ) *POS< 1 . 7 . IN ) 

APOSON = ( APQSON* TT) +TPP 

T I NF= ( P3 D ( I >- AP AR A > *( l.-SGGFRT( IN) J + APOSON 
: TOT AL-< I .8 . IN)TOTAL NUMBER OF EGGS IN POPULATION ON DA Y I . 

TO TAL ( I ♦ 8.1 N) “T CTA L( I .8 1 1 N) +TI NF 
IF( AP.LT . .000 1 ) GO TO 1110 
X - 1 

: ZZ-NO OF EGGS THAT HATCHED ON DAY OF SIMULATION. 

ZZ=SRV IVE*DISTR (DISTS )/AP 
GO TO 1111 

1110 ZZ = SR VIVE 

1111 CONTINUE 

AP =1 • -0 ISTS 
POP( I )-POP( I )-ZZ 

; SLPOP-NQ. OF 1 DAY OLD SMALL LARVAE* 

SL POP ( 1 + 1 >“SLP0P< I Fi >+ZZ 
SRVI VE=SRVIVE-ZZ ^ 

C TO TA L- ( I ♦ 1 • I N) T OT AL NO. OF EGGS IN POPULATION 
I INCLUDING INFERTILE AND PARISITIZED EGGS. 

TOTAL (I » 1 « I N ) =TCTALC 1.1. I N ) + POP < I ) —A PAR A+T P 
100 CONTINUE 
1100 CONTINUE 

; CHECK TO SEE IF LARVAE IS DEVELOPING ON EARS. 

IF(KODE{ 11 ) .EQ. 1 .AND. K.GE.IFRUIT) F3DDFK-.65 
SMALL LARVAE ARE EXPOSED TO PARASITES AND PREDATORS 
2DAYS WHEN SILKS ARE AVAILABLE FOR SITES 

IF (KODE(li). EQ .1 ♦ AND . K.GE.IFRUIT) NN-1 
PP- TEMPORARY STORAGE GF SMALL LAEVAE 0 TO l 
LX-CALCULAT ION OF APPROX. DAYS REQUIRED r OR 

larvae. 

PP=SLPOP(K) 

. LX-IGENC K,2. I N > *FOODFK+ .5 


ON DAY I. PARENT HIS 


FOR ONLY 


DAY OLD ON DAY K. 

development o- small 


c 

c 

c 

c 


c 

c 


L-K 

IF (SLPOP (K ) .L T. I . )GO TO 1201 
XBAR- MEANS DAY NUMBER WHEN SMALL 


LARVAE BECOMES LARGE LARVAE. 


XSD-SDE V( K. 2> STA NDARD DEVIATION FOR SMALL 
D ISTS -PERCENT AGE OF SMALL LARVAE THAT HAVE 
LARVAE. 

XBAR=K4*L x- 1 
XSD=SDEV (K . 2 > 

D I ST S~ 0# 

SR V I VE=SLPOP ( K) 

SRVI VE-TEMPORARY STORAGE OF SMALL LARVAE SURVIVING. 
A P— 1 • 

AP-PERCENT AGE OF SMALL LARVAE THAT 
LAR VAE. 

N=0 


LARVAF DEVELOPMENT. 
DEVELOPED INTO LARGE 


HAVE NOT DEVELOPED IN TO LARGE 
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338 
3 39 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 
3 54 
355 
3 56 
357 
3 58 
359 
3 60 
3 61 

362 

363 
164 

365 

366 

367 
36a 

369 

370 

371 

372 

373 
3 74 

375 

376 

377 
3 78 

379 

380 
3 81 
362 
383 
3 04 
365 
3 86 
3 67 
388 
3 69 
3 90 
391 
5 92 
393 
3 94 
5 95 
3 96 

3 97 

398 

399 

4 00 
401 
4 02 

403 

404 
4 05 
4 06 
407 
4 0 8 

409 

410 

411 

412 

413 
4 14 
415 
4 16 

417 

418 

419 
42Q 
421 
4 22 


C N— AGE QF SMALL LARVAE. 

XPOP =0. 

C XPOP-TGTAL NO. OF SMALL LARVAE CANNABILlZED BV LARGER LARVA-. 
CANLAP=0. 

C C AN LA R— NO . OF SMALL LARVAE C ANN AB ALL ZED BY OTHER LARVAE OP THE SAME 
C SIZE. 

C START OF SMALL LARVAE LOOP. 

DO 101 i=k,tlast 

IF ( SPVI VE.LT* 1 « )G0 TO 1200 

N- N + 1 

C CHECK TO SEE IF SMALL LARVAE ARE ON CORN SILKS. 

IF ( NN.EQ. t.AND. N. GT.2) GO TO 14 
D =PARS IT ( U 2> *< 1 ,-POS(I t 1. IN> ) 

C EQUATION WHICH REDUCES PREDATION PERCENT LEVEL DUE TO A GEL OF SMALL 
C LARVAE. 

B ” ( E X P < — *42 *N ) / .657 ) *PREDQR ( I, 2) *( 1 *-P3 5( I ♦ 2. IN ) ) 

GO TO 1 3 
14 0=0* 

E-0. 

13 T - { 1 • - D ) * ( 1 .-E ) * ( 1 . - DBA TH (2.1 N) ) *< 1 .-PQS ( I ♦ 3 * I N ) > 

SRVIVE=SRVIVE*T 
KK-I 

C JJ-CODE TQ OPERATE CANNIBALISM SUBROUTINE. 

J J=365 

C CANNIBAL- SUBRQUTI NE THAT CALCULATES CANNIBALISM FOR 7 E A LARVA^. 
CALL CANBAL 
1301 X = I 

IFtAP *LT . .00 0 1 ) GO TO 12 10 

C ZZ“ NUMBER OF SMALL LARVAE THAT HAVE DEV ELOPED INTO LARGE LARVAE ON 
C DAY OF SIMULATION. 

ZZ-SPVI VE*DI STR {DISTS )/A P 
GO TO 1211 

1210 ZZ=SRVIVE 

1211 CONTINUE 

AP =1. -DISTS 

C LLP OP— NO OF LARGE LARVAE 0-1 DAY OLD. 

LLFOP( 1*1) =LL PQP( I +1 ) -f ZZ 
SRVI VErSRVl VE“Z Z 


5 ( 1 T0TAL NUMBER OF SVALL LARVAE IN THE POPULATION ON 

C DAY I i 

TOTAL < I .2 .1 M=TOT AL< 1.2 , IN J +SRVIVE 
101 CONTINUE 

1200 CONTINUE 

1201 CONTINUE 

C LLP OP —l AP G 5 LARVAE POPULATION (0-1 DAY OF AGE). 

SR VI VE=LLPOPf K> 

L=K- LX+1 

_ IFIL.LT. IFIRST) L=K 

C L ESTIMATES DAY WHEN EGGS WERE LAID ON FRESH SILKS FOR CA'JADOLISM 
C OF LARGE LARVA. 

J J — O 

CALL CANBAL 

WR ITEC6 . 600 )L .KK.R »PRQuAR •PR6NE W 
600 FORMAT (2 17 .FI 0.0.2 F7 « 3 ) 

IF ( LLPDP ( K) .L T# l)GO TO 13 00 
C I GE hH (K i 3 » I N) DEVELOPMENT TIME FOR LARGE LARVAE, 
r ** E FCR LARGF LARVAE ADJUSTED FOR NUTRITION. 

C XBA R— OAY NO-. WHEN THE MEAN NUMBER OF L43GE LARVAE WILL PUPATE. 

C XSD-STANDARD DEVIATION FOR PUPATION. * 

£ ? l ST 5™ E RC “ NT 0F LARGE LARVAE THAT HAVE PUPATED. 

C AP- PERCENT QF LARGE LARVAE THAT HAVE NOT PUPATED. 

LX-IGEN< K « 3 « I M *FOCDFK+.5 
XBAR = K +L X- 1 
XSD=SDEV<K ,3 ) 

Dl S T 5—0. 

AP — 1 . 

: START OF LARGE LARVAE LOOP. 

DO 102 I=K,ILAST 

IF (SRVIVE.LT. 1 . )GO TO 1300 

T= ( 1 # — 05 ATH ( 3 .1 N) > * ( 1 *-PQS( I ,4 ♦ I N) > 

: T-MORTALITY DUB TO INSECTICIDE AND NATURAL MORTALITY. 

SRVIVE-NO. GF LARGE LARVAE SURVIVING* 

SR VI VE-SRVI VE 4T 
X= I 

IF ( AP .LT*. 0001 ) GO TO 1310 
ZZ = SRV IV E*OISTR< DI STS )/AP 

Z Z— NO OF LARGE LARVAE THAT HAVE PUPATED ON DAY OF SIMJLATION. 

GO TO 1311 

1310 ZZ — S R V I V E 

1311 CONTINUE 


AP -1 • — D ISTS 
C PFQF-NO.OF PUPAE* 

PP OP< 1+1 J =PPO P ( 1+ 1 ) + ZZ 
SR V I VE =S RV IVE — ZZ 
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423 
424 * 
425 
4 26 
427 
42 ^ 
4 29 
4 30 
431 
4 32 
433 
4 34 

435 

436 

437 
4 39 
4 39 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 
4 SI 

452 

453 
4f 4 
4 55 
4 56 

457 

458 
4 59 

460 

461 

462 

463 

464 

465 

466 

467 

466 
469 
4 70 
471 
4 72 
473 
4 74 
4 75 

476 

477 
4 78 
4 79 
4 80 
481 
4fl2 
4 83 
4 84 
485 
4 66 

467 
480 
489 
4 90 

491 

492 

493 

494 
4 95 

4 96 

497 

498 

499 

500 

5 01 
5 C 2 
5 03 

504 

505 

506 
5 07 


C TOTAL— ( I 1 3 • X N ) TOTAL NO* OF LARGE LARVAE IN THE POPULATION* ALL AGES 
C ON CAY OF SIMULATION* 

TOT AL ( I < 3 » I N ) =T OT AL< 1*3* IN) +*5RV I V E 
102 CONTINUE 

C SRV I V E— NO ♦ OF PUPAE (0-1 DAY OLD) 

1300 SR VI VE=PPOP< K > 

IF (PPOP( K) .LT .1 >GO TO 1400 

C IXBAR-X8AR- DAY NO. WHEN THE MEAN NUMBER OF PUPAE BECOME ADULTS. 

C XSD-STANDARD DEVIATION OF EMERGING ADULTS. 

C DlSTS-PERCENT OF PUPAE THAT HAVE EMERGED. 

C AP- PERCENT OF PUPAE THAT HAVE NOT EMERGED. 

I X BAR=K + I GEN ( K, 4. I N)- 1 
XBA R-IXBAR 
XSD = SDEV ( K .4) 

D I STS — 0 . 

A P~1 ♦ 

C START OF PUPAE LOOP. 

00 1 03 l -K .1 LAST 

IF{ SR VI VE. LT* 1. >G0 TO 1400 

C SRV IV E— NO . OF PUPAE SURVIVING M OR TAL I TYS . < INSEC TI C I DE + NATURAL) 

SR VI VE^SRVI VE*( 1 .-DEATH < 4 , IN )) * ( 1 .-POS < I • 5t IN) ) 


□ N DAY OF SIMULATION. 


■1 DAY OF AGE ) . 


X=I 

IF (AP.LT . .0001 ) GO TO 1410 
C ZZ- NO. OF PUPAE THAT HAVE EMERGED 
ZZ-SRVIVE*DISTR(D I STS >/AP 
GO TO 1411 
141 0 ZZ = SR VI VE 
1411 CONT INUE 

AP -l.-DISTS 
SRVI VE=SRVI VE-Z Z 

C POAPOP-NO. OF PRE-OVIPOS IT IONING ADULTS. < Q- 
POAPOP ( I +1 > “PDA POP ( I + 1M-ZZ 
TOTAL( I. 4. IN)=TOTAL(I *4. IN) +SR VI VE 
C TOTAL-! I . 4f IN)TOTAL NG. OF PRE OV I PQS I T I 0 N I NG ADULTS A- L AGES IN 
C POPULATION. 

103 CONTINUE 
1 400 CONTINUE 

C JJ-CAY NO. WHEN THE MEAN NUMBER OF PRE-OVIPOSI TIONED MDTHS BECOME 
C OVI POSI TI ONED ADULTS. 

JJssK+I GEN( KfSf IN) 

IF (POAPOP<K) .LT .1 . IGO TO 20 2 
IF! JJ.GT.l LAST) JJ^ILAST 

C SRVIVE-NO.OF PR E- O V I P OS I T I ON I N G ADULTS! 0-1 DAY OF AGE). 

C L- AGE COUNTER OF ADULTS. 

C DEA THY-MGRTALI TY OF ADULTS FROM PREVIOUS DAY. 

SRVIVE^POAPOP <K > 

L = 1 

CEATHY=0 .0 

C ADTS- NO. OF ADULTS 
C OF ADULTS THAT HAVE 
C THIS PARTICULAR CHORT . 

ADTS! K -1 ,1 N) =0.0 

C START OF PR E-OVAPOSI TIONING LOOP. 

DO 1 04 I -K, J J 

C ADULT- SUB RO UTI NE USED FOR CALCULATIONS FOR NATURAL MORTALITY DUE TO 
C TEMPERATURE. 

C SURV-PERCENT SURVIVED TO THE NEXT DAY. 

CALL ADULT 

SRV I VE — S RV I VE *SURV *! 1 .-DEATH! 5. I N ) ) *< 1 .-PDS! I . 6* IN) ) 

IF!KODE< 6) *£Q *0 )GO TO 1041 
IF(KODE( 11 ) .EQ. 3) GO TO 1041 

C THIS REDUCES THE NUMBER OF ADULTS IN CROP EACH DAY BY THE PERCENT 
C PR0BA8LITY THAT THEY HAVE MIGRATED AWAY. 

ADTS ! 1 • I N ) = ADTS ( 1-1 t IN ) *5URV *( 1. -DEATH! 5 • I N ) > f SR VI V£ 

1 R T ( I .1 N) 

C ADTSMG-TOTAL NO. OF ADULTS THAT HAVE MIGRATED FROM A PARTICULAR CROP 
C ON DAY OF SIMULA TI CN. ( ALL AGES) 

ADTSMG! I .KCRQPf IN>=ADTSMGU >KCROP»IN) + ADTS! I *IN) 

SR VI VE-SRVI VE*< 1 - P R MG RT ! I • I N ) ) 

CONTINUE 

TOT AL-( I f5 f IN ) TOTAL NO OF PRE - OV IP05IT IONED ADULTS ON DAY OF 
C SIMULATION. (ALL AGES) 

TOTAL ( I * 5. IN ) -TOTAL! I * 5 * I N ) ♦ SR VI VE 

104 L-L+l 
LL-L-1 

C T-TEMPERATU RE SUM OF RUNNING AVERAGE. 

C AGE-LENGTH OF TI ME MOTH IS BEING OV I P05 I T I ONED . 


M IGR AT ING ON 
MIGRATED AND 


DAY O- SIMULATION 
SURVIVED FROM THE 


PLUS THE 
PREV IOUS 


NUMBER 
DAY FOR 


* PRMG 


1041 


T=0. 

AGE=0 . 

KK - J J+ 1 

IF(KK .GT .ILASDGQ TO 202 
C JJ-MEAN LENGTH OF CV I POS I T I ONI N G PERIOD. 
JJ = JJ + I3EN!K, 6. IN) 

1F( JJ .GT . I LAST ) JJ-ILAST 
C START OF OVAPOSI TI CNED LOOP. 
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5 0a 
5C9 

510 

511 

512 
5 13 

514 

515 

516 

517 
5 16 

519 

520 

521 
5 22 

523 

524 

525 
5 26 

527 

528 

529 

530 

531 
5 32 
533 
5 34 
53 5* 
5 36 
5 37 
53Q 
5 39 
540 
5 41 

542 

543 

544 
5 45 

546 

547 

548 

549 

550 

551 

552 

553 

554 

555 
5 56 
5 57 

558 

559 

560 

561 

562 

563 

564 

565 

566 

567 

568 
5 69 
5 70 
571 
5 72 

573 

574 
5 75 
576 
5 77 
578 
5 79 
5 Q0 
581 
5 82 

583 

584 
5 85 
5 86 

587 

588 
5 89 
5 90 

591 

592 


DO 105 I=KK«JJ 
AGE=AGE+1 .0 
T-T+T5M4 VG ( I ) 

C T RU N— RU NNI NG DAILY AVERAGE TEMPERATURE# 

TRUN^T /AGE 

IF (TRUN.LT .55 .) GO TO 1052 

C IPEAK-NO. OF DAYS IT TAKES AN OV I POS I T I ONI N G MOTH TO REACH MAX. EGG 
C PRODUCTION. 

IP E AK= < 3 0 2920 2. 9 *( TRUN ) * +! - 3. 1 11 )> + • 4 
IF(L.LT# IPEAKJGC TC 1053 

C I TOP- NO . O p DAYS MOTH WILL 9E AT MAX PRODUCTION OF EGGS. 

IT OP- ( 2 3 9741 • 3*< TRUN ) * * ( -2 • 61 540 2 ) >> . 5 

I PEAK— I PEA K+ I TOP 

IF(L .LE. IPEAK )GO TO 1054 

C I END- NO. OF DAYS AN OV IPOS IT ION ING MOTH WILL PRODUCE EGGS. 

IEND = ( 91 C99 8. 7* ( TR LN ) * * ( -2 . 53 ) ) * • 5 
IFCL.GE. IEND ) GO TO 10 52 
XS=L-I PEAK 
XB= IEND- [PEAK 

C ELPOTA-PROB ABIL ITY THAT A MOTH IS AT MAX. EGG PRODUCTION. 

ELP3TA=1 .-XS/X0 
GO TO 10 60 
1 052 ELPDTA-0 . 0 
GO TO 1060 
1 053 XS=L-LL 

XF=IPEAK-LL 
ELPDTA=XS/XF 
GO TO 1060 
1054 LLPDTA-1 . 0 
GO TO 1060 

C REDUCES PROBABILITY OF MAX PRODUCTION OF EGGS EFFECTED 8Y MOONLIGHT. 
1060 £LPOTA=ELPDTA 
CALL ADULT 

SR VI V£ - SR VI VE * S UR V * ( 1 .-DEATH (5 .IN) > * U .-POS { 1 .5 . INI) 

IF ! KD DE! 6) .EQ .0 ) GO TO 10 51 
IF (KCROP.GT .2 ) - GO TO 1051 

ADTSI I .1N)=ADTS( I - 1 *IN) ♦ SUR V * { 1 . -DE ATH (5 . I N ) ) +S R V I V E + PRMGRT ( 
2 I , IN) 

ADTSMG ( I § KCROPt I N) =ADTSMG(I * KCROP# IN) +ADT5 < I. IN ) 

C TOTMGS-TDTAL NO. OF 0 V IPOSI T I ONI NG ADULTS THAT HAVE MIGRATED AND 
C SURVIVED. 

TQTMGSU . I N ) ~ TO TMGS( I ,IN)+ADTS( I ♦INI *E'LPDTA 
SR VI VE^SRVI VE*( 1 • - FRMGRT ( U IN) > 

1051 CONTINUE 

C EQELAD-EQUI V IL ANT EGG LAYING ADULTS < TOTAL NO* OF OVl PO SI T I ONI NG 
C ADULTS. 

C REMAINING IN CROP* AD J LSTEO FOR EGG LAY POTENTIAL AND MOONLIGHT) 

EGELAD { I * I N ) = E QEL AC ( I * IN ) +5 RV I V E *ELP DT A 
C TOTAL-! I *6 * IN) TO TA L NC. OF 0 VI P CS I T I ONI NG ADULTS IN POPULAT I ON - 
C ALL AGES ON DAY OF SIMULATION. 

TOTAL! I * 6*IN)-TCTAL( I *6 * I N)+SRV IVE 
105 L=L+1 


202 IF <K0DE<7 ) .NE .K > GC TO 200 

605 FORMAT!' * .2 5 X . • B 0 LLW ORM HELI0THI5 ZEA* ) 

606 FORMAT ( * *• 25X • TOBACCO BUDWORM HELIOTHIS VIRESCENSM 

200 CONTINUE 

c PRINTOUT OF ALL POPULATION TOTALS. 

DO 1070 IN = I H * IV 
WRI TE( 6.608) 

608 FORMAT (• 1 * .25X, 'PR INTQUT OF TOTAL POPULATION OF VARIOUS LIFE STAGE 
* S ON ANY DAY' ) 

IFON.EQ.2) GO TO 1071 
WR IT E (6 * 605 ) 

GO TO 1072 
1071 WR I T E ( 6 * 60 6) 

1 072 WR I TE (6 *610) 

610 FORMAT ( * DAY NO. EGGS 1 5T-3RD 4TH-5TH PUPAE P ADULT O AOULT 
EQELAD MIGRATS IN EGGS PAR EGG TOTMGS TRICHS*) 

IF< KCROP *G T. 2 I KCR C3P=2 
FACT B8-KRUN- ! KCROP+10.)-1 . 

FACT AA -F ACTS B/ ( FACT BB4-1 . ) 

FACTBB^F AC TBB+1. 

00 107 K-I FIRST . IL AST 
IFCKRUN.LT.il) GO TO 1076 

^ T ^ G ^4 KCR0Pt IN )=ADTSMG !K . KCROP .IN) * F ACTA A+TOTMGS < K , I N > /F ACTBB 
GU TO 1077 

1 076 ADTSMG! K .KCROP * I N) = TOTMGS (K . IN) 

1077 CONTINUE 

IF ( KODE( 16) . GE .200 ) GO TO 1073 

Im ) oV T .T°J AL( )S( INI *TOTAL(K.l,INJ 

}£<f { I .LT .TCT AL(K t7 > IN) ) S ( IN + 2 >= TOT AL ( K , 7, IN ) 

N) * S( IN+2 )=EQELAD(K, IN) 

107 3 CONTI NUE S * LT * ACT SMG< X . KCROP . IN ) ) S( I N +2 ) =ADTSMG< K > KCROP » IN) 

T0T9 -( TOTAL! K,9 t IN) /( TOTAL! K, 1 . I N ) +1 . ) >*100. 
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WRITE(6,602)K, ( TOT AL ( K, J , IN) , J = 1 , 7)»EQ£LAD(K»IN) , ADT SMG{ K , KCFOP, I N 
1 ) • IOTA L< K.8.U) *TQT9*T0TMGS(K»IN) ♦ PA RS I T { K » 1 ) 

60 2 FORMAT ( « ' , I 6 * 1 2F 8 . 0, F 0 . 1 , F 8 . 0 ) 

TOTAL ( K » 8 ♦ I N) CTS M G< K , K CRQ P » I N ) 

107 CONTINUE 

IFUN.EQ.2) GO TO 1070 
I P X” 1 3 

C 5 MO TH-5UBR0 UT I N F FOR APPLYING 3 DAY MOVING AVERAGE TO ADULTS THAT 
C HAVE MIGRATED. 

ILASTC=I LAST 
CALL SMOOTH 
1070 CONTINUE 

IF (KDDE(16) • GE *20 0)G0 TO 1075 
C SETS FLAGS FOR PLOT OF DATA. 

IF ( KOD t ( 16).GE.100) NPL OT -K 0 DE ( 1 6 ) -l 00 
IF(KODE< 16) ,LT. 100 ) NPL3T =KODE(16) 

DO 1 074 1=1 , NFLCT 

RE AD (5,510 ) ( ( IGKODEI J ,K ) , J= 1 .1 2 ) .K = l ,3) 

510 FORMAT (3612) 

K1=U 
K2=l 2 

DO 1082 K= 1 ,3 
DO 108 1 J=KWK2 

IF< IGKODEI J, K> . EG.O ) GO TO 1081 
IF ( IDX.EQ. 0) GO TO 552 

DO 1083 KJ=IF IRST, ILAST 

C DAT X-ARR AY FOR STORING FIELD DATA OF INSECT POPULATION WANTED ON 
C PLOT TO COMPARE TO SIMULATED WEATHER- 
1083 DATX(KJ) =0 . 

IDX-1 

552 RE AD < 18, 551 1 END= 1 081 ) ID, X 
551 FORMAT ( 6 X » I 6 ,F6 ,0) 

I F ( ID. EQ * 0 ) GO TO 1081 
0 ATX ( ID) =X 
GO TO 55 2 
1081 CONTINUE 

C ^GRAPH-SUBROUTINE FOR PLOTTING VARIOUS STAGES FOR INSECT POPULATION. 
CALL GRAPH 
1074 CONTINUE 

C^EcTJ^sIe IF THIS IS LAST SIMULATION WITH INSECT POPULATION. 

RE AD ( 5*550, END=108E) KRUN 
550 FORMAT (I 3) 

IF (KRUN. EQ.O) STOP 
DO 1080 1=1,4 

1 080 S(I)~0.0 
GO TO 34 
1085 STOP 
END 


C THI !*"* SUBROUTINE ° CALCULATES MORTALITY RATES FOR ADULTS AS DESCRIBED 

C IN COMMON K0DE( 20 > .DUM60 (1 I 788 ) .TEMAVGI 370 >,DUM61< 16282) >L , I .DEATHY. 
A SURV. SAD ULT< 228 6) 

IFITEMAV G( I ) .LT .68 . ) GO TO 103 
F ACT K=— 1 .291 9 + . 02 7 2 4 * TEM A VG < 1) 

GO TO 102 
103 FACTK=.56 

C FACTK=t. 22-92, 1.00-86. .B06-77. .560-66 

102 OEAT HX=1 .7724+EXP< -11 . 44 57/ ( P ACTK*L) ) 

IF(L.EQ.l) GO TO 101 

DEATH Y= 1 . 7724*E XP ( - 1 1 . 445 7/ ( F ACTK* ( L-l ) > ) 

101 SURV =( 1 .-DEATHX )/( 1 .-DEATHYJ 

RETURN 
END 


C T Hi| U SUBROJt1nE A CALCULATES THE PROBABILITY OF CANNIBALISM OF H. ZEA 

° ^COMMON 0UM48C 32B1) ,EARS( 370) ,0 UM4 9 ( l 0747 ) . T OTAL ( 370 . 9 .2 ) . DU M50 (740 
A2) «KK. DUMB (74 4 ) .KCRUP, DUM51 , SITES! 369) , S.DF , DUM5 2( 74 0) .CANLAR.JJ.I 
B .NN.XPOP ,N,PP.DUM54(2) ,P OP , S CA NB 1 425 ) . P RBLAR » PRBNEW 
C NN EQUAL l WHEN CROP IS CORN AND SILKS ARE PRESENT 
S= SI TES ( KK) 

IF(NN.EQ.1)S=EARS<KK ) -E AR S ( KK- 1 5 ) +L . 

C ESTIMATrON e OF°PROEleiLITY QP A LARGER(4-5 INSTAR) LARVAE BEING PRESENT 
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12 

C ON 

A FEEDING SITE ON DAY I 

13 


PRBLAR-TCTALC 1,3,1 )/S 

14 


IF ( PRBL4 R.GT. 1, )PR BLAR -1 

15 


AT E= ,05+ .05*N 

15 


ATEUP-PRBLAR* FOP* A TE 

17 


POP-POP- ATEUP 

IB 


XPQP=XPOP+ATEUP 

19 


IF ( I .LT. JJ) RETURN 

20 

200 

I F < PQ P ,LT . 1 . ) GO TO 1 

21 


F-PDP 

2 2 


*8AR=POP /S 

23 


OF= ,36 

24 


PRBNEW-1 . — { OF / ( DF + XBAR) ) 

25 


C A NLA R -P OP- S* PR BNE W 

26 


POP— POP — C ANL A R 

27 

4 

CONTINUE 

2B 

C ASSUMES THAT ONLY 1 LARGE LAI 

29 


PR8NE W-C ANLAR/P 

30 


GO TO 5 

31 

1 

PR BN EW — 0 • 

32 

5 

IF(PP.LT#1 *) GO TO 6 

33 


PR B LAR ~ XPOP/P P 

34 


GO TO 2 

35 

6 

PRBLAR-0. 

35 

2 

CONTINUE 

37 


RETURN 

28 


END 


WILL SURVIVE PER SITE 


1 

2 

3 

4 
k 

6 

7 

Q 

9 

10 

n 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 
22 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


SUBROUTINE CROP 

Scj 1 fT?w? R ?yi r NE CALCULATES THE VARIOUS PHENOLOGICAL STAGES, THE 
RELATIVE SEARCH AREA, AND NUMBER 3 F SITES (E.G.,EARS)F OR CORN. 0 EG RE E 
-DAYS AS Or, SCR I B ED ARE USED TO MAKF THFSF r At nil ATTnw<; . 



*A22/ — 40. 6667, 83. 3333, -41. 6667/ 

DATA A 30, A 31 , A 3 2 i A 33/ .2.1 .522, —2 . 8977 , 1 .4394/ 

d d a at t£ ^io*. amIi^mII;^?^ :!; • 13822 ” 1 

KcJoP=KS^lS*r 2/4 * 3433l, - 5 * 4l664 ‘ 1 - 6927/ 

IF ( KCROP ,E0, 1 ) DD = 1 24 0 
IFCKCROP .EG. 2 ) DD= 1 390. 

CD A=20 5. 

IF (KCROP .NE. 3) GO TO 66 
CD A=63 1 • 

DDB “57 0. 

COM= 1130 . 

VI ELD-1 , 5 

A5 SU MS 44X |_0S S OF SQUARES AND 1 Q\ LOSS OF BOLLS 
Y 1 tLD- 1*5 

YI ELD S = YIELD* 7 0 00000. 

YIELD8=Y IELD* 76 00000, 

YIELDT-^ 00000 *Y IELO* .8 
CONTINUE 

IF (KRUN, EO.O ) GO TO 9 
IF( KOD 5 ( 14) .E 0. 0) RETURN 

rp d for 5i?j??2!^:^ , il^S:^ ANTS,NATRAt 

FOR L SORGHS M “^l"DElRle E DAVS : AF?E^7E IE MlxtMIN FaR 3 0RG ”U M-N AT R AC < I A > 

I^HF* -AGE-1735 CEGR EE DAVS 
ASSUME 1 05,000 BOLLS/ACRE=lBALE/^RE AC 

g®rflin! T iSWg?. , JI L |S^S 8 (!I? B r 0 S L^ S 5 D 7^° F ° R C0T ™”* 

«tp t Sc 

I iF( ( 7p?E 0 o p :o) 0 ;p=7T ,IP ‘ EO * 0,IP=go 

IF(IE.fcQ.0IIE=Ip«.8 

IF ( NA TRAC. NE . 0) GO TO 60 

IA= IE+1 

IF ( KCROP ,EQ . 3 (GO TO 65 
IF(IF.NE.O) GC TO 63 
SUMrEDCORN C I A )-EDCORN ( I E > 

IF(SUM.GE.CDA)GC TC 62 
l A=I A+ 1 
GO TO 61 


66 


8 

500 


:.<PCTCRP(I> ,1=1 ,5) 


61 
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55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 
63 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
1 00 
1 01 
102 
1 03 
1 04 
1 05 
1 C6 
1 07 
1 00 
1 C9 
1 10 
1 11 
112 
113 
1 14 
1 15 
1 l£ 
1 17 
1 10 
119 
1 20 
X 21 
1 22 
1 23 
1 24 
1 25 
126 
1 27 
1 20 
129 
1 30 
131 
1 32 
1 33 
1 34 
1 35 
1 36 
1 37 
l 38 
i 39 


65 SUMMED l I A ) -L D (I E) 

IF<SUM #GE .DDB )G0 TO 62 
I A — 1 A+l 
GO TO 65 
63 I A — I F— 50 

DDF^DD*. 84 

67 IF ( EDCOR N( IF ) — EDCQRN< IA ) .GE • DDF I GO TO 62 

I A — I A— 1 
GO TO 67 
62 NATRAC— 1 A 
60 KFF = 1FIRST-1 

i a-natrac 

SEXRTO < KFF .1 )=0.0 
SEXRTO < KF F ♦ 2 ) - 0 • 0 
EGLYTM(KFF)-EGLYTM u FIR ST ) 

DO 5 K=KFF,ILAST 

C ADJUSTE5 PROBABLE EGG LAY =QR TEMPERATURE 
X~ A8 5 (75 *-EGLYTM(K ) ) 

I F ( X • G T • 2 1 • ) X- 21 . 

I F ( X .LT *3* ) X-2. 

P-1. 1-. 05146* X 
9 DO 5 IN= IH, I V 

IF (K.GE* NATRAC) GO TO 110 
EGLYPR ( K .1 N> 1 *P 
EGLY CR < K .IN)- .333 
X — 0 • 


GO TO 111 

UC GO TO < 102. 105. 103.100.5 ) .KCR0P 
C NATRAC IS DATE CROP BECOMES ATTRACTIVE TO HELtOTHlS 
C COFN-NATRAC WHEN 6-12 INCHES TALL-1210 DEGREE DAYS 
C FROM NATRAC TO MAXIMUM S ILKING- ( 50-66 DEGREES) 

C NATRAC ALSO DATE WHEN SEARCH AREA-1 UNIT 

102 X- (EOCQRNt K )-EDCORN(NATRAC) )/DD 
IF ( X.GT* 1.60) GO TC 112 

120 IFCK .GT. ILASTC) GO TO 104 
IF< KODE< 1 > .NE .3 ) GO TO 104 

C REPLACE GO TO 104 WITH GO TO 106 IF TRAP SEX RATIO TO 0E USED 
GO TO 10 4 

105 X— { E DC OR N ( K ) -ED CORK ( NAT R AC) )/DD 
IF (X.GT. 1. 2) GO TO 112 
IFOC.GT. ILASTC) GO TO 107 

IF ( KOD L ( 1 ) .NE »3 )GO TO 107 
GO TO 10 6 

103 X— ( ED ( K) -ED< N AT RAC ) )/26 
IF(K .GT. ILASTC) GO TO 111 

IF( KODE< 1 ) .NE .3 .ANC.KODE< 1 > . NE .5) GO TO 111 

106 EGLYCR( K,tN)-(SQRT(SEXRTO(K.IN)))/3. 

IF(EGLYCR(K. IN) .GT .1)EGLYCR<K.IN)-1.0 
IF(EGLYCR(K. I Is) .LT . .166 ) EGL YCR ( K • IN )= .166 
EGLYPR ( K .1 N) = EGLYCR(K , I N> *P 

C— SB XR TO EQUATIONS NEED TO BE PUT HERE OR CAL. OF EGGLAY PROOF- 

104 IF(IH. EQ .1 .AND. IN. EQ. 2)GO TO 2 
IF(X .GT. 1 .12 > GO TO 2016 

IF( X.GT. .B6846) GO TO 2015 

SEXRT O ( K . IN)— A30 F < A 3 1 + ( A 3 2+ A 33* X ) *X ) *X 

GO TO 2019 

2016 SE XRTO(K . IN)-A40 + ( A41+442+X) *X 
GO TO 2019 

SE XRTO( K .1 N > — A2 0+ ( A2 1 + A2 2 * X ) *X 
EGLYCR(K ,IN)-SQRT( SEXRTO(K.IN) ) 

I F ( EGLYC R( KtIK) .LT *.166) EGL YCR (K » IN > = .166 
EGLYPR( K .IN) -EGLYC R< K .IN) *P 
GO TO 111 

IFMH.EQ.l .AND. IN.EQ.21GO TO 2 
IF( X .GT . .035 ) GO TO 2025 

SE XR TO ( K . I N) - AM30 + ( AM31 + { AM 3 2 + AM 33 *X ) *X ) *X 
GO TO 20 26 , 

SEXR TO ( K . IN)- AM20 + ( AM21 *AM2 2*X ) *X 
C REDUCTION OF EGG LAY DUE TO CROP FACTOR (.5 
202 6 EGLYCR (K , I N) = SORT < SEXRTO < K. IN) ) * .5 
EGLYPR ( K • I N) = EG LYC R< K . I N ) *P 
GO TO 111 
EGLYPR (K , IM> = 0 • 

EGLYCR < K .1 N> ==# 1 66 
GO TO ( 10. 11. 12.13).KCROP 
CONTINUE „ 0 

IF( IH.EQ *1 .AND. IN .EQ. 2> GO TO 2 
IF(X.LT . .8 4 ) GO TO 303 
IF ( X.G T. 1. 16) GO TC 303 
I M=K 

TP — ( 1 .-X )*1 6. 66667 
S=.4*EXP<-(TP*TP>/2.) 

OELT AX— ( X— X 1 )/ .06 
ASSUME 1 EAR PER PLANT 


201 5 
2019 


107 


2 02 5 


112 

111 

10 


FOR SORGHUM) 



*S*D£LTAX 


1 40 
1 41 
142 
1 43 
1 44 

145 

146 

147 

148 
1 49 

150 

151 

152 

153 
1 54 
1 55 

156 

157 

158 
1 59 
160 
161 
162 

163 

164 

165 

166 
1 67 
168 
1 69 

170 

171 

172 

173 

174 
1 75 
176 
1 77 

178 

179 
1 80 
181 
182 
1 63 
1 84 
1 65 
186 
187 
1 88 
I 89 
1 90 
191 
1 92 

193 

194 
I 95 
196 
1 97 
l 9 a 

1 99 
200 
2C1 

2 02 

203 

204 

205 

206 
2 07 
2 08 

209 

210 
211 
212 

213 

214 

215 

216 
217 
216 

219 

220 
221 
222 
2 23 
224 


302 EARSCK >= EARS( K- 1 )+PLANTS 
GO TO 304 

303 EARS < K)=EARS(K-1> 

304 XI =X 

IF ( K.LE* NATRAC) GO TO 300 
I F < X *G T • 1) GO TO 301 
SA REA { K ) -X*5 . +1 . 

GO TO 33 

300 SAREA ( K ) =1 .0 
GO TO 33 

3 01 SAREAC K) = SARE A< K-l ) 

C ASSUME 6 SITES PER PLANT AT A SEARCH AREA OF 1 

33 SI TESCK) =SAREA( K) *PLANTS + EARS ( K ) 

IF(KOOc( 6) .EQ .0 )GO TO 5 
CALL MIGRAT 
GO TO 5 

11 CONTINUE 

IFCIH.EQ.1 * AND* IN* EQ. 2 i GO TO 2 
IFCX.LT. .84) GO TO 404 
I F < X . GT . 1 .16 )GO TO 404 
I M =K 

TP— ( 1 .-X >/ .04 

S=.4*EXP(-(TP*TP>/2 • ) 

DELTA X={ X-Xl ) /. 04 

HEADS <K ) -HE ADS C K-l )+*PLANTS*S*DELT AX 
GO TO 405 

404 HEADS(K )=HEADS( K~ 1 ) 

4 05 XI =X 

IFCK.LE. NATRAOGO TO 400 
IF { X .GT ♦ 1 ) GO TO 401 
SAREAC K ) =X*5. +1 
GO TO 43 

400 SARE A ( K )— 1 .0 
GO TO *3 

401 SAREA { K ) =SAR EA< K- 1 ) 

4 3 SI TES ( K) =SAREA( K) * PL ANT S + HE A OS <K ) 

IF< KOD E( 6) .EQ.O)GO TO 5 
CALL MIGRAT 
GO TO 5 

12 CONTINUE 

C I S=F I RST SQUARE. 18 = FIRST BOLLt I E= EMER G ENCE DAT E i IP = D LA NT I NG DATE 
IM=NATRAC+88 

IF< IH.EQ .1 .AND. IN.EQ. 2>G0 TO 2 
IFCK.LT. NATRAOGO TO 21 
IF<X .GT .88* )GO TO 22 
IF ( K. EQ. KFF) X 1=X~1 . 

TP = < 40.- X) /I 3 . 3 
S-.4*EXPC — <TP*TP )/ 2 . ) 

DELTAX-C X-Xl ) n 3.3 
SQUARE ( < ) = S*D EL TAX ^YIELDS 
GO TO 24 

22 SQUARE ( K ) - 0. 0 

24 IFCX.LE.21 . ) G 0 TO 23 
TP^( 66. -X) /I 6 . 

S- .4*EXR { — ( TP *TP ) / 2 . ) 

DE LT AX — ( X— X 1 )/l 6. 

BOLLSC K ) -$*DE LT A X* YlELDB 
GO TO 26 

23 BOLLS ( K) ”0.0 

26 SAREA(K) =X/20*+1.0 

FRU I T-SQUARE < K) +BOLLS <K ) 

2 7 SI TES( K> = SAREAC K) * PLANT S + FRU IT 
XI =X 

IF (K.GT. ILASTC) GO TO 28 

IFCKODEC 1) .EQ .3. OR .KODEC 1 ). EQ. 5) GO TO 5 
C REDUCTION OF EGG LAY DUE TO CROP FACT OT( ■ 7 FOR COTTON) 

28 EGLYCRCK ,IN)-(FRUI T/YlELDT+.l ) *.7 
EGLYPR ( K . IN )=EGLYCR<K . IN ) 

GO TO 5 

21 SQUARE ( K ) = 0. 0 
BOLLS ( K ) =0 .0 
SA RE A ( K) =1 .0 
3ITESCK ) =PLAN 7 S 
XI -X 
GO TO 5 

2 EGLYPR(< « IN ) = EGLYPR ( K » I N - 1 ) 

EGLYCRCK. I N)=EGLYCR(K» IN-1 > 

IF (KODEC 6) .NE. 1 )GO TO 5 
CALL MIGRAT 
GO TO 5 

108 EGLYCRCK • IN >= .5 
EGLYPR (K?IN) — .5+P 
5 CONTINUE 
13 RETURN 
END 
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SUBROUTING EGGINP 

c THIS SUBROUTINE HANDLES THE INPUT OP HELIOTHIS POPULATIONS FOR 
C INITIATING THE MODEL* 

COMMON K CDE (20 ) , F A CT D , D J M 1( 4 74 4 > , IFIRST* ILASTC, IH,IV*DUM2< U 10> * 

A IL A ST * KR UN » SE XR TO ( 370 ,2) .EG PEMT ( 2 ) , PCTC RP< 5 > ,DUM3 ( 5 2 90 ), A Dt S MG 
13(370*2,2 > * TOT AL (370*9*2) • SUB EGG ( 9692 ) 

DIMENSION EQE LA D (370*2) 

DIMENSION S< 1 0) 

I END— 3 70 

IF< KRUN. EQ. 0) GO TC 20 
DO 21 1-1*370 

DO 21 J— I H * I V 
EQELAD < I * J) - 0 • 

21 TOTAL (I*7*J)=0.0 
IF(KODE< 6) .EG. 2 ) GC TO 22 

20 CONTINUE 
L=KODE ( 1 ) 

X=1 . 

I F !L • EQ * 4 ) X= • 3 1 
IP(L*£Q*6) X— • 31 
GO TO ( 1 *2*3* 2* 22, 30*30) , L 
C READS IN MOTHES TOTAL NO. BOTH SEXES COMBINED 

1 DO 19 J=IFlRST*ILA$TC 

U READ ( 1 1* 502) I • BMAL E , 0FMALE, CMAL E • CFM AL E 

I F ( I • NE • J) GO TO 1 1 
TOTAL ! I , 7* 1) = < BMALE+3FM4LE) 

TOTAL! 1*7*2)= ( C M AL E+CFMAL E ) 

19 CONTINUE 

GO TO 5 

C READS IN 0-1 DAY CL D EGGS/ACRE 

2 DO 29 J=IFIRST, ILASTC 

12 READ(1 1, 501 ) I *B,C 

501 FORMAT (9 X, 13 ,2F6. 0 ) 

IF{ I .NE.J ) GO TO 1 2 
IF<L.EQ.4 ) 1=1-2 

TOTAL! 1,7,1) = <B/EGPEMT!1 > )* X 
TOTAL! 1*7*2)= <C/EGPEMT( 2 ) )*X 

29 CONTINUE 
GO TO 5 

C READS IN MOThES SEX SEPARABLY 

3 00 39 J=IFIRST,I LA STC 

13 READ ( 1 1, 50 2) I ♦ BMAL E* 0FM AL E* CMALE * CFM A LE 

502 FORMAT! I 6,4F6.1 ) 

IF(I.NE.J) GO TO 13 

TOTAL! 1*7,1 ) = E3MALE *2. 

TOTAL! 1*7,2) =C MALE *2 . 

SEXRT 0 ! I ,1 )=( EMALE + 1 *0 ) / ( BFMAL.E + 1 * ) 

SE XRTO ! I *2) * < CMALE Pi .0 )/ (CFMALEhl .} 

39 CONTINUE 

GO TO 5 

30 ILASTC=1 
IFIRST = 3 65 
READ(1 1 ,504)8 ,C *1 

504 FORMAT! 1 0X.2F6. G,6X,I 3) 

IF! I « E Q • 0 ) GO TO 5 

IF !L*EQ* 6) 1 = 1-2 

TOTAL! I. 7*1 ) = <0*X)/E3PEMT !I) 

TOTAL! I ,7 *2 >= <C*X )/ EG PEMT (2 ) 

IF( IFIRST.GT* I ) IFI RST=I 
IL=ILAST C 

IF! I LASTC.LT. I) 1 LA STC =1 

IF! IFIRST .EO . ILASTOGO TO 31 

IDD=ILASTC~IL 

IF! I DD ♦ E Q. 1)G 0 TO 31 

XDLl=TOr AL! ILASTC, 7,1 ) -TO TAL ( I L , 7 , U 
XDL2 = TOT AL! I LASTC , 7 ,2 ) -TOTAL ! I L ,7 , 2 > 

XDL 1 = XD*_ i/IDD 
XDL2 =XDL2/ I DD 
IDD= IDD- 1 
DO 32 J=1 , IDD 

TOTAL! IL + J,7,I) =T0TAL!IL*7,1 )+XDLl*J 
TOTAL! I - fj,7* 2)=T0TAL{ I L , 7, 2>FXDL2*J 
32 CONTINUE 
GO TO 30 
31 ! L= I LAST C 

GO TO 30 

22 DO 24 1= IFIR ST, ILA ST 

DO 24 I N= I H • 1 V 

IF! IN.EQ.2) GC TO 27 
DO 20 J«l, 2 

C MOTHS EMERGING FRCM CORN LAY 3 TIMES A3 MANY EGGS AS MOTHS EMERGING 
C FROM SORGHUM. 

XEG G = 1 .0 

IF! J.EQ. 2) XEGG = .333 

TOTAL! I, 7,IN)=( ADTSMG! I, J ,1 N ) APCTCR 3 ( J ) /PC TCRP! KODE < 11 )) )*XEGG + TO T 
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66 

87 

88 
Q9 

90 

91 
93 

93 

94 

95 

96 

97 


27 

26 


2 ALU *7 ,1 N) 

28 CONTINU- 
GO TO 24 

££l£?? T * ILASrc) GO TO 24 

READ ( 1 1 , 503)K *CMALE »C FM4 LE 
F0RMAT(I6U2X # 2F6U L£ 
IFCK.NE* I > GO TO 26 

CONTI^r ,IN,aCMALE * a * 
CONTINUE 
RETURN 
END 


503 


24 


1 " 
1 1 
l ? 
1 3 
1 4 

1 5 
1 6 
! 7 
1 8 
10 
20 
21 
22 

2 3 

24 
2 5 
£6 
27 

25 
29 
3c 
31 
.32 
31 

34 

35 
15 

37 

38 

39 

40 

41 
4 2 

43 

44 

45 

46 

47 


RLDUCT ION OF CATCH OY LOW temperatures* 


SUBROUTINE FLDPOP 

/ ACRE ) AS OFSCP I RED J RY 'hARTST ACK T £t P AL* A ( fq? ? i T ? I £ Ln pap ULATION< MOTHS 
HOLLINGSWORTH 1T4 the . T?ut A N0 HARTSTACK AND 

NIGHT TEMPERATURE ALSoI Hh LIGHT Tf AP CATCHES APE ADJUSTED FOR LCW 
REAL NT TWIN 

T „ . .T*™ ? " • >0/ . -U'H ,<<M . B , 1 0S60 . / 

TPAP FACTt ^r, FOR 40 WATT iilLACKL I GHT T PAP . 

IF (Knni : ( 9) , FiO « 2 J f ’L TURN 

FAC?0=A.(S6C V/C *** * ( ALr,G( f FD+R ) /B ) - ( FD/( FD+R ) ) ) 

NOT COUMtrnf K,N r ' PDUCES CATCH BY 30 PERCENT if center cylinder is 
F ACTD =FACTO * 1 • 3 C 
ADJUSTMENT HF TRAP CATCH Fnp 
I HH~ I H 

IF CKliDF ( 1 > . FO .7 ) I HH= 2 
on 100 IDY= [FIRST • 1LASTC 
DO ICO I N = I H H i I V 
IF { 1 DY • GC *15 4) GO TO 102 
F = l 36.9-. 16 66 7 * I D Y 
A 1 =31 2 .4-5, 5* JOY 
A 2=- 2 1 . 6 76 J + « 2 0 2 P 1 *1 DY 
A3 = ,252667- ,('01 762* I DY 
TFNP = NITMI N( l DY } 
l p < I I T M ?* LT#r<5 * > 1 L- MP = 55 , 

F t =A 1 + ( A 2 + A 3 * temp ) * Tl r M p 
TRAP T N=FACTD*F/F \ 

IF (T* API N.LT •FACT!.)) TRAP TN=FACTD 
GU Tf) 103 

TRAPIN=FACTD 

CONTI* NUF^ *T t IN) = ( T °T AL( I D Y » 7 * I N ) * T R AP I N ) 

IF (KC10E( 3 ) , FG ,0 ) PE TURN 

AND 

J p X = 1 2 
r n = i 

IF(I,F0.2)IM=io 
CALL SMOOTJ-l 
CC NT I MUF 
r'FTUPN 
FND 


1C? 
1 0 3 
1 r o 

TF 


10 1 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

u 

12 

13 

14 

15 

16 

17 

18 


SUBROUTINE GENLN 

OEvLoPMeNT^VAR IOUS A I T Nl-CT E S^GirApnN? A ^r REQU 1 RED F0R 

FROM MEAN DEVELOPMENT TIME* STAGES ALONG WITH EXPECTED DEVIATIONS 
WTE MAVG (370) * SGENLN < 18572 ) 1 1 0> ILAST * KR UN • D UM1 4( 747) *1 GEN ( 370 , 1 4) . 

dimInIion 5fu)?P<£? > * At<4, ‘ A2 ‘ 4, ’ A3(A > 

DIMENSION F ( 4 ) 

DATA F/I ?0»lf 73^1^731 if* l 3 /’* 2235 *“ 555 /,A 3 /3 ** 00135 **000 95/ 

DATA ItODE/O/ *' 

REMDvi^r fi 1 J^o? Q ;. IK0DE, RETURN 

BE CHANGEdI ARD <IFr{,< R | JN.NE.O) RETURN) IF DEVELOPMENT T IME IS TO 
JF(KRUN.NE.O) RETURN 

IF(K0DE(U).EQ.2.AND.IK0DE.EQ.l ) RETURN 
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19 

20 
21 
22 

23 

24 

25 

26 
27 
23 
29 
JO 

31 

32 
23 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 
40 

51 

52 

53 
U 

55 

56 

57 
>0 
59 

r 0 
n 
n 


l F < K0 DE ( 11 } • EQ # 1 .AND. IKODE.E Q.2 > RETURN 
IKODE=KODt < 1 1 ) 

GO TO ( 1 » 1. 2« 1 ) t IKODE 

2 CONTINUE 

1 DEGDA Y{ 1 * 1 ) “ 7 2. 9 
DEGDAY < 1 . 2) = 69. 

DEGD A Y < 2 .1 ) =1 56.8 
DEGDAY ( 2,2)=:1£6.8 
DE GD AY < 3 ,1 >=216 .5 
DEGDA Y< 3 .21=2 16.5 
DEGDAY (4,1 ) = 323 .1 
DEGD AY ( 4 .2 >=32/3 .1 
DEGDA Y( 5* 1 ) = 1 COO. 0 
OEGD AY < 5 » 2 ) = 1 000.0 

3 CONTINUE 
IN=7*< I H— 1 ) 

EDI CK = ED ( I F I RST-1 > 

DO 100 I=IFI RST,I LAST 

IGENX=0 

DO 1 0 A L=1 ,4 

J = L +IN 

SOEV ( I #L )=0 .0 
DO 1 01 K=I ,365 

X= A1 (L >M A2(L )+A3(L) A VG ( I ) > * TEM4 VG < I ) 

SOE V ( I ,L)=SDEV( I.L)+X*X 

IF(ED< K) -EDI CK, LT. DEGDAY (L, I H) ) GO TO 101 
IX =K — I + 1 
IGENX=IGENX+1 X 
IGEN ( I ♦ J )= I X 

SD E V < I ,L ) =SGRT ( SDE V ( I , L )/ IX ) *F ( L I 
GO TO 104 
101 CONTINUE 
1 04 CONTI NUE 

IGtfNf I 9 I N+5) = < 532 7* 35*TEMAVGU >**<-L .8238 }) + . 6 
IGEN( I . I N + 7 ) = IGENX + I GEN ( I , I N +5 > 

IGEN< I .1 N+6)=28 .2 5- .2 403 7 *TEMAVG ( I > 

ED I CK= ED < I ) 

100 CONTINUE 

IF ( IH.EO .2) RETURN 
IF ( TV .EQ .1 >R E TURN 
DO 202 I =1 .4 

202 P{ I >=DLGDA Y{ I ,2) /DEGDA Y( I >1 ) 

DO 200 K = I F I R ST , I L AST 

DO 203 I =1 ,4 

J=I+7 

IGEN(K,J) = IGEN(K»I )*P(I)+.5 
20 3 CONTINUE 

DO 201 I =5 ,7 
J=I + 7 

IGEN < K . J )= IGE N( K, I ) 

2 01 CONTINUE 
20 0 CONTINUE 
RETURN 
END 


1 SUBROUTINE GRAPH 

2 C THE USE OF THIS SUBROUTINE IS OPTIONAL DEPENDING ON THE NEED FOR A 

3 C GRAPH OF THE OUTPUT AND OR INPUT OF THE MODEL. 

4 COMMON K ODE ( 2 0) .DUM66( 4001) t S( 4) ,DUM67< 740) . IFIR ST. D UM68( 1113) ♦ 

5 A I LAST , DU M69 <6 29 0) , EQEL AD { 37 0 1 2 ) , DUM7 0 ( 1480),T0TAL{ 37 0, 9 * 2) ♦ DUM 7 X 

6 B{ 9272) . I GKODS <12.3),DATX< 37 0> ♦ SG R APH (14) 

7 DIMENSION SCALE<5) 

8 INTEGER VFA(5)/U1H ',»,I4t'.' S»X,!1SM '/ 

9 INTEGER VFB( 5 )/ « < 1 H+ ■ , ■ , 4 X, * ♦ • 't'X. !!•.•) • / 

10 INTEGER DIGIT (127) .DlGITl (70). DIGIT2<57> 

11 EQUI VALENCE (DIG I T( 1 ) . D I G I T1 ( 1 )) i(DIGIT(7i ) .0 IGI T2 (1 ) ) 

,2 OAT A D I GIT 1 /<1 ' , *2 % *3’ • '4 ' , ’5* , »6‘ » <7' U8'» • 9S 1 10' ,M 1 1 

13 X 1 2 1 , 1 1 3* » ' 1 4 » . ' 1 5 f * f 1 6 ' f f 17 1 1 20 # • 1 2 1 • • * 22 • § ■ 2 3 1 t 1 24 • t 9 2 

4 X5* ♦ *26 ■ . *27 • . *2 6* . «29< t * 30* ♦ « 31 • ♦ • 32 • . * 3 3* 34' 35 1 36 1 37 ■ , • 39 

,5 X•♦ f 39 , • , 40^,*41♦♦ , 42 , ,*43•,^44 , ,•45 f f , 46■f , 47', , 40«♦♦49 , |'50 , • , 51 , 

16 X,»52 t f'5 3 i .*54«.»55 # . f 56'.«57',*5a',‘59 , .'60*, i 61*.» 62*, , 63»f f 64*. 

7 X'65' ,»66'« »67‘, *69 ' , * 69 1 , *70 V 

18 DATA D I G I T 2 /» 71 • • ' 72 ' , • 73 • . *7 4 ' . * 75* , * 7 6' . ■ 77 % » 79 • , 1 79 * • * 8 

9 X0 , .«81*,*82 4 »'83*, , 84 f ,"85 l . , 86 , i , 87*. 1 88* . ' 89* .* 90* ,• 91 • ,* 92 » » • 93 

10 X 1 , f 94 • . ' 95 ■ . 1 96 • . • 97* . ■ 9 8 ' . * 99’ . 1 100 * . • 101 • , * 10 2’ . • 103* , • 104 ». • 1 05 

l\ XS'106*, •107* .* 108*. *109* .•110* 1 IP ,« llfi ‘ 113* ,<1 14' , 1 115* , • U6 

!2 X • . M 17* . •! 18* , • 1 1 9 • » 1 120 * , ■ 1 21 1 . M 22 • . • 123* , • 124 ' , ■ 1 25* » ■ 126' . * 1 27 

>3 X'/ 

!4 DIMENSION IN SEC T( 3 > . I ST4 GE< 1 0) 

!5 DATA INSECT/' ~ZEA # . * -V I R « . • - Z+V •/ 

!6 DATA I STAGE / "EGGS* * 1 SLA R 1 . ■ LL AR f .• PUPA* , • POAD* ,'OPAD* INPT*t ■ IN 

!7 2 MG • « 1 PAEG* . • E QQ A ' / 

!8 DI ME NS I 0 N SS(2) 


87 



29 

30 

31 

32 

33 

34 

35 

36 

37 

38 
3 3 

40 

41 

42 

43 

44 

45 

46 

47 
43 

49 

50 

51 
£2 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 
67 
63 

69 

70 

71 

72 

73 

74 

75 
75 

77 

78 

79 

80 
ai 

32 

83 

84 

85 

86 
87 
83 

89 

90 

91 

92 

93 

94 

95 

96 

97 
93 
99 

l 00 
101 
1 02 
1 03 
104 
1 05 
1 06 
1 07 
103 
109 
1 10 
1 11 
112 
1 13 


SXZ=0. 
sx v=o • 
sxz v=o. 

N-i 

HR I TE ( 1 *61 1 ) 

61 1 FORMA T ( MM 
N5C ALE =0 
ISfG=l 
KSTG=4 

24 DO 2020 I S- I S TG »KST G 
DO 2020 1N=1 *3 

IF< IGKODEC IS* IN 1.EG.0 IGO TO 2020 
IF ( 1 N« £ Q « 1 ) 5 X Z = 1 
IF( IN*EQ*2)SXV«1 
IF( IN.EQ *3)SX ZV~1 

WRITE< 1*610) N.I STAGEC IS) • INSECT (IN) 
610 FORM AT ( • *,20X* II* * = '.A4«A4> 

N-N*l 

2020 CONTINUE 

IF( IS *GT «6 ) GO TO 25 
I S TG = 9 
KSTG=9 
GO TO 24 

25 CONTINUE 
NN = N 


iriiw.tU* UbU 1U d U2I 

IF (KOOE ( 16) . GE. 100) GO TO 2024 
WRITEU *600) 

600 FORMA T ( • LOS SCALE PLOT 1 ) 

WRITE( 1.604) 

O04 V23^;c;00f5M 00 ,0,,25X ’ ,l000 - 0,<25X * , 10000 - 0 *-23X..10O00< 
GO TO 2021 
2 024 WRI TE( 1, 60 3) 

603 FORMAT ( * LINEAR SCALE PLOT') 

iF(!«;:io l .nSo*fo V u 0 * 1 * AND,sxzv * N -* I,GCI TQ 10 

0,1 >SX1=S( 1+NSCALE) 

GO ( T0^I|°* l> SX1=St 2+ NSCALE) 
to SX 1=S ( 1 + NS CAL E ) 

Go'fo 2 ^ SCALE> ,GT * S( 1+NSCALE ) ) SX1 =S( 2f NSCALE ) 

1 2 IcAClJ \ 1 = 0 S “ LE > +S < 2+ »SC A LE) 

DO 101 1=2,5 

lot SiirlL n =SCALE < I-l)tSXl/4. 

WRITE ( 1 i 606 ) ( $C AL E ( I ) . I “ 1 e i 

606 ^r T iV '^ToCilFilloTFilloM 

IF (NSCA_ E.EO . C) SS{ 1 1 = sxl / 1?7 ' 

IF (NSCALE. EQ.21SSI 21=5X1/ 127* 

2021 &S=i ALE * EQ ‘ 2,G ° ™ l023 127 * 

KSTG = 8 

27 DO till \m S ,l 6 ' K * 1G 

IF ( IN^EQ^l )SXZ='l* E GO TO 2022 

2£ < I N# EQ.2 )SX V-l 
IF ( IN ,EQ . 3) SXZV-I 

WR f ITfcn, 61 o)N, I S T AGE(IS). I NSECT( IN) 

2022 CONTINUE 

ISTG = 1*0 T ,G ° T0 26 

kstg=io 

GO TG 27 
26 CONTINUE 

IF ( N*EQ* NN ) GO JO 2023 

«■ «»■ 

' RITE » '6051 

N GH 0 lki’ 5X,,1,0,,a7X,,10 ‘ 0, ,^7X ’ MOO * 0 ' ,25X ’ , 1000 • 0 ** 25X • 

2025 WRITE! 1 , 603) 

NSCALE=2 
„„ GO TO 13 

2023 55 c 58gf c 5£.VtS5«?iSr TO 3000 

.d L " T , , 

0^3002 ,-».!* ' ' 30 X , • •'•31x*'.',3ox, 31X, '•') 

DO 3002 I N=1 ,3 


0 * 
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114 
1 15 
116 
117 
1 Id 
1 19 
1 20 
121 
122 


IFdGKODE! I. IN) *EQ*0)G0 TO 3002 
IF ( I N ®EQ *3 ) GO TO 3003 

GO TO ( 4 .4.4 .4. 5.5 *5*5*4 , 6 * 3 004 *3 0 05 ) * I 
K^TOT AL< Jo I# I N )/SS ( 2 ) +1 * 

GO TO 302 

K=TOT AL ( J* I* IN) /SSI 1)41. 

GO TO 302 

K-EQELAD( JoIN)/SS( 2) + 1 « 

GO TO 302 


123 

3005 

K=DATX( JJ/SSC 2)41 . 

124 


GO TO 30 2 

1 25 

3 004 

K*DATX< J)/SS< 1)4-1 • 

126 


GO TO 302 

127 

3003 

GO TO <7 .7.7. 7. 6. 8.0. 8.7. 9. 3004» 3005) 1 1 

126 

7 

K*( TOTAL ! J#I . 1 ) 4T0TAL! J * I «2 ) )/SS< l ) >1 . 

129 


GO TO 302 

130 

6 

K=! TOTAL! J.I .1 I4T0TAL! J. 1.2) )/SS<2) *1 * 

131 


GO TO 302 

132 

9 

K=MEQELAD< J. 1 > + EQELAD( J. 2>)/SS(2J + l« 

133 

3 02 

IF ( K* GT* 12 7) K=127 

134 


VFB! 3>=DIGIT< K) 

135 


WRITEU . VFB)N 

136 


N=N+1 

137 

3002 

CONTINUE 

136 

3001 

CONTINUE 

139 


GO TO 5000 

140 

3000 

DO *000 J=IFIRST. ILAST 

141 


WRI T El 1 .601) J 

142 


N=a 

143 


DO 4002 I“1 .10 

144 


DO 4002 IN-1.3 

145 


IF ( I GKOD E< I. IN) *EQ *0) GO TO 4002 

146 


IF { I N • EQ* 3 ) GO TO 4003 

147 


GO TO ( 14. 14. 14.14. 15«15. 15. 15.14916 ). 

148 

14 

B= TOTAL! J.I »I Nl/100* 

149 


GO TO 207 

150 

15 

B=T OT AL ( J. I. IN) 

151 


GO TO 207 

152 

16 

B=EQELAD( J. IN > 

153 


GO TO 207 

154 

4003 

GO TO ! 1 7*17 .17 .1 7 .18 .18.18. 1 B .1 7 .1 9 ) . I 

155 

17 

B~ (TOTAL! J. I. 1)+T0TAL( J.I .2) )/100* 

156 


GO TO 207 

157 

10 

B=< TOTAL (J.I , 1>+T0TAL< J.I .2) ) 

158 


GO TO 207 

159 

19 

B=EQELA9 ( J . 1 ) +E QE LAD! J.2) 

160 

207 

IFCB.LT • 1 *1 >G0 TO 203 

161 


K= ( AL0G1 0!B)/4* 1*127 

162 


IF! K *GT • 127) K — 1 27 

163 

205 

VFB(3 J~0 IGIT (K> 

l€ 4 


WRITE(1 . VFB) N 

165 


GO TO 20 8 

166 

2 03 

K=1 

167 


GO TO 205 

168 

208 

CONTINUE 

t 69 


N=N4l 

170 


GO TO 40 02 

171 

4002 

CONTINUE 

172 

4000 

CONTINUE 

1 73 

5000 

CONTINUE 

174 


RETURN 

175 


END 


1 

2 

3 

4 

5 

6 
7 

a 

9 

10 

11 

12 

13 

14 

15 

16 
17 


SUBROUTINE MIGRAT 

C THIS SUBROUTINE CALCULATES THE PROBABILITY OF MOTHS MIGRATING 
C FROM A PARTICULAR CROP* 

COMMON K CD E ! 2 0 ) .DUM55C27698 ) .PRMGRT! 370 »2 ) . DUM56! 747 ) .KCRPI 
A 1 1 19 ) ♦ IN t K • DUM58* X. SMI GR A < 42 1) 

DIMENSION PCTMGI2) 

DATA PC T MG <11/* 6/ • PCT MG ( 2 ) / • 6 / 

DATA XM/ 1 • 10/* X D/ 1 *45/ 

IF! KODE( 15)*EQ*0) RETURN 
PRMGRTCK •IN)=0.0 
IF!X*LE*XM) RETU RN 
IFC X*GE* XD ) GO TO 1 

PRMGRTCK « IN Is (( X-XM)/( XD-XM) >*PCTMG( IN) 

RETURN 

1 PRMGRT < K*IN)sPCTMGCIN) 

RETURN 

END 


AWAY 

*DUM57< 
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1 

? 

7 

4 

5 

7 

3 

9 

10 
1 1 
12 
1 3 
1 4 
1 5 
1 6 
1 7 
1 8 
1 9 
2 r 
21 
22 
23 


r 

c 

r 


SUBROUTINE MOON 

FOR EACH SIMULATION DAY THIS SUBROUTINE CALCULATES THE NUMPFP OF 
OAYS FPOM FULL MOQN r WHICH IS USED TO SELECT THE PPDPtR MOONLIGHT 
ADJUSTMENT * ACT HP FOP OV I POST I ON * 

COMMON K ODE (20 ) . DUM3 1 ( 4749) , AOMOONi 370 ) , 0UM32< 37 0)iT CLOUD ( 370 ) 
ZS MOON ( 2 4 ,971 ) 

DIMENSION F M ( 29 ) 

DATA FM/1 ,9 1 • , . 9 9 , . 99 , • 99 , • 9 7 * . 0 4 * ♦ 9 1 * * B 6 » . 8 2 , *73, • 6 ] * • 4b * i?4i 
2*53, , 75, .99, ,96, 10*1 , / 

A = K0Or: ( 3 ) 

IFfA.GT. IS. ) A =A — 2 9,5 
A= A - 1 5. 

DO IC O T =1 ,365 
B=T 

A 

IF ( J ,LT. 30 ) GOTO A 

A=A+9Q .5 

J=1 


A CONTINUE 

ADMOON ( I>=FM(J)f(l,-FM(J>) +TCLOUD ( I ) 

100 continue 

RETURN 


FNO 


? 4 , 


3 

4 

5 

6 
7 
9 
9 

10 
1 1 
1 2 
1 3 

14 

15 

16 
17 
16 

19 

20 


23 

24 

25 

26 
27 
26 

29 

30 

31 

32 

33 

34 

35 

36 

37 
36 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 
56 

59 

60 


SUBROUTINE PAR TNP 

C THIS SUFI PHUT I Nfc: CALCULATES PARASITISM UP EGGS AND LARVAE AS DESCRIBED 
C. IN TFXT 

CCLMCN KOOK ( ?0 ) ,DUM21 (4 745), I P I \< S T , DUM 2 2 ( 3 7 3 ) ,SAPEA( 370 ) ,OUM?3(;j7C 
A ) , T L AS T , K a J N , DU M 2 4 ( 2 C 3 5 7 ) , P A P 3 1 T ( 3 7 G * 2 ) , DUM 2 b ( l 4 8 0 ) , I A , 1 M , 5 P A R ( 2 29 
30 ) 

D I MENS inN P A< ? ) 

IF (KPIJN, EQ, C ) GO TO 1C 
Ofi 11 I = IFIRST, ILAST 
DO 11 J- 1 ,2 
1 1 PAPS! T{ I ♦ J) =0 .0 
1 0 l. = KQDL (5) 

C IF KODF 5 = 4 PEAL) PA { 1 ) = C 

GO TQ (1,2,301,1 
1 PMTNl = I A 
PM 4X1 = I A4*°C 


PM IN 2 = T A + 2 p 
pwAX?= I M — 20 
PDDWN= I M 
PFNO= I MF25 

PEAD(5,S00) P A ( 1 ), P A ( 2 ) 

SCO FORMAT (2F6, 3) 

00 4 1= I FIRST * I LAST 

1 F ( I • L I: , PM I N 1 • (J f-> ■ I # G F. . P F. NU ) GG T ij 
f F ( K O 0 L ( 1 l ) • L O • 3 ) GO TG 44 

f f ( I ,ge<pma xn r,n to 7 

FX = ( ( t — P M I N | )/( PMAX1-PMINI > ) * # 5 
GO TO 6 

7 IF ( T ,GT .PMIN2 ) GO TO 8 
FX = ,S 
GO TO 6 

3 IF( I .GF.PMAX? ) GO TO 9 

FX= ,5+ ( ( I -PM I M2 ) /( PM A X 2 — PM I N 2 ) ) *. 5 
GO TO f, 

9 IF ( I ,GT • POOWN ) GO TO 9 9 
c 8 F X = 1 ,0 
GO TO 6 

99 / x=( PENO-I )/( PL ND-F'DOWN ) 

GG r n '6 

44 ! F( I , GF , OMAX 1 ) GG TO 45 

FX = < I-PMIN 1 )/ (PMAX1-PM I N1 ) 

GO TO 6 

4 5 IF { I , 6 T .PMAX? ) GO TO 46 
GO TO 98 


4 6 
6 

4 


f X^=( PFNO-1 ) /{ PEND-PMAX2 ) 
PARS! T( I , 1 )=PA( J )*FX 
PA9SI T ( I ,2) =PA( 2 ) *FX 
CONTINUE 

I F ( L . F 0 , 4 ) G n TO 12 
GO TO 6 


4 


I *?J 5C| ‘ END ~ 5 > 1 * TPICH , OTHERS, OLAPP A 
59 1 FORMAT ( I 6 , 3 F 6 * 0 ) 

IF (I •FOtC) GO TO 5 

c egg par as i t r zm 

N = 4576 9 


• 69 3/ ( N* SAP E A ( ! ) ) 

PT R I CM = 1 , — (FXP(TRT CH*f3 ) ) 
N=500 0 

B=-,6*93/(N*SARFA ( I ) ) 
FOTHPS=l 4 - ( KXP ( OTHERS *13 ) ) 



61 

62 

63 

64 

65 

66 
67 
6R 
60 
7C 
7 1 

72 

73 

74 

75 

76 
7 7 
7H 
70 


f- A PS I T { I ♦ l ) = 1 < (1 .~PfR IC H) * < 1 .-PH THRP > ) 

r ] -3 I NS T A P I APVAL n APAS IT 1 7 M 

M * - • 6 9 3 / ( N * S A U K A { T ) ) 
r HTH Rr» = i * - ( f x p < n i a o p a * i * ) ) 

PAOS I T ( I t 2 ) -PPl MPS 
GO TO P 

i pi: ao < i:j, sc 2 . lnd= c: ) 1 , x, y 

5'' 3 HIRM AT ( 1 6t2f* 6. ) 

IP ( I .FO.O) GO TO G 
PAPS IT( [ ,1 )=X 
PAPS I T ( 1 , ? ) = Y 
GP TO 3 

1 ? f 1 'P AO < j 3 , SO 1 * R N1 ' = S ) I 1 T k 1 CH 
].= ( I .1:0." ) G 1 j TO r - 
PAPS T T { I , 1 > = TR TCH 
r,n to \p 
5 PE TUP N 
fc NO 


t 

5> 

3 

4 

5 

6 

7 

8 


9 

10 

11 

12 

13 

14 

1 5 
16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 

2 B 

29 

30 
3! 

32 

33 


34 

15 

36 

37 
3B 

39 

40 

41 

42 

43 

44 

45 
4 6 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 
50 


59 

60 


SUBROUTINE PPDINP 

C THIS SUBROUTINE CALCULATES PREDATION OF FGGS AND LARVAE A$ DESCRIBED 
C *IN TE XT* 

COMMON K 0 DE.( 20 > , DUM 26 ( 47 A 5 > , I F I R 5 T , DU M 2 7 < 37 3 > , SAREA< 370 ) ,DUM 2 8(370 
A ) , I LAST , KRJN, DUM 29 ( 21097 ) , PPE DO P ( 3 70 , 2 ) * DUM 30 ( 740 ) » I A , I M, SPRDt 229 0 
B ) 

OIMFNSIGN P A < 2 ) 

nr mens icn pdfrte<7 ) , pdlrte < r > ,preo<7 > 

DATA PDERTF /I * 25 , t . 5 ,* 25 ,. 5 , 1 * 0 , * 75 , 1 . 75 / 

DATA PDLRTE/. 5, 1 . 0 , . 25 , . 75,1 * 25,1 . 75 , 1,5 0 / 

C POL RTF = PR E CAT OR RATING FOR LARVAE 
C PDE 9 T£=PPEDAT 0 R RATING FOR EGGS 

C 3 =SEYMINUS , 2 = 0 THER COCCI NELL IDAF , 3 “ OT HEP COL EC PT ERA , 

C 4 — ORT US , 5 =GEOCQR I S , 6 - OT HEP HEM I PTFR A ♦ 7 = CHRY 50 P A IMMATURES 
IFOCRUN.EO.C ) GO TO 10 
DO 11 I=IFIRST, ILAST 
DO U J=> 1,2 
1 1 PREDDR < I , J )=0 .0 
1 0 L-KODE ( 12 ) 

FX= 1 . 0 

GO TO (l , 2 , 3,30 ) * L 

1 FMIN!=IA 
FMAXl=tA +20 
PMIN 2 =I A +20 
PMA X 2 s I M -20 
FO OWN = I M 
PFND=IM +25 

READ( 5 , 500 ) PA(l)iPA( 2 ) 

500 FORMAT ( 3 F 6 * 3 ) 

DO 4 IsslFIRST, ILAST 

IF( I .LF.PMIN 1 * 0 R* I *GE .PEND) GO TO 4 

IF(KODE( 11 ) * E 0 • 3 ) GO TO 44 

IF( I .GE.PMAX1 ) GO TO 7 

FX = < < I-PMIN 1 ) /( PMA XI -PM INI ) ) * ,5 

GO TO 6 

7 IF(I .GT.PMIN 2 ) GO TO 8 
F X = * 5 
, GO TO 6 

4 I F ( T * GE , PMA X 2 ) GO TO 9 

FX~ • 5 + ( (I~PMIN 2 )/< PMAX 2 -PMIN 2 ) )*,5 
GO TO 6 

9 IF( I *GT .PDOWN) GO TO 99 

98 F X = t *0 
GO TO 6 

99 FX“< PEND- I ) /( PEND-* PD OWN) 

GO TO 6 

44 I F ( I • GF * FMA X 1 ) GO TO 45 
FX=( I-PMIN 1 )/( PMAXl-PMINl ) 

GO TO 6 

45 IF( I*GT • PMAX 2 ) GO TO 46 

go rn 

46 FX=<PEND-I )/<PEND-PMAX 2 > 

6 PRFDOR < I , 1 ) =PA < 1 )*FX 

PREOOR ( I , 2 ) =PA ( 2 )*FX 
4 CONTINUE 
GO TO 5 

2 PE ADC 1 4 ,501 ,FND~ 5 ) I , ( PRED ( J ) , J= 1 , 7 > , SPID,CHRYSP 

501 FORMAT ( 9 X, I 3 , 6 F 6 . 0 , 6 X, 2 F 6 . 0 , 6 X,F 6 # 0 ) 

I F < I • E Q * 0 ) GO TO 5 

(ITHEGG - 0 • 0 


41 



CTHLAP=0 .0 
DO HO K=1 ,7 

OTHEGG=QTHFGG+PRED ( K ) *PDERTE (K) 

GTHLAP=OTHLAR + PRED ( K ) +PDLPTF (K > 

20 CONTINUE 
C EGG PREDATION 
N- 44 R2 7 

G=-.693/( M*SAREA< I ) ) 

PCHRYS-1 .-(EXP< C HRYS P*B ) ) 

N = 1 436 8 

B=-,693/<N*SAREA < I ) ) 

POTHRS^? 4-(FXP( 0 T HEGG * B ) ) 

PR- DOR ( I *1 >=1 .-< J .-PCHRYS>*(1 . — PPTHRS ) 

C i-3 INSTAP LARVAE PREDATION 
M= 1 3051 

B-- * 693/ ( N *SAPFA ( I ) ) 

PCHRYS= 1 .-(EXP(CHRYSP*B> ) 

U-P. 929 

e-- 1 693/(N*5AREA ( I ) ) 

POTHPS -1 ( E X P { OTHL A R * B ) ) 

PR~DOR ( I 1 2) - 1 •— ( 1 .-PCHRYS)*( 1 • — POTHPS ) 

GO TO 2 

-1 PEAD(14,50?.END=5) i*x*y*z 
502 FORMAT { I6.3E6. 3 ) 

IF< I «EQ • 0 ) GO TO 5 
PFEDOR ( I * i )=X 
PREOOR ( I,2)=Y 
GO TO 3 
30 I F p -365 
TLL-t 

3T PE AD U 4 « 50 3 « E Np =3 A ) SP » ALB . HE MF • CHL » C HA , CALL OP, I 
E03 fqpmAT (6X*6F5 .0.13) 
tF(I.FQ#0) GO TO 34 

PREDOR ( I f \ )=ALB*1 .5+HFMP*! .0+CHL*I . 75 +CH A+C ALL OP 
PREDORd ,2>=ALB+HEMP*1 .25+CHL*l • 5+C HA +C ALLO P 
I F ( I F F . GT .1 ) IFF = | 

IL=ILL 

IF ( ILL iLT*I > I LL= I 

I F { IFF. FO. ILL >GO TO 32 

IC'D=ILL-TL 

TF( IDD .EQ.t ) G 0 to 32 
XDU-PRFDOR ( ILL . 1 )-PPEDOR( IL , t > 

XDL? = PPEDOR( ILL t?)-PREDOR< IL,2) 

XDU = XDLt/It.O 
XOL2-XDL2/I DD 
inD=IDD-l 
DO 33 t I DD 

PREDOP ( IL + J.J )=PPLOOR< I L * I )+XDLl*J 
PRFDOR { IL+J . 2 )=PRFDOP( 1 Li2)+XDL2*J 

33 CONTINUE 
GO TO 31 

32 TL-ILL 
GO TO 31 

34 NF =14 36 3 
NL =89 29 

DO 35 J=IFF, ILL 
S=-.693/<MFYSARFA( J > ) 

PREDOP( J.1 )=! FXP< PPEDOR( J * 1 ) *B) ) 

P = -.eP3 /{NL ySAREA( J > ) 

PRFDGR( J.2)=l .-{ FX P( PREDOP ( J , 2)*R) ) 

35 CONTINUE 
5 RETURN 


23 


END 


1 SUBROUTINE SMOOTH 

2 C THIS SUBROUTINE WHEN CALLED APPLIES A DOUBLE 3-DAY MOVING AVERAGE 

3 C TO INSECT POPULATION DATA. 

4 COMMON DUM62< 4765) *1 XF *1 XL* DUM63 (7411 ) *SM<370.24) ,DU M64 ( 9 267 ) • IN. 

5 A DU M65 (3) * IPX* SS MOO (420) 

6 C IXF-FIRST DAY OF SMOOTH IXL-LAST DAY OP SMOOTH IP- STAGE IN- INSECT 

T DIMENSION S( 5 ) 

3 IP^IPX+IN 

? DO 100 I «3*5 

) 100 S(I)-0.0 

l IXXF-IXF+2 

> S( 1 )-SM ( IXF+1 .IP) 

3 S(2 ) -5M( IXF * IP) 

% DO 101 I«I XXF *1 XL 

> K=0 

5 DO 102 J-2 *5 

r S< 5-K)~S(4-K) 



IS 

102 

K-K+l 

19 


SC 1)=SM( I .IP) 

20 


SM i 1-2* IP )=<S { 5t+2»*S( 4 ) + 3* *S{ 3)+2.*S< 2) + S< I ) ) 79. 

21 

101 

CONTINUE 

22 


SM i I XL- l*IP) = (S(4)+2. *50)4-3 «*S(2) +2 o*S(l))/8* 

23 


SM < I XL . I P )« ( S (3 ) + 2 **S<2 ) 4-3 * * S( 1 ) ) / 6* 

24 


RETURN 

25 


END 


I 



SUBROUTINE SPRAY 

? 

r 

READS IN T NSECT IC IRF MOP T A L I T I ES - J- 1 THRU l 2 

3 

c 

THIS SUBROUTINE HANDLES ALL INPUT OF MORTALITIES CAUSED RY 

4 

r 

I NSEC TIC FOES* 

5 



COMMON KODF ( HO ) • DUM33 ( 4 745 ) , I F I P S T , D UM 3 4 , I H . T V , 0ljM35( 1 1 10 ) . I LAST 

A 



AKRUN, 0UM3M l 5 1 77 ) , PdS { 3 70 • 7* 2 ) * SSPR A Y( 4512 ) 

7 



D I ME NS TON p ( 7 ♦? ) 

rt 



I F ( K 00 E ( 4 ) • EO « ? I RETURN 

9 



IF (KODE ( 4 ) • EO *0 » ANR.KRUN.EQ* 0 ) RETURN 

10 



nu 2 j=i ,7 

1 1 



on 2 T- IF IRST * I LAST 

1 2 



DO 2 IM-IHtIV 

13 



POS( I. J,TM)=O.C 

14 



IF (KCJDF ( 4 ) . FO .0 ) RETURN 

1*5 

1 CO 

R E AO < 1 2 . 50 0 . END = 5 ) I . ( ( P ( J . IN ) ,J = 1 , 7) , TN^IH, I V ) 

16 


rc o 

FORMAT ( 15, 1 4F4«3 ) 

17 

* 


IF ( t ,FQ ,C > RETURN 

1* 



DO 10 1 rM=IH, IV 

19 



nn 101 J-l . 7 



1C 1 

POS( I , J , TN) =P< J , IN) 

21 



GO TO 1x00 

?? 

5 


RF TURN 

23 



END 


1 


subroutine temp 

■o 

C ALL OF THE TEMPERATURE AND CLOUD COVER TNPIJT AND THE CALCULATION OF 

3 

C DAYS ARF HANDLED BY TEMP* 

4 


REAL NITMTN 

5 


COMMON KHDE ( ? 0 ) , DU M5 ( 1 4 8 1 ) . N I TM I N ( 3 7 0 ) * FD ( 3 70 ) . F£OCOPN< 30 0) .DUMft 

6 


A ( 1 854) . EGL YTM( 370) » I F IP ST • !.HJM7( 743 ) , T CLOUD ( 37 0 3 * I L AST . DUMA ( 5 9 28 ) 

7 


OTCMAVGC 370) .STEMP( 185 72 )» TFT STT * I LA STT 

n- 


DIMENSION CL OUR < 24 ) 

9 


DIMENSION TMFA(BOO) 

l 0 


0 I MENS T ON TM ( 24 ) 

1 1 


KKK=K OOF ( 2 ) 

12 


KK-KOni-( 10) 

13 


FO ( IFISTT-1 ) = o * c 

14 


FOCnRNf TFISTT-1 )-g*o 

15 


FT = 54 ♦ 7 

K> 


T 1 -5 0 « 

17 


T2~86 . 

1 8 


TQ~T2-T 1 

19 

50 3 

FORMAT < l 3, 2F3 *C ) 



DO ?04 IRY- IF I STT* TLA STT 

21 


TF(KKK.NF*3) GO TO 2 C 0 

22 


1 F ( I 0 Y « FO • IF ISTT ) GO TO 404 

23 


GO TO 410 

24 

BCD 

K PUN" 0 

PS 

POO 

PE AD ( 1 f> ,50? ,FNR-402 ) !.(TM(J).J=1.24) 

26 

EC 

2 FORM AT ( * T 3 * 3X , 24F3 . 0 ) 

27 


IF (I .NF. IDV) GO TO 200 

28 


M~ I 

20 


GO TO 490 

30 

4 02 

CONTINUE 

3 1 


WPT TF ( 6, 1 00 0 ) I F TSTT, M 

32 

l ccc 

FORMAT ( 1 1 1 ♦ • HOURLY JFMPFRATUPF INPUT FROM DAYS TD M3) 

33 

4 04 

RE AD ( 16, 50 3 ) I , THP.TL.P 

34 


IF( I .NF.inY) GO TO 4C4 

35 


WRITE(6«?09f! ) IDY.ILASTT ^ ^ . T „ , __ . » 

35 

20CO 

FORMAT ( 1 1 . 1 MA X / M 1 N TFMPIRATURF INPUT FROM DAYS M3, 1 TO **I3/) 

37 


KKK = 3 

3 A 

4 1 0 

CONTINUE 

39 


THO=THP 

40 


TLO-TLP 

41 


READ ( 1 6, 503 ) I *THP,TLP 

42 


TAfl= ( THO + TLD ) /?. * 

43 


TPO= ( THO-TLO) /? * 

44 


TAM-(TH04*TLP)/2 * 

45 


TPM=( THO-TL P ) /? * 

4 6 


TM(1 )=TAO+TPG*( -*96593) 

47 


DO 4 1 7 J=2, 1 4 


43 



46 


49 

4 17 

50 


51 


52 

4 1 9 

53 

4 90 

54 

20 i 

55 

504 

56 


57 


58 

2 06 

59 


60 

30 5 

61 

21 3 

62 


63 


64 


65 


66 


67 


68 


69 


70 

585 

71 


72 

35 

73 


74 


75 

3? 

76 

33 

77 

34 

7fl 

30 

70 


BO 


ai 


BP 

56 3 

8 3 


84 

561 

85 

56? 

86 

3 

87 


88 


80 


90 


91 


92 

3C4 

93 

20 ? 

94 



TM( J ) = T An+TPn*COS( . 26 1 0 *( 1 4- J ) ) 

CP NT I NUF 

nn 4 \ 9 Jr: 15,? 4 

TM( J)=TA M+TPM ❖CO S ( . 26 ! 8 *( J- 1 4 ) ) 

CONTINUE 

TF(KK.EO.O) GO TO 105 

RF ADM 7 ,50 4 , END =30 6 ) I , ( CLOUD (J ) , J= 1 , 24 ) 
FORMAT ( 2X , I 3, 3X , 24F3, 1 ) 

I F ( I « NF • I O Y ) GO TO 20 1 
GO TO 30 5 
PONT I NUE 

KK=0 

CONTINUE 

ED ( IDY )=E0( IDY - 1 ) 

ED CORN ( I DY ) =EOCORN< I DY- 1 ) 

T F M A V G ( IDY)~0.0 
NlTMINf T OY ) =9 • 9 


no 1 1=1 , ?4 

TF (TM< I ) .GT.Q2 ) GO TO 565 
T C M“(TM( T ) - FT ) / 2 4 « 

TF ( T F M * L T ■ 0 • ) TFM = 0 • 

GO TO 35 

TFMOVP=TM( T ) -92 . 

TFM={ 3 7,3-TEMOVR )/2 4 , 

CONTI NUE 


X = TM( T )-T 1 
IF( X )39 ,39, 3? 

IF(X-TD) 34.34 ,33 
X = TO- 1 . 423* (X-TD ) 

FDCORN( I Oy ) =EDCOPN ( TOY ) 4-X/24 * 

CONTTNUF 

TF(t»GT # 6) GO TO 563 

NTTMIN ( I DY- 1 ) = N I T M I N ( I D Y— 1 ) + TM( I )/l2. 
GO TO 561 

I F < T , L T *19 ) GO TO 561 

NlTMTN (I DY) =NI TM IN( TOY ) 4-TM( I ) /I 2, 

E0( IDY )=ED( I D Y ) +TEM 

TFMAVGI IDY) =TEMAVG( IDY ) fTM( I )/24, 
CONTINUE 


OA YLN= 12* 14 + 1 ♦ 94* COS ( • 0 l 72 1 * ( TOY- 1 75 ) ) 

TFGG-1 3 ■ 5 + D A YL.N/2 • 

FGLYTMf TDY) = (TM( IEGG)+TM( IFGG + 1 ) +TM ( T EGG 4*2) )/3« 

I F f K K tFQ.O) GO TO 304 

T CLOUD ( I Dy ) - ( CLOUO( TEGG ) +CLOUO( IFGG+1 ) 4- CLOUD ( IFGG+2 ) )/3 . 

CONTINUE 

RFTUPN 

END 


1 

2 

3 

4 

5 

6 
7 
5 
Q 

10 
1 1 
12 
1 3 
1 4 
t 5 
16 
1 7 
IB 
1 9 
20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


C 

C 


C 

C 


SUBROUTINE TRKf.PM 

THIS SUM/? OUT INF PROVIDES A DYNAMIC POPULATION MODEL 
FOR TPICHOGPAMMA . 

COMMON KOUF (2l ) . DUMA C ( 51 1 9) , S AREA ( 3 70 ) , DUMA 1 ( 370 > , I L AST , DUM4 2 (7 46 ) 

4* 1 r, n K i n M ?. ! n * ' '« U i 4 f l . 3 5 rY-S 4 ) \ 1 0 T AL ( i 7 1; ’ 9 • 2 ) ' f> 05 ( 3 7 0 » 7 t 2 ) , P A PS I T ( 37 0 , 
y l ?n ’ ' > ><>UM45( > , I ,0UM*6( ) I I N , l)UM47( *1 C ) , TP , !.)FATH ( tj 

( . t 2 ) , S T P K ( 3 ) 


D I vitiNS ION F r W A T F ( 2 C ) 

DATA TKSUPV/1 • / 

LATA F E.P AT { . / • 5 i , F , . 8 7, , 9 6 , ,95,. 95 , « 9 3 , 

? i . 4 , - 3 , , 2 i . 1 , c • / 

IF(TM.LT, 1 )M:T URN 

TOT AL < I , 9, IN ) = TOTAL! I , 9, I N ) +TP 
T 5 = T P 


» 9 , • 8 2 , ,74 , • 66 , #6 


, • 55 i • 5 » • 4 5 


K = T 4-TGENf I , 1 , IN) +IGLNI I , 2 , I N ) - 1 

I F ( K • G T , I L A 5 T ) r>f • T 1 J I j i • J 

11=14-1 

DO ion J = I I ,K 

X— PRFDOP ( J ,1)^(1 ( J , 2, IN ) ) 

X= U • “X ) 1 ♦ - DEATH ( 1 * IN ) ) *( l , -PCJS ( J , 7 , I N) ) 


TOTAL ( J, O, I N)=T OTAL { J,5 , I N) 4-TS 
TOTAL (J , 1 , I N)=TOTAL( J , 1 , I N) 4-TS 
100 CONTINUF 

AS FUNCTION OF SAPr.A 
EACH EGG * 


IF ( K 0 D f ( 5 ) • NE « 4 ) PR TURN 
AS3UMF EGGS A°F HARD TO FIND 
ASSUME 2 FEMALES EMERGE FROM 
T$ = TS*2 .*SAROA( 1 I ) 


K = KFi 
KK=K+ 1 0 

IF(KK ,GT, ILA.ST) KK=ILAST 
jagf:=i 


44 



33 DO 200 J = K » KK 

34 T S = T S * T K 5 iJ ^ V 

35 PARS I T C J ♦ 1 ) =TS*FERATE < I AGE ) + PARS I T ( J , 1) 

36 I AG!£ = t AGF + 1 

37 pro CONTINUF 

3A RETURN 

39 END 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 


FUNCTION! DISTR(DISTS) 

C THIS SUBROUTINE CALCULATES THE DEVELOPMENT DISTRIBUTIONS OF COHERTS 
C OF EGGS 9 LARVAE* AND PUPA* A NORMAL DISTRIBUTION IS ASSUMED (ST INNER 
C 1974)* THE ARRAY GAUS IS AN ACCUMULATION OF THE NORMAL DISTRIBUTION 
C PROBABILITIES FOR + OR - 4 STANDARD DEVIATIONS FROM THE MEAN* 

COMMON DUM72< 30328 ) . X * D UM 73 < 41 8) * XBAR *XSD « DUM74 
DI MENS ION GAUS1 (40 I * G AU5 2 < 4 0 ) * GAUS < 0 I ) 

EQUIVALENCE IGAUSHl) *GAUS<2>) • ( GAUS2 (1 >*GAUS<42>> 

DAT A GAU S< 1 ) / .0 00 0 3/ 

DATA GAU SI / *00005 • *0 000 7* . 0 0 01 * *0 0 0 1 5 • .0 0 0 20 * *00 03* *0 005* *0 007**00 

1 10* * 0 0 1 3 * • 001 9* *0 026** 0035 *• 0047** 006 2 »* 00 0 2 • • 01 07 * • 01 39 • *0 1 79 • #02 
220 » *0287 » *035 9 « *0 4 46* *0 54 8* *066 8* * 08 0 8* . 0 96 8» * 1 1 51 . • 1357* • 1 587* . 1 0 
341 «* 21 19 * *242 0. *2 743 *.3085* *3446 « *3021 * .4207* *4602* *5000^ 

DATA GAUS 2/ *5388* • 5793* * 6 179. . 6554, • 691 5* * 7257 * . 7580 **7881 **8159.* 

2 8413 • .864 3 *. 8 84 9# .903 2* *9192* *9 33 2*. 9 45 2. *9 554* .9 641 * .9713* * 9772 *. 
39621 * .98 61 * . 9893* . 991 8* .9938*. 9953 *• 9965 * .9974 * .9901 * . 9907 * *9990* . 
49993* .9995*. 9997* .9998* .99985* . 9999 * .99993* . 99995 *• 9 9997/ 

Y— < X- XBAR+.5) / X SD 
IF ( Y . GT .—4 * ) GO TO 1 
DI STR-0. 

RETURN 

1 CONTINUE 

IF (r.LT.4. )GO TO 2 
DI STR— 1 • — DIST S 
DISTS^l • 

RE TURN 

2 CONTINUE 

INDEX=IFIX(<Y+4 .1 1410 *1 
Y~GAUS< INDEX) 

DIST R~ Y— DI STS 
DI STS-Y 
RETURN 
END 


EXAMPLE INPUT 


1 

3 1 

27 

0 1 

1 

160 1 

2 

45 4 0 

40 

30 150 

900 

300 45 

3 

82 146 





4 

59 236 





5 

56 

71 

140 


1 4000 

6 

060 

01. 0 




7 

150 

200 




8 

1 1 1 





9 



1 1 1 



10 






It 






12 

82 

10 




13 

83 





14 

P4 

02 




15 

Q 5 





16 

6 6 

08 

02 



17 

37 





18 

88 





19 

89 

02 




20 

90 

02 

04 



21 

91 

1 2 

02 



22 

92 

0 2 




23 

9 3 





24 

94 


01 



25 

95 

02 




26 

96 

06 

06 



27 

9 7 

1 0 

07 



28 

98 

81 

37 



29 

99 

09 

1 0 



30 

1 0 0 

Ot 

04 



31 

l 01 


02 



32 

102 

01 




33 

103 






1 205 11111 102 1 3 

40 40 30 150 900 400 


2 1 8 1 8 59 1 


ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 
ALL CROPS LT 1975 


45 



7 A 

15 

3 6 
7 7 
7.5 
?q 
AO 
A 1 
4? 

4 ? 

44 

45 

46 

4 7 
45 
4<5 
50 

5 i 
5? 

53 

54 
5 C 

56 

57 
5 5 
5Q 
60 
61 
65 

63 

64 

65 

66 
67 
65 

69 

70 

71 

72 
77 

74 

75 

76 

77 
7 8 
79 
BO 
R! 
ft? 
83 
ft 4 
85 
R6 
87 
8P 
BO 

90 

91 

92 
97 

04 

05 

96 

97 

06 
90 

100 

}0J 

10? 

1 03 

104 

105 
1 06 
107 
1 OR 
109 
UO 

m 

112 

113 

114 

115 

116 

117 

118 


105 
I 0 5 
1 0 6 
i r 

1 C n 
i 0 0 
1 1 0 
1 1 1 
1 t 2 

1 ! 7 
1 1 4 

J 1 - T 

1 \ 6 
117 
1 1 7 
i 1 9 
120 
1 2 
1 22 
1 2? 

* 2 
3 25 
12 5 
127 
* 2 Q 
i ? o 

1 30 
1 31 
1 3 > 

1 33 
1 34 
1 35 
3 36 

! 37 

i 7ft 
1 39 
1 40 


^ 8 
0 2 

i 1 


1 0 
;? 5 
» 6 
1 0 
0* 
20 
r o 
.1 B 
75 
2 “* 
1 0 
55 
105 
95 
6 0 
1 0 R 

1 90 

* 7 
9 7 
59 
05 
2 0 

03 
0? 
0 ? 
O'* 
O c 
2 3 
2ft 
31 


r 4 

C c 

03 

5 


OP 

’ 1 

nc 

?ft 

00 

v- 0 


1 4 
50 

07 

0 1 
1 n 


05 
9 5 
Cl 

ro 

3 7 
1 O 



1 4 

\ 

-1 

0 

f. 

0 


1 4 

p 

- 

B 

c 

fl 


14 3 

1 > 

0 

78 


1 44 

65 


3 


1 45 

65 


<r 


1 4 

6 

30 



1 

F ? 

75 

*6 

*4 

5 3 

52 

\ 

( 0 

75 

56 

34 

54 

57 

l 

61 

75 

59 

66 

6 6 

54 

1 

6 2 

7 C 

5ft 

c 7 

56 

56 

1 

6 3 

75 

40 

AC 

40 

40 

1 

64 

">5 

4 1 

At 

a r 

40 

1 

6 5 

75 

5? 

54 

54 

64 

1 

£ 6 

75 

69 

69 

69 

6Q 

1 

67 

75 

47 

46 

45 

a ? 

1 

6 8 

75 

35 

55 

<? } ~ 

5? 

1 

6 9 

75 

64 

ez 

62 

62 

1 

70 

75 

5 3 

54 

56 

55 

1 

71 

7 5 

7? 

72 

7? 

71 

1 

71 

75 

42 

41 

4 2 

40 

1 

7 3 

75 

32 

30 

29 

30 

1 

74 

7ft 

45 

47 

4 6 

44 

1 

7ft 

75 

60 

6C 

5n 

5 7 

1 

76 

75 

5t 

51 

51 

FO 

1 

7^ 

75 

55 

55 

53 

54 

1 

7j 

7« 

47 

45 

an 

4 7 

1 

-r } 

76 

53 

54 

64 

53 

1 

PO 

75 

61 

63 

63 

62 

1 

ft t 

75 

59 

60 

69 

69 

J 

ft 2 

75 

70 

70 

7Q 

7! 

r 

83 

7 e 

56 

54 

5? 

51 

1 

ft 4 S75 

51 

5° 

51 

47 

1 

ft 5 

75 

55 

52 

53 

52 

1 

86 

7? 

.7 0 

70 

70 

70 

1 

ft 7 

7 C 

71 

71 

71 

7! 

t 

ft ft 

75 

45 

45 

4 4 

42 

1 

ft 3 

75 

35 

75 

3F 

35 

1 

90 

“>5 

43 

4 1 

it 

43 

1 

91 

75 

4ft 

46 

46 

48 

1 

o> 

7-5 

62 

62 

6 3 

64 

i 

93 

75 

41 

A! 

39 

39 

1 

94 

75 

46 

45 

44 

4 3 

1 

95 

75 

53 

5 0 

51 

ft 1 

1 

96 

76 

60 

60 

61 

e 1 

1 

97 

75 

60 

59 

5ft 

57 

* 

98 

75 

7Q 

65 

VI 

62 


ALL CPOPS LT 1975 
ALL CPOPS L T 1975 
ALL CHOPS LT 1975 
CORN LT 1 975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1 975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN L T 1975 ACTU 
CORN LT 1 975 ACTU 
CORN LT J 975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN L ■** 1 975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1 975 ACTU 
CORN LT 1975 ACTU 
CORN LT \ 975 ACTU 
CORN LT 1 975 ACTU 
CORN LT 1975 ACTU 
CORN L T 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1 975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT i 975 ACTU 
CORN LT 1 975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1975 ACTU 
CORN LT 1 975 ACTU 
CORN LT 1975 ACTU 


51 4P 

5ft tl 

52 51 
56 *5 
40 40 

4 0 4 1 

5 5 5 7 

67 67 
42 41 

5 0 60 
6? 54 

55 57 

71 70 

39 39 
31 30 

44 4 5 

56 55 

49 ^0 

57 57 
4? 4 2 

5 0 4Q 

6 2 t 2 

6 9 69 

72 7? 

50 50 
4 6 4 7 
52 53 
70 70 

7 1 71 

4 6 1 o 
34 33 
4? 4 5 
47 47 

64 65 

37 34 
43 4 3 
50 50 
eo 60 

57 58 

63 61 


4ft 53 
62 63 
47 52 
54 54 

40 40 

41 42 
5ft 50 
67 68 
41 4 5 
50 50 
50 51 
5ft 56 
70 70 
39 3ft 
30 35 
46 48 

35 57 
50 52 
84 58 
4 3 53 
4 9 56 
63 63 

69 70 

70 70 
50 56 
46 54 
57 59 
7 C 70 
72 72 

4 0 40 
34 ?8 
46 50 
48 50 
66 67 

36 41 
43 47 

5 3 57 
60 61 
59 61 
59 59 


58 60 
65 6 8 
57 61 

54 54 
40 41 

4 3 4 3 
62 65 
12 75 
49 5? 

r, p 2 r> 

5 3 ft ft 

59 61 

70 73 

3 8 3 8 

4 0 43 

5 0 54 
59 63 

55 60 
62 66 

5 9 6 4 

62 7 1 
65 67 
72 74 

71 71 

61 64 

6 4 71 

e 1 63 

71 73 

72 73 

4 0 4 0 
42 46 
52 52 
52 87 
6 0 56 
46 50 

5 3 65 

6 1 62 
6 3 65 

63 66 
60 62 


6? 

72 
65 
54 
4 3 
43 
67 
78 

56 
53 
60 
6 1 
7ft 

40 
48 

57 

64 
61 
59 

69 

72 

70 
75 

73 

65 

74 
65 
74 
73 

41 
40 
5>4 53 

63 
59 
53 
59 
65 
70 
67 

64 


€6 65 67 

74 76 79 

67 69 71 
57 57 59 

45 46 
49 51 

69 70 77 

77 80 74 
5o 61 63 

56 66 58 

63 63 64 

57 65 67 

77 71 
44 4 7 

_ 5? 55 

6 0 61 6 1 

64 66 67 

64 *5 61 

70 72 72 

75 77 7 ft 
74 77 79 

70 72 

78 81 
73 74 7ft 
6 B 70 73 
77 78 80 
66 68 70 
77 8J 83 

60 56 52 
42 43 42 

52 55 55 
- - 54 5 6 

66 70 73 
62 67 61 
57 59 61 

61 63 64 

67 67 67 

71 75 75 

68 72 70 
64 66 67 


43 

44 


81 

4? 

51 


74 

82 


72 74 
7ft 80 

73 73 

59 50 

46 47 
54 54 

77 78 

72 70 
63 64 

62 65 

63 61 
77 78 
66 61 
44 46 
57 58 

60 61 
66 64 
59 58 

72 71 
79 78 

79 78 
76 76 
81 81 

75 75 

73 71 

80 79 
68 68 
84 82 
49 46 
40 39 
56 57 
59 60 

74 75 
62 61 
62 64 

65 65 

68 70 

76 -*4 

69 69 

66 65 


73 71 66 

79 74 67 
71 6ft 63 

57 54 50 

47 46 42 
55 54 53 
77 75 72 
67 63 59 
64 61 

71 71 

60 58 57 

76 74 7 2 

58 55 54 

44 44 40 

59 57 53 

60 60 57 

63 61 59 

58 5ft 57 
70 67 63 

77 73 67 

77 75 70 

74 73 71 

78 77 75 

75 7? 72 

70 67 63 
77 74 69 
67 66 67 
81 78 75 

45 44 44 

38 39 39 

57 55 53 
60 57 53 
74 69 65 

59 57 53 
62 60 56 

64 61 59 

65 64 62 

73 70 68 
69 68 66 
64 64 64 


64 61 
63 63 
61 59 
49 46 
4 3 4 3 
53 53 

70 68 

56 64 
57 57 56 
70 69 68 

57 56 

71 70 
52 52 
39 37 
51 50 

57 58 
56 55 

56 5B 

58 57 
63 61 
68 67 
70 70 
73 72 
67 65 

59 54 

65 63 
70 71 
73 73 
4 5 45 
38 37 
43 44 
50 49 
62 61 

50 48 

51 49 

57 56 
61 60 
67 65 

66 67 
63 6 3 


46 



U 9 

120 

121 

122 

1 

1 

l 

1 
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